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We predict cross sections in deep inelastic scattering (DIS) for the production of two jets — 
one along the proton beam direction created by initial state radiation (ISR) and another 
CO | created by final state radiation after the hard collision. Our results include fixed order 

corrections and a summation of large logarithms up to next-to-next-to-leading logarithmic 
' (NNLL) accuracy in resummed perturbation theory. We make predictions for three versions 

of a DIS event shape 1-jcttiness, each of which constrains hadronic final states to be well 
collimatcd into two jets along the beam and final-state jet directions, but which differ in 
their sensitivity to the transverse momentum of the ISR from the proton beam. We use 
the tools of soft collinear effective theory (SCET) to derive factorization theorems for these 
j-^. three versions of 1-jettiness. The sensitivity to the ISR gives rise to significantly different 

f\j ■ structures in the corresponding factorization theorems — for example, dependence on either 

the ordinary or the generalized fc^-dependent beam function. Despite the differences among 
1-jcttiness definitions, we show that the leading nonperturbative correction that shifts the 
tail region of their distributions is given by a single universal nonperturbative parameter 
fit, even accounting for hadron mass effects. Finally, we give numerical results for Q 2 and x 
Q-r values explored at the HERA collider, emphasizing that the target of our factorization-based 

analyses is to open the door for higher-precision jet phenomenology in DIS. 
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I. INTRODUCTION 

Deep inelastic scattering (DIS) of an energetic lepton 
from a proton target at large momentum transfer probes 
the partonic structure of the proton and the nature of 
the strong interaction, and was an important ingredient 
in the development of the theory of Quantum Chromody- 
namics (QCD) [IHi] . Modern DIS experiments at HERA 
and Jefferson Lab continue to illuminate the internal par- 
tonic structure of hadrons, yielding information on par- 
ton distribution functions of all types, as well as the value 
of the strong coupling a s itself (see e.g. Q). The pre- 
cision of a s extractions from DIS jet cross sections is 
currently limited by the availability of theoretical predic- 
tions only at NLO Q. 

Predicting the dependence of such cross sections on 
jet algorithms, sizes, and vetoes to high accuracy cur- 
rently presents a formidable challenge. The dependence 
on more "global" observables characterizing the jet-like 
structure of final states can often be predicted to much 
higher accuracy. Indeed, some of the most precise ex- 
tractions of a s today come from hadronic event shapes 
in e + e~ collisions, for which theoretical predictions in 
QCD exist to N 3 LL accuracy in resummed perturbation 
theory matched to N 3 LO fixed-order results j8l4l4|. along 
with a wealth of data from LEP. Using event shapes to 
describe jet-like final states in QCD in a global manner 
holds the potential to improve the description of jet pro- 
duction in DIS to the same high level of precision. 

Thrust distributions in DIS were considered in fl5| and 
calculated to NLL accuracy in resummed perturbation 
theory, matched to 0(a s ) fixed-order results analytically 
and compared to O(oq) fixed-order results calculated nu- 
merically Since then the improvement of theoret- 
ical predictions for DIS event shapes beyond these orders 
of accuracy has not received much attention. The intro- 
duction of soft collinear effective theory (SCET) [H-[22| 
has brought about a revolution in methods to achieve 
higher-order resummation in a variety of applications in 
QCD, leading, for example, to the N 3 LL resummation 
of thrust 12 1 and heavy jet mass [l3| in e + e~ collisions 
mentioned above . SCET has been used to predict a wide 



variety of event shapes in e + e~ collisions [23H25| and 
pp collisions [26l - [29j . going beyond the resummed accu- 
racy previously available. A wealth of data now exists on 
event shapes in DIS from measurements at HERA by the 
ZEUS and HI collaborations (30M35| . To take advantage 
of these data, for instance to achieve high-precision ex- 
tractions of a s , requires commensurate accuracy in theo- 
retical predictions. Thanks to advances already made in 
tools and calculations for e + e~ and pp event shapes, the 
time is ripe to extend the accuracy of DIS event shape 
predictions beyond NLL. (DIS in the endpoint region, 
x 1, has been studied with SCET in HISS].) 

Traditional ways to define jet cross sections involve the 
use of a jet algorithm (such as bp-type recombination al- 
gorithms or infrared-safe cone algorithms (4l| - |46j ). and 
often a jet veto as well. Predicting the dependence on jet 



algorithms, sizes, and vetoes to high accuracy is currently 
a formidable theoretical problem in QCD. In particular, 
non-global logarithms (NGLs) [13, EH can arise and com- 
plicate resummation beginning at NLL order for observ- 
ables that probe soft radiation with different measures 
in sharply divided regions of phase space, as occurs with 
some jet vetoes, for instance [49l - l55j . Similar clustering 
logs due to the way algorithms cluster soft gluons can also 
spoil resummation beginning at NLL order 
[59j . NGLs and clustering logs limit the precision one can 
achieve in theoretical predictions for jet cross sections in 
QCD. A great deal of progress has been made to resum 
NGLs numerically in the large- Nc limit [13, E3 , to un- 
derstand the origin and structure of NGLs in the frame- 
work of effective field theory [EH, Io0l - l62l ] , and to find ways 
to minimize their numerical impact (e.g. [29l. l63j ) . but a 
generic approach to obtain NNLL and higher order pre- 
dictions does not yet exist. These complications due to 
non-global methods of measuring jets provide a strong 
motivation to use global measurements of hadronic final 
states that still probe their jet-like structure and are re- 
summablc to arbitrarily high accuracy in QCD pertur- 
bation theory. The first steps needed for higher order 
resummation in DIS are the derivations of appropriate 
factorization theorems. 

Precisely such a global measure of jet-like structure of 
hadronic final states is the iV-jettiness introduced in 27 1. 
./V-jettiness Tjy is global event shape that is a generaliza- 
tion of thrust [64| and can be used in any type of collision 
to constrain the final state to contain N + Nb jets, where 
Nb is the number of initial-state hadronic "beam" direc- 
tions. In e + e~ collisions, events with small tn contain 
N jets in the final state; in pp collisions, they contain 
N + 2 jets, with two along the beam directions from ini- 
tial state radiation (ISR). In DIS, small r^r constrains 
events to have N + 1 jets, with one jet along the beam 
direction from ISR from the proton. 

In this paper we will predict a special case of N- 
jcttiness cross sections in DIS, the 1-jcttincss. We define 
a whole class of DIS 1-jcttiness observables by 
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Tl = Q2 mm {9B ' Pi, Q-J ■ Pi} > 



(1) 



iex 



where qb is a four-vector along the incident proton beam 
direction and qj is another four-vector picking out the 
direction of the additional final-state jet we wish to mea- 
sure. Particles i in the final state X are grouped into 
regions, according to which vector qs,j they are closer 
to as measured by the dot products in Eq. ([TJ. Different 
choices of qs.j give different definitions of the 1-jettincss. 
In this paper we consider three such choices: 



r{ 1 : 
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= xP, 


Q a j 


= jet axis 


(2a) 






= xP, 


<A 


= q + xP 


(2b) 






= p, 


</j 


= k, 


(2c) 



where P and k are the initial proton and electron mo- 
menta, and Q and x are the usual DIS momentum trans- 
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fer and the Bjorken scaling variable. The three versions 
of T\ in Eq. ^ are named for one of their distinctive 
properties: rf aligns the vector with the physical jet 
axis as identified by a jet algorithm or by minimization of 
the sum in E q. (CD ) over possible directions of see for 
example Ref. [65j . This jet axis is almost but not quite 
equal to q + xP, which is used as the vector q b j in t\. The 
measurement of rf groups final state particles in Eq. (CD) 
into exact back-to-back hemispheres in the Breit frame. 
Finally, rf groups particles into exact back-to-back hemi- 
spheres in the center- of -momentum frame. 

Note that the three T\S in Eq. are physically dis- 
tinct observables. Each one of them can be defined in any 
reference frame, but the definitions may be simpler in one 
frame versus another. The DIS 1-jettiness rf coincides 
with the version of 1-jettiness recently considered in (6(| 
at NLL order, and is closest in spirit to the original N- 
jettiness event shape in (27|. No factorization theorems 
so far exist for cither rf or rf. 

There arc in fact a number of DIS event shapes that 
have been measured by experiments at HERA. Two ver- 
sions of thrust [U were measured by the HI Collabo- 
ration H3-H1], and by the ZEUS collaboration [Hfsi]. 
The DIS thrust variables t 7 \e are all based on hemi- 
spheres in the Breit frame where the thrust axis h is 
cither frozen to z (along the virtual 7 or weak boson), or 
determined from a minimization. They have been com- 
puted to NLL+0(o; s ) [HI, |48|. The t^e measure par- 
ticles from only one hemisphere, and the choice of nor- 
malization E determines whether they are global or non- 
global [48| (where the non-global variables were used for 
the experimental measurements). Our 1-jettiness event 
shapes defined in Eqs. (CD) and ^ arc global variables, 
avoiding NGLs by including information from all parti- 
cles in the final state. We will demonstrate that our DIS 
1-jettiness variable rf actually exactly coincides with the 
DIS thrust tq = tzq, computed in [15| at NLL+0(a s ). 

It would be interesting to re-analyze HERA data to 
measure global 1-jettiness or thrust variables. For such 
measurements, one may be concerned about the contribu- 
tion of the proton remnants to Eq. (CD)- However, the rem- 
nants remain close to the q B axis, so their contributions 
to the sum giving n are exponentially suppressed f67j . 
It is only the larger-angle soft radiation and ISR in the 
beam region and the collision products in the qj region 
that need to be measured. In fact, we will show below 
that one can measure r°' b ' c only from the products in 
the qj region, obtaining the g^-region contributions by 
momentum conservation (however for t° this is true only 
in the two-jet region rf <C 1). 

We will give predictions for cross sections in the three 
versions of t\ in Eq. © accurate for small t\. We will 
also prove factorization theorems for all three variables 
rf ,b ' c . The structure of these factorization theorems will 
differ because r°' b ' c each probe initial- and final-state ra- 
diation in DIS differently. Besides grouping final-state 
hadrons into different regions, each version has a differ- 
ent sensitivity to the transverse momentum of ISR. For 



rf and rf , the nonzero k± of ISR causes the final-state 
jet momentum to deviate from the qj axis by an amount 
~ due to momentum conservation. This affects the 
measurement of rf or rf at leading order. For rf , q°j is 
always aligned with the physical jet momentum and so is 
insensitive to the k± of ISR at leading order. This leads 
to different structures in the factorization theorems for 

T l • 

We will prove that the cross sections in all three vari- 
ables factorize as special cases of the form: 



da 



dcr Q 



2 ^H K {Q 2 ,n) Jdtjdt B dk s d 2 Pl _ 



dx dQ 2 dri dx dQ 2 

x J q (tj - (q± + p±) 2 , /i) B K/p (t B ,x, p 2 x , n) 

x S hcm i(k s ,n) s(ti - — - — - , (3) 
V sj s B Q R J 

where k runs over quark and antiquark flavors, sj, s B , Q B 
are normalization constants given in Eqs. (|54p and f|58[) 
that depend on the choice of observable t\ in Eq. ©. 
(To is the Born cross section, H K is a hard function aris- 
ing from integrating out hard degrees of freedom from 
QCD in matching onto SCET, J q is a quark jet function 
describing collinear radiation in the final-state jet, and 
B K / p is a quark beam function containing both pertur- 
bativc collinear radiation in a function X K j as well as the 



proton parton distribution function (PDF) /. 



j/p- 



B K / p (t,x,p 2 L ,n) =y^[ —X K j(t, -,p±,fi) f j/p (z,fj,). 
j 

(4) 

This beam function depends on the transverse virtuality 
t B of the quark k as well as the transverse momentum 
of ISR. Shomi m Eq. © describes soft radiation from 
both the proton beam and the final-state jet. Despite 
the fact that the 1-jettiness Eq. (CD) may not divide the 
final state into hemispheres, we will nevertheless show 
that the soft function for any 1-jettiness in DIS is related 
to the hemisphere soft function Shomi- Finally, is the 
transverse momentum of the momentum transfer q in the 
DIS collision with respect to the jet and beam directions. 

We briefly discuss differences in the factorization the- 
orem for r°' b,c . For rf , the jet axis is aligned so that 
the argument of the jet function tj — (q± + pjj 2 — > tj 
with zero transverse momentum, and p^ then gets av- 
eraged over in Eq. (C|]), removing the dependence on this 
variable in the beam function and yielding the ordinary 
beam function of Ref. 



67] 



For r°, 



rf, the convolution 
over pj^ remains and thus they are sensitive to trans- 
verse momentum of ISR. Thus for Tf ,rf results depend 
on generalized p^-dependent beam function introduced 
in Ref. [6|| . The final difference is that is identically 
zero for , while it is nonzero for rf, causing these ob- 
servables to differ and inducing additional complications 
in the convolution over p^ for rf . In particular the cross 
section for rf docs not start at rf = 0, but rather at 
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T i = qI/Q 2 due to the nonzero q_i_ injected into the 
collision and the choice here for the jet axis. 

The ingredients in the factorization theorem Eq. ([3]) 
depend on an arbitrary scale [i that arises due to inte- 
grating out degrees of freedom from QCD, matching onto 
a theory of collinear and soft modes, and then integrat- 
ing out collinear degrees of freedom and matching onto 
just soft modes. The resulting hard, jet, beam, and soft 
functions each depend on logs of \i over physical vari- 
ables. Renormalization group (RG) evolution allows us 
to evolve each function from a scale hh,j,b,s where these 
logs are minimized to the common scale fj.. This evolu- 
tion resums logs of r°' b,c to all orders in a s , to a given 
order of logarithmic accuracy determined by the order 
to which we know the anomalous dimensions for the RG 
evolution. We will use this technology to resum logs of 
1-jettiness in DIS to NNLL accuracy for T^' b ' c . 

The factorized cross section in Eq. © accurately pre- 
dicts the Ti distribution in the peak region and for the 
tail to the right of the peak, where t\ -C 1 and logs of 
Ti are large. To be accurate for larger n, the prediction 
of Eq. ((3]) must be matched onto predictions of fixed- 
order QCD perturbation theory to determine the "non- 
singular" terms. In this paper we do not perform the 
matching onto the 0(a s ) and C(a 2 ) tail of the t\ distri- 
butions, leaving that to future work. However, by com- 
paring the unmatched predictions of Eq. ([3]) integrated 
over n to the QCD total cross section at x, Q 2 we can es- 
timate the small size of these missing corrections at large 
Ti. We emphasize that Eq. Q accurately captures the 
distribution for smaller ti near the peak region. 

The factorization theorem Eq. (|3|) also allows us to ac- 
count for nonperturbative effects — not only in the parton 
distributions f(x,/j,), but also through a shape function 
that appears in the soft function S. In e + e~ collisions, 
the leading nonperturbative corrections from this shape 
function have been shown to be universal for different 
event sha pes and collision energies [69l - l72j (for earlier 
work see |73l - t75| ). The same conclusions hold for the 
soft shape function in Eq. Q , endowing it with real pre- 
dictive power. We will analyze the dominant effects of 
the nonperturbative soft shape function on the DIS 1- 
jcttincss. For the peak region we include a simple non- 
perturbative model function to show the impact these 
corrections have and how they modify the perturbatively 
calculated distribution. For the tail region the leading 
shape function power correction is a simple dimension- 1 
parameter fi°' b ' c that induces a shift to r"' ' c , and is de- 
fined by a matrix element of a soft Wilson line operator. 
For our observables we will prove that there is universal- 
ity for this correction, namely that ft® = fl\ = OJ. This 
follows from a general analysis we carry out for how the 
direction of axes affect nonperturbative matrix elements 
for two-jet soft Wilson line operators. 

The paper is organized as follows: In Sec. [TT] we re- 
view the kinematics of DIS in several commonly used 
reference frames, laying out the notation for our subse- 
quent analyses. In Sec. IIIII we define the three versions 



of 1-jettiness in DIS that we will use in this paper and 
consider their physics in some detail. In Sec. II VI we follow 
the usual formalism for calculating the DIS cross section 
in QCD, and introduce an additional measurement of the 
1-jettiness into the hadronic tensor that appears therein. 
Scc.|V]is the technical heart of the paper. Here we present 
the elements of the SCET formalism that we need and 
give a detailed proof of the factorization theorems for the 
generic DIS 1-jettiness in Eq. (fTJ) and the three special- 
izations we give in Eq. © . In particular we derive in 
each factorization theorem how the observable depends 
on the transverse momentum of ISR through the beam 
function, and also show that by rescaling arguments we 
can always use the hemisphere soft function for each ver- 
sion of 1-jcttiness. 

In Sec. I VII we use the factorization theorems from 
Sec. [V] to give predictions for the singular terms in the 
T\ distributions at fixed order 0(a s ), and also enumerate 
the results for the hard, jet, beam, and soft functions that 
we will need to perform the RG evolution in the next sec- 
tion. In Sec. IVIIl we perform the RG evolution and give 
our resummed predictions to NNLL accuracy. We com- 
pare our predictions for t\ to those of [l5[ at NLL. We 
also explain the "profiles" for the individual hard, jet, 
beam, and soft scales which we use to perform the RG 
evolution [TJ, [28|, [76[ . These profiles allow for a smooth 
transition from the tail region into the peak region where 
the soft scale becomes nonperturbative, and into the far 
tail region where the resummation of logarithms must be 
turned off. Then we explain how we incorporate nonper- 
turbative hadronization corrections into our predictions 
through a soft shape function and discuss the Oi param- 
eters. We show that the shifts f2°' b ' c to the tail region 
of all three versions of the 1-jettiness distributions obey 
universality. 

In Sec. I Villi we present numerical results for our pre- 
dictions to NNLL for the T-f' ' c cross sections, including 
also their x and Q 2 dependence. We consider both in- 
tegrated (cumulant) and differential cross sections. The 
particular results we present are for x, Q 2 values studied 
at HERA [H, [35|. However, the analytic results we give 
in Sec. IVIIl can just as easily be used for other experi- 
ments at different kinematics, such as at Jefferson Lab 
(JLab) (77[, or for nuclear states other than the proton, 
such as those at the future Electron- Ion Collider (ETC) 
[Zll and Large Hadron Electron Collider (LHeC) [zl]. 

In Sec. IIXI we conclude. In several Appendices we col- 
lect various technical details that are used in the main 
body of the paper. In particular, in App. [D] we collect 
the anomalous dimensions we need to get to NNLL ac- 
curacy in the Ti cross sections, and in App. [F] we give 
the resummed cross sections in an alternative formalism 
[H, El to that used in Sec. EH] [z|. In Sec. IVHl 
we use a formalism that expresses the result of the RG 
evolution of Eq. (J3)) entirely in momentum space, while 
in App. [F] we use a formalism that expresses the RG 
evolution through Laplace space objects. These two ap- 
proaches give identical analytic results at each order in 



rcsummed perturbation theory, but since both are com- 
monly used in the SCET literature we provide both re- 
sults for people who prefer one or the other. Indeed, all 
of our numeric results have been cross checked between 
two codes which each use one of these two approaches. 

The reader mainly interested in the phenomenology of 
DIS 1-jettiness and our numerical predictions may read 
Sees. [T MlIII and then skip to Sec. I Villi For those inter- 
ested in details of the factorization and resummation, we 
provide these in Sees. HVIIVIII and the Appendices. 



II. KINEMATICS OF DIS 

In this section we define the kinematic variables in DIS 
that we will use throughout the paper. We also consider 
three reference frames — center-of- momentum (CM), tar- 
get rest frame, and Breit frame — and describe the picture 
of the events in each of these frames. 



A. Kinematic variables 



where the total invariant mass s = (P + k) 2 = 2P-k and 
y ranges from < y < 1. The variable y measures the 
energy loss of the electron in the target rest frame. For 
a given s Eq. © relates x, y, and Q 2 to one another, 
allowing one of the three variables to be eliminated. The 
invariant mass of the final state in terms of the above 
variables is 

V 2 x = ^—^Q 2 = {l-x)ys. (10) 

x 

In the classic DIS region one has p 2 x ~ Q 2 for generic 
x. In the endpoint region 1 — x ~ Aqcd/Q, the final 
state is a single narrow jet with momentum of order Q 
in the virtual boson direction (and studied with SCET 
in Refs. [36l - l40l |). The resonance region where 1 — x ~ 
Aq CD /<5 2 cannot be treated with inclusive perturbative 
methods. 

In this work we are interested in the classic region 
where 1 — x 3> Aqcd/Q i-e., x ~ 1 — x < 1. In this 
region one can have more than a single jet. Below, we 
will make an additional measurement that picks out two- 
jet-like final states. 



In DIS, an incoming electron with momentum k and 
a proton with momentum P undergo hard scattering by 
exchange of a virtual boson (photon or Z) with a large 
momentum q, and outgoing electron k' . The boson mo- 
mentum q can be determined from the initial- and final- 
state electron momenta, 



q = k — k . 



(5) 



In inclusive DIS, the final states from the hard scat- 
tering are inclusively denoted as X and their total mo- 
mentum is denoted as px- Using Eq. ([5]) momentum 
conservation k + P = k' + Px can be written as 



P = Px- 



(6) 



The momentum scale Q of the hard scattering is defined 
by the virtuality of the exchanged gauge boson. Because 
the boson has a spacelike (negative) virtuality, one de- 
fines the positive definite quantity Q 2 by 



(7) 



where we will be interested in Q > Aqcd- Next one 
defines dimcnsionlcss Lorcntz invariant variables. The 
Bjorken scaling variable x is defined by 



Q' 2 



2P-q 2P-q 

where x ranges between < x < 1. 
invariant quantity y is defined by 



(8) 



Another Lorentz- 



B. Center-of-momentum frame 

A two-jet-like event in the CM frame is illustrated in 
Fig. [T] An incoming electron and proton collide and pro- 
duce in the final state an outgoing electron and hadrons. 
The hadrons, mostly collimated into two jets with addi- 
tional soft particles elsewhere, are grouped into two re- 
gions "Kb and Hj, and pj and pb are the total momenta 
of particles in the each region. The regions Hb,j are not 
necessarily hemispheres in this frame, though we drew 
them as such in Fig. [1] The definitions of the regions are 
described in Sec. IIII Al As shown in Fig. [TJ the electron 
direction is defined to be the +z direction and the proton 
direction to be the —z direction. In the CM frame the 
initial electron and proton momenta are 



2 



(11) 



where the light cone vectors are 



n z = (1,0,0,1), n z = (1,0,0,-1). (12) 

They satisfy n z ■ n z = 2 and n z ■ n z = n z -n z = 0. An 
arbitrary four vector can be written as 



y» = y+_JL + y-_±_ + V g ^ 



(13) 



where V+ = V-n z and V~ = V-n z and V$ = -Vf, < 0. 
In this frame x, y take the values 



PJ 



FIG. 1: Two-jet like event in center-of-momentum frame, in 
which one jet is produced by initial state radiation from the 
proton, and the other by the hard collision with the electron. 
Particles are grouped into two regions Hj : b with total mo- 
menta pj,s in each region. Different choices of "1-jettiness" 
observables will give different boundaries for the two regions. 



and so q is given by 



(15) 



which satisfies Q 2 



xys. Here qt is a four- vector 



transverse to n z ,n z and satisfies = — q^ < 0. 



C. Target rest frame 

The same two-jet like event in Fig. [T] is illustrated as 
it would appear in the target rest frame in Fig. [5] The 
proton is at rest. The regions in Fig. [2] are transformed 
from those in Fig.[T]because of the boost along the proton 
direction. In this frame, the initial electron and proton 
momenta are 

k » = JLH^L pn = M n * +U * , (16) 
M 2 2 v ' 

satisfying 2k- P — s. Here M is the proton mass. We 
reach this frame by a boost of momenta p 11 in the CM 
frame along the z direction, 



M 

n z -p -> —= n z -p, 



n z -p 



M 



n z -p. 



(17) 



Therefore, in this frame, g M in Eq. (|15p is boosted to 
become 



— z - - xM 1 - -% — + dt 

xM 2 \ Q 2 2 lT 



and x, y are given by 



E-E' 



V = 



E 



(18) 



(19) 



2M(E - E') 

Here E and E' are the energies of the incoming and out- 
going electron, respectively, measured in the target rest 
frame. Here y is the fractional electron energy loss. 



FIG. 2: Two-jet like event in target rest frame. The regions 
T-Lj.b and directions of the total momenta pj,b in these re- 
gions are boosted from the CM frame in Fig. [1] Both jets go 
forward, but those in Hj are more highly collimated. 



D. Breit Frame 

In the Breit frame, the virtual boson with momentum 
q^ and proton with momentum collide along the z 
direction. This frame is useful because the proton initial 
state radiation moving along the proton direction can be 
relatively well separated from other scattering products. 
One might worry that an ISR jet, which we want to mea- 
sure in this paper, could be contaminated by the proton 
remnants which arc difficult to separate from ISR. How- 
ever, the 1-jettiness observable in Sec. IIII that we use 
to measure the jets in the final state is actually insensitive 
to this contamination since contributions from the region 
of the beam remnant give exponentially suppressed con- 
tributions to the variable. The contributions from the 
beam region are by far dominated by the initial state ra- 
diation at larger angles. The picture of the two-jet like 
event in the Breit frame is similar to Fig. [T] with incoming 
electron replaced by virtual boson and with the outgoing 
electron removed. 

The Breit frame is defined as that in which the mo- 
mentum transfer q is purely spacelike: 



(20) 



where we align n z to be along the proton direction: 

P^9.<, (21) 
x 2 

The incoming electron has momentum 

k , = Q^ + Q l^l^_ + k , (22) 

2/2 y 2 

where k§, = Q 2 (l — y)/y 2 - The outgoing electron then 
has momentum 

k>» = Q^t'4 + -^- + k T- (23) 
V 2 y 2 

Unlike the CM and target rest frames, where for a fixed 
s the incident momenta are fixed, in the Breit frame the 
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incident momenta are functions of x,y. Thus each point 
in the differential cross section in x, y corresponds to a 
different Brcit frame. 



III. HADRONIC OBSERVABLES 

A. A-jettiness 

To restrict final states to be two-jet-like, we must make 
a measurement on the hadronic state and require ener- 
getic radiation to be collimated along two light-like di- 
rections. An observable naturally suited to this role is 
the A-jettiness [27j |. In our case, with one proton beam, 
1-jettiness t\ can be used to restrict final states to those 
that have two jets: one along the original proton direc- 
tion (beam) from ISR and another produced from the 
hard scattering. Recalling the definition of T\ in Eq. (TTJ : 



of freedom that describe fluctuations collimated in the 
beam and jet regions. Refs. [13, [65| discussed the pos- 
sibility of also minimizing over possible vectors Qr. j to 
give the smallest possible t\ in Eq. (|24p , and Ref. [65( de- 
veloped a fast algorithm to carry out this minimization. 
Here we will take qs,j to be fixed vectors. We will dis- 
cuss several possible choices for qb,j below, each giving 
a different definition of t\ . 

Measuring T\ to be small means the final state has at 
most two collimated jets, one in the qs direction and one 
in the qj direction (irrespective of the exact definition 
of qs and qj). For power counting purposes we will use 
Ti ~ A 2 which defines a small parameter A <C 1 in which 
we will perform the expansion to obtain the leading-order 
factorization theorem for DIS 1-jettiness cross sections. 



1. Ti : 1-jettiness aligned with the jet axis 



Tl = 7p E min {<lB-Pi,qj-Pi} , (24) 

where qB,qj are massless four- vectors chosen to lie along 
the beam and jet directions. 

The minimum operator in Eq. (|24|) groups particles 
in X with the four-vector to which they are closest (in 
the sense of the dot product). We will call the region in 
which particles are grouped with the beam Hb and the 
region in which particles are grouped with the jet Hj. 
We denote the total momentum in the beam region as 
Pb and total momentum in the jet region as pj: 



Pb 



= E 



Pi ■ 



Pj = 



E 



Pi 



(25) 



These regions are illustrated for two examples in the CM 
and target rest frames in Fig. Q] and Fig. [2l 

The 1-jettiness t\ can be expressed as the sum 



Tl 



TB + Tj , 



where tb and tj are defined by 



tb 



2gB-PB 
Q 2 



Tj 



2q.J-PJ 
Q 2 



(26) 



(27) 



The variables tbj are projections oips,j onto the refer- 
ences vector qB,j- They can be thought as two indepen- 
dent obscrvables, and T\ is one possible combination of 
them. Another combination gives a generalized rapidity 
gap and is discussed in App. 

The reference vectors qs and qj can be expressed as 



1b 



uj b 



Lt 

1j 



(28) 



for light-like vectors ns,j given by ns.j = (l,ns,j), 
where iib,j are unit 3- vectors satisfying n| j = 1. Below 
we will use the vectors Hb j to define the directions of 
the collinear fields in SCET which we use for the degrees 



The first version of 1-jettiness that we consider is rf, 
which is defined by choosing the beam reference vector q B 
in Eq. (|24l) to be proportional to the proton momentum, 
and the jet reference vector q'j to be the jet momentum 
as given by a jet algorithm such as anti-kx [Hf: 



■Pi,Qj-Pi}' 



These reference vectors are given by the values 



q%» = xP^ 



1j 



(29) 



(30) 



where qj- is 0(Q\). This is because xP is the longitudi- 
nal momentum of the parton that hard scatters from the 
virtual photon of momentum q, which would produce a 
jet of momentum q + xP, but the colliding parton may 
also have a transverse momentum of order QX. It can- 
not be larger, otherwise it would cause T\ to be larger 
than C(A 2 ). Various jet algorithms give the same value 
of q a j up to negligible power corrections of 0(<3A 2 ), and 
the cross section does not actually depend on which of 
these algorithms is used. Here it would also be equiva- 
lent to leading power to define r° by minimizing the sum 
in Eq. (|29|) with respect to nj in q a j. The total momen- 
tum of particles in the jet region T-Lj is pj = q a j + k for 
a soft momentum k of 0(QX 2 ). Thus, to the order we 
are working, the sum over particles in the jet region T-Lj 
in Eq. (|29|) gives the total invariant mass of those parti- 
cles, 2q°, -Pj = p 2 j = m 2 (for more discussion of this see 
Refs. [29|,l8l|)- 

We will show below that in deriving the correct fac- 
torization theorem for the rf cross section, we must use 
the fact that qj- is chosen to make the relative transverse 
momentum between q°j and the actual jet momentum pj 
be zero (technically the dominant 0{Q\) part must be 
zero and a small 0(QX 2 ) part is still allowed). That is, 
q°j is aligned with the jet, hence the name rf. This is 
also important for experimentally measuring rf. Never- 
theless, once this factorization theorem is known, qj is 
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FIG. 3: (a) 1-jettiness rf measures the small light-cone component of the momentum in the jet region Hj along the "true" jet 
axis q 1 }, which is proportional to the jet invariant mass and is thus insensitive at leading order in A to the transverse momentum 
Pb of ISR. Thus pg gets averaged over in calculating the rf cross section, (b) 1-jettiness Ti measures the small light-cone 
component of pj along the fixed axis q b j — q + xP. This projection is sensitive to and balances the transverse momentum pg 
of ISR. The transverse momenta of pg and pj get convolved together in calculating the cross section. Both rf and r\ divide 
the final state into hemispheres in the Breit frame, (c) 1-jettiness rf divides event into back-to-back hemispheres in the CM 
frame and projects beam and jet momenta onto n z ,n 2 axes. These projections are sensitive to the transverse momentum pg 
of ISR. The momentum transfer q has a nonzero transverse component in these coordinates, and the jet and beam momenta 
are convolved in pg in calculating the cross section. 



not directly required for calculating the objects such as 
hard and soft functions that appear in the factorization 
theorem. For the other versions of 1-jettiness we consider 
below, the reference vector qj is not aligned exactly with 
the jet, and the transverse momentum between qj and 
the jet momentum pj will be nonzero, as illustrated in 
Fig. [3] This will change the structure of the correspond- 
ing factorization theorems, introducing convolutions over 
the transverse momenta of radiation from the beam and 
from the final-state jet. 



2. t{ : hemisphere 1-jettiness in the Breit frame 



A second way to define 1-jettiness in DIS is 



^2 m[n {lB -Piilj 'Pi} : 



where 



hi* 

Qb 



if 



xP^ . 



(31) 



(32) 



In this case, q b j is given exactly by the quantity q + xP 
which can be constructed from the electron and proton 
momenta k, fc', P, and needs no information about the 
jet momentum given by any jet-finding algorithm. Thus 
in general q b j differs by a transverse momentum qj- ^ QX 
from the vector qj used in the rf definition of 1-jettiness 
we introduced above in Eq. (|2"9"]l . Note that since q = 
9 j ~ Qb' Q itself has zero tranverse momentum q± with 
respect to the directions n b j, n b B of q b ,q B . 

This choice of vectors is natural in the Breit frame 
(hence the name rf), in which it divides the final state 
into back-to-back hemispheres. In the Breit frame, 



b Breit 1 . r _ 

7-j = — 2_^mm{n z -pi 



n 



■Pi} ■ 



(33) 



This definition directly corresponds to the thrust tq in 
DIS defined in ^ . 

We will often work in the CM frame in intermediate 
stages of calculation below. Expressing q b B j in the CM 
frame, we find 



b A* r' v z 
q B =Xy/s — , 



(34) 



where = (1 — y)Q 2 and q b is a massless vector. q b ir 
Eq. can also be written in the form 



q b f = P T e Y ^-+P T e- Y ^ + P T n^, (35) 



where the jet transverse momentum and rapidity are 



i 



Y = - In 



2 x(l-y)' 



(36) 



and fiT is a unit vector in the direction of qx- These 
relations can be inverted to give 



Pre 



-Y 



V 



P T e> 



(37) 



Equating the Oth components of Eqs. (f2"5)l and (|35[) . we 
find that 

Lu b j = 2P T coshr = [y + x(l -y)]Vs. (38) 

Calculating t\ in the CM frame groups particles into 
non-hcmisphcrc-likc regions. Particles with momenta p 
are grouped into the beam or jet regions according to 
which dot product is smaller: 

xy/s n| -p uj b jn b j-p 

li-B ■ „ < -Z , 



n ^ _ Xy^n b B -p > uj b jn b j-p 



(39) 
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Using Eq. (|38|). we can write these conditions as 



Hi 



n b R -p y 
rij-p x 



x 



n b j -p 



(40) 

In order to understand the regions defined by Eq. (1401) . 
let us consider simple case y ~ 1 and x < y. For this 
case q b j in Eq. (|35|) is n z -collinear because in Eq. (f36|) 
Pt and Y are small and large, respectively. We can 
replace rij and n b B in Eq. (|4"0"|) by n z and n 2 and set 
n z -p/n z -p = l/(tan 2 6*/2) where 9 is the polar angle of 
masslcss particle p. Then, the jet region is a symmetric 
cone around the n z direction of opening angle given by 



tan 



R 



x 

y 



(41) 



and the beam region is everything outside. For generic x 
and y, the jet region is not symmetric around the n b j. 

As mentioned above in the description of rf , the vec- 
tor q b = q + xP is the 4-momentum of a jet produced by 
scattering at momentum q on an incoming parton with 
momentum exactly equal to xP. In general the colliding 
incoming parton will have a nonzero transverse momen- 
tum due to ISR, causing the produced jet momentum 
to deviate by 0(QX) from q b . The scale 0(QX) is per- 
turbative and this transverse momentum is much larger 
than the intrinsic transverse momentum of partons in 
the proton. The observable rf differs from rf in that 
rf measures the true invariant mass mj of the jet while 
rf simply projects the jet momentum onto the fixed axis 
q b = q + xP which does not vary with the exact direc- 
tion of the jet. The jet axis varies from q b due to ISR 
from the beam before the hard collision. This subtle dif- 
ference leads to a different structure in the factorization 
theorems for rf and rf . 

For the 1-jettiness for DIS studied in [f36j], the proce- 
dure for determining the qj was described as determin- 
ing the jet axis from a jet algorithm. This makes their 
q,j correctly correspond with our q a j. However, they also 
used the formulas Eqs. f|35[) and ([36]) to describe their qj, 
which yields q,j = q + xP, and this would correspond to 
our r b . This choice neglects the 0(QX) transverse mo- 
mentum between qj and the jet momentum pj, which 
taken literally would lead to an incorrect factorization 
theorem for the observable rf. However, after the cor- 
rect form of the factorization theorem for t° is known 
(which was written in [6(1), this approximation is valid 
for calculating the objects in that theorem to leading or- 
der in A. Thus, the t\ in [66| is the same as our rf defined 
above in Eq. (|29|). where qj is aligned along pj. 



3. rf; hemisphere 1-jettiness in the CM frame 

A third way to define the 1-jettiness in DIS is with the 
proton and electron momenta 



1-jettiness 


axis qj 


axis qB 


generic n 


nj 

LOT 

2 


UB — 


rf 


xP + q + qj 


xP 


rf 


xP + q 


xP 


rf 


k 


P 



TABLE I: Reference vectors qj and qB defining the axes for 
various versions of 1-jettiness. For rf the q.j axis is defined to 
be the jet momentum qj given by, e.g., the anti-kT algorithm. 
This axis is given by q + xP up to transverse momentum 
corrections of order qj ~ O(QX). The exact value of qj will 
not be needed for our calculation, only the fact that there is no 
relative transverse momentum larger than 0(Q\ 2 ) between 
the momentum pj in the jet region "Hj and the axis q 1 }. This 
is in contrast to rf , for which the cross section will depend on 
the transverse momentum between pj and q b j=q + xP, but 
where q± = 0. Finally for rf we also have q± 7^ 0. 



We use the superscripts c because this choice naturally 
divides the final state into hemispheres in the CM frame, 
mimickin g th e thrust defined in the CM frame for e + e~ 



64|. 



collisions 

In the CM frame the momenta k and P are along the 
z and —z directions as in Eq. ([lip . In this frame the ref- 
erence vectors qj t B are given by the light-cone directions 
n'i 



l b j and normalizations ui B f. 



and 



(43) 



(44) 



In this frame, ri is then given by 

rf = ^ min{n 2 -p u n z -p,} . (45) 



<f B " = P" 



q7 = k" . 



(42) 



The minimum here assigns particles to either the hemi- 
sphere containing the proton or electron. States with 
small t\ thus have two nearly back-to-back jets in this 
frame. 

The essential differences among rf , rf , rf are illus- 
trated in Fig. |3] drawn in the CM frame and summa- 
rized in Table |U t\ and rf project the jet momentum 
onto a fixed axis, and are sensitive at leading order to 
the transverse momentum of initial state radiation from 
the incoming proton, while rf always projects the jet 
momentum onto the axis with respect to which it has 
no transverse momentum, and so measures the invari- 
ant mass of the jet which is insensitive at leading order 
to the transverse momentum of ISR. Table U summarizes 
the choices of reference vectors qj : B for the three versions 
of 1-jcttincss defined in this section. 
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B. Versions of DIS Thrust 

Several thrust DIS event shapes have been considered 
in the literature [82| . and some of them have been mea- 
sured by experiments. One version, called tq in [l5| but 
not yet measured, is defined in the Breit frame by 



Broit 2 ^— \ 



(46) 



where He is the "current hemisphere" in the direction 
set by the virtual boson q. We will show below in sec- 
tion UTTD] that tq is equivalent to our t\ f. 

Another version of thrust, used in [30, HH and called 
r t E in jUj, is defined using a thrust axis whose definition 
involves a maximization procedure over particles in the 
current hemisphere He = Hj in the Breit frame: 



tie 



(47) 



The maximization aligns the vector n with the direction 
of the jet in the current hemisphere, just like the q°j vec- 
tor in our definition of t°. However, because the sums in 
both the numerator and denominator arc limited to He, 
the observable is actually non-global [48[ , cutting out ra- 
diation from the remnant hemisphere. 1 Thus it differs 
from our r° which sums over both hemispheres. It can- 
not be simply related to a global version of 1-jettiness as 
above. A global thrust event shape, r t Q, can be obtained 
by replacing the denominator in Eq. (|47|) by Q/2, but 
this version of the thrust event shape is also not related 
to our t°. 

Yet another variation is t z e [13, S| which is like 
Eq. (|47p with the same normalization, but with respect 
to the z-axis in the Breit frame. It is also not global 
[48| . HI and ZEUS have measured t zE and r tE [32l.[35j. 
It would be interesting to reanalyze the data to measure 
the global observables r°' b ' c we predict in this paper at 
NNLL order. 



Jet and Beam Momenta 



Jet and beam contributions to 1-jettiness 



the three different versions of 1-jettiness T^' b ' c . The ex- 
pression n in Eq. (|26p can be written in terms of nj ■ pj 
and riB-pB as 



Tl = 



nj-pj , n B -p B 



Q 



J 



where Qj and Qb are given by 



i - Q2 o ~ Q2 

\J — , HB — • 

UJj U B 



(48) 



(49) 



Table HTl lists explicit expressions for Qb,j m the CM, 
Breit, and target rest frames for the three versions of 
1-jettiness rf,Ti,rf. 



For the three different cases t°' ,c of Eq. (|4"8")) . the con- 
tributions nj - pj will be with respect to different vectors 

a,b,c 



"J 



, and nj ■ pj,iiB ■ Pb will include momenta of parti- 
cles in different regions Hj : b in the three cases. For rf , 
the differences between energies ui b and u° and between 
unit vectors n b j and n°j are of order A since the vectors q b j 
and q a j differ due to the transverse momentum of ISR of 
order QX. So using the same expression uij in Eq. (|38|) for 
ujj is correct up to corrections suppressed by A that can 
be neglected in computing t". Nevertheless, the values of 
nj-pj in the equations for t\ and rf do differ at leading 
power (nj-pj ~ QX 2 ) because the 0(A) difference in the 
axes n b j and n a T is dotted the into transverse momentum 



in pj which is of 0(QX). This difference is reflected in 
the different factorization theorems for rf and t\ . 

The discussion on the jet and beam regions Hj.b in 
Sec. IIII A 21 can be done for a generic t\. For particles 
with momenta p grouped into the beam or jet region, 
the criteria qj ■ p < qs • p and qs • p < q.j • p that define 
the regions Hj.b, respectively, can be written 



^ v UJ 'P ^ u B nj-n B 
P € Hj : < 



peH! 



nj-p 
n B -p 
fiB-p 



< 



2ojj 
uJjnj-riB 

2ujb 



(50a) 
(50b) 



Here nj and % are the normalized conjugate vectors to 
nj and ns, respectively. Their definitions are 



n 5 - ^- 
nj-n B 



-a _ ^ n j 
n B = 



nj-n B 



(51) 



The cross sections for the different versions of 1- 
jettiness in Sec. IIII Al will all be expressed in terms of 
beam, jet, and soft functions that depend on the projec- 
tions of the total momenta in the regions Hb and Hj 
onto the reference vectors qB,j in the definition of the 
1-jettiness Eq. These vectors point in the direction 

of light-cone vectors ns = n z and nj, which varies for 



1 The variable TtE is also not IR safe without a minimal energy 
constraint on the T-L c hemisphere. 



chosen so that nj ■ nj = ub ■ Ub = 2. The parameters 
Rjb characterize the sizes of the regions Hj,b into which 
the 1-jettiness Eq. (|24|) partitions final-state particles. 
The variables on the left-hand sides are analogous to the 
ratio of momenta related to rapidity: n-p/n-p = e~ 2Y for 
back-to-back directions n, n. They can be interpreted as 
a generalized rapidity, e~~ or e~ 2Yn B™B as defined 

by Eq. (|A2[) . These rapidities are defined in terms of 4- 
vectors hj^b and uj^b, which are not in general back-to- 
back. Rj : b in Eq. (|50p characterizes the range of these 
generalized rapidities that are included in each of the 
regions Hjb- 
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_1_-J c L Llllcbb 


frame 


D T 




B 7 




^ J? 


SJ 


SB 




















generic n 








/uj B nj-n B 


ju)jnj-n B 


Q 2 


qB ' q Q 2 


qrq Q 2 








1 


V 2a; s 

1 


\/2qj-q B 


qB ■ qj 


qB ■ qj 




CM 




\£° 














y + x(l — y) 




y + x(l — y) 












Breit 


Q 


Q 


1 


1 


Q 


Q 2 


Q 2 




Target-rest 






iM 












CM 




^fxyQ 


l 


1 










Breit 




xQ 






y/xyQ 


vQ 2 


xyQ 2 


2-2/ 


2-y 






Target-rest 




Q 2 

M 


M 

7^ 


M 









TABLE II: Kinematic variables characterizing 1-jettiness. Normalizations Qj and Qb in the expression Eq. (|48[l and sizes 
Rj,b of the jet and beam regions Hj, B in Eq. (|50p for the different versions of 1-jettiness, in three different reference frames 
described in Sec. ITT] and the Lorentz invariant combinations Q B . = Qj/Rj = Qb/Rb in Eq. (|54[) and s./,s given in Eq. (|58ll . 



Invariants for 1-jettiness 

For later purposes we will express Eq. (|48l) in terms 
of separate nj-collincar, ns-collincar, and soft contribu- 
tions: 



Tl 



n j ■ (p j + k j) n B - {p c B + k B ) 



Qj 



Qb 



(52) 



where p c j is the total momentum of all nj-collinear 
modes, p c B is the total momentum of all ne-collinear 
modes, and kj^ B are the total momenta of soft modes 
in regions rlj, B , respectively. These modes are defined 
by the scaling of their light-cone components of momen- 
tum: 



nj-collinear : 
ns-collinear : 
soft : 



(nj-p,nj-p,p±) ~ Q(A 2 , 1, A) 
(n B -p,n B -p,p±) ~ Q(A 2 , 1, A) 

k ~ QX 2 . (53) 



The normalization constants Qj tB in Eq. (|52[) are not 
Lorentz invariant (which for SCET corresponds to a repa- 
rameterization invariance [83l HH), but by combining 
them with other kinematic quantities we can form in- 
variants in terms of which we can express Eq. (|52l) . One 
set of such combinations uses Rj, B in Eq. (|50l) . The 
sizes Rj >B of the regions are not Lorentz-invariant — 
they depend on the choice of frame. However, the ratios 
Qj/Rj and Qb/Rb are Lorentz/reparameterization in- 
variant and, in fact, are equal: 



Qj_ 

R.i 



to 
Rb 



\/2q. r q B 



(54) 



Expressions for Rb.j and Qr for each case t"' & ' c are given 
in Table [TTJ (Strictly speaking, dot products with are 
not Lorentz-invariant due to dependence on the jet algo- 
rithm, but for calculating Qr and sj^ B we can use the 
approximation q j = q j = q + xP to leading order in A, 
which does give Lorentz-invariant dot products.) 

It is useful to re-express the soft contribution in 



Eq. ([52]) by rescaling the vectors nj yB by 
which gives us 



L J,B 



nj-kj n B -k B n'j ■ k j + n' B ■ k B 
T s = — 1 = — 7^ — ~ • (55) 



Qj 



Qb 



Qr 



This relation will help us simplify the soft function in 
the factorized t\ cross sections later on. This is because 
rewriting the particle grouping in Eq. (|50p in terms of 
n'j B absorbs the factor R,j tB giving n'j ■ p/n'j ■ p < 1 
and n' B ■ p/n' B ■ p > 1. Hence with these variables the 
hemispheres "Hj,b are symmetric, which makes it possible 
to connect our soft function to the usual hemisphere soft 
function. 

We can also re-express the rij B collinear contributions 
to Ti in Eq. in terms of another set of Lorentz- 

invariant combinations involving Qj, B . In the T\ factor- 
ization theorems we derive below, the arguments of the 
collinear jet and beam functions appearing therein will 
naturally depend on "transverse virtualities" n ■ p n ■ p of 
the nj-collinear jet and of the struck parton in the pro- 
ton, respectively. Relating the rij-collincar contribution 
to Ti to the transverse virtuality tj of the jet, 



nj-pj _ nj-pjnj-p'y 



t., 



nj-pjQj nj-qQj 



0(A 4 ), (56) 
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where in the middle step we simply multiplied top and 
bottom by the large component fij ■ pj of the total 
collinear momentum in region Hj, and in the last step 
we used in the denominator hj ■ pj = fij ■ q + 0(QX 2 ). 
The large component of the jet momentum can only come 
from the momentum transferred into the collision by the 
virtual boson of momentum q — the proton with which it 
collides only has a large component in the nj ■ p compo- 
nent. Similarly, the ns-collinear contribution to T\ is 



Qb 



-n B -Px n B -p x 
n B -Px Qb 



t B 



-riB-qQB 



(57) 

where p x is the momentum of the parton that is struck 
by the virtual boson of momentum q. In the middle step 
we used that riB-p% = —nB-Px since the struck parton 
recoils against the ISR and balances the small component 
of momentum in the ns direction. In the last step, we 
defined the positive virtuality ts = —ub ■ Px n-B • Px of 
the spacelikc struck parton and in the denominator used 
that fis-Px = ~ns-q + 0(Q\ 2 ). This is because the 
collision of the virtual boson and struck parton is the nj- 
collinear jet which has no large momentum in the ub • p 
component. Thus momentum conservation requires that 
the large components of ns ■ q and Ub ■ Px cancel. 

The quantities in the denominators of the relations 
Eqs. (f5T>f and (|5"7| are Lorentz invariant: 



Sj = nj-qQj 



SB 



-n B -qHB 



qB-qj 

-Qj-Q 

qs-qj 



Q 



(58a) 
(58b) 



where the minus sign in sb makes it positive since ns-q < 



0. For the cases r^' b ' c , sj and sg take the special values 
given in Table UB 

Using the definitions of Qb an d sj,s in Eqs. (|5"4")l and 
(|58f these factors can be combined to give the transverse 
virtuality of the exchanged boson q: 



sjs B 



Q\ 

where we used 



riB-qnj-q 



^ = Q 2 (l-qi/Q 2 ), (59) 



q = n B -q- 



n B 



nj 



(60) 



and q 2 = —Q 2 . The transverse momentum q± is orthog- 
onal to n B< j. The relation Eq. (|59|) will be useful in 
evaluating the fixed-order n cross section in App. [G] We 
will use that q^/Q 2 ~ A 2 when 1-jcttiness is measured 
to be small, t\ ~ A 2 . A larger qj^ cannot be transferred 
into the final state for this to be true, since particles have 
to be collimated along qj^s or be soft. 



e~' Ay l of their rapidity with respect to q B . Only the 
energetic ISR and soft radiation at larger angles in Hb 
contribute to n . Although these contributions are easier 
to measure, one may still prefer to measure particles only 
in the Wj jet region in the direction of qj. In general, 
such a restriction in the final state is non-global, and 
leads to NGLs. However, by momentum conservation, 
we can show that each of the global r"' b,c observables we 
consider can be rewritten in terms of momenta of parti- 
cles only in the Hj region (for case a this is true only in 
the 2-jet region rf <C 1). 

First, consider t\. In the Breit frame, 



t Broit 1 v . r 
Ti = jz nnn{n z • p it n z ■ pi} 



i E i-P*i)+ ( E i+P*i) 



(61) 



1 

Q 



iex ie u b j 



where X 
that in the Breit frame 



H b j + H b B denotes the entire final state. Note 



(62) 



where p^ x = Y^izxPi- Thus, E x + p z x = Q, and we 
obtain 



b Breit 1 2 ^ _ 



(63) 



where in the last equality we recall that Eq. is pre- 
cisely the definition in Eq. (|46| of the DIS thrust variable 
called tq in [l5|, where the hemisphere T-L b j in the Breit 
frame was called the "current hemisphere" He- We will 
comment further on the relation between the results of 
[HI for tq and our results for rf in Sec. IVII Bl below. 
Eq. (|63[) shows that t\ can always be computed just in 
terms of the measurements of momenta of particles in the 
current hemisphere He = Wj- 

The same arguments as for in the Breit frame apply 
to rf in the CM frame. In the CM frame, 



CM 1 



xyy/s 



^ min{n z -pi,n z -pi} 



(64) 



iex 



xyyfs 



;.ex 



In this frame, we have that 



D. Momentum Conservation and the Beam Region 

We noted earlier that the contribution of proton rem- 
nants to n is exponentially suppressed, by a factor 



Px 



P 



(65) 
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Inclusive DIS cross section 



CM 



1 - 



Pz; 



(66) 



Thus, rf also can be measured just from momenta of 
particles in the Hj hemisphere in the CM frame. 

Finally, the above argument can be extended to apply 
also to the 1-jettiness rf , but only for the region where 
rf <C 1. rf can be written 



2 



E ft'* 



E 



(67) 



Now, q% = q b B , while q°j = q b j + 0(QA). Thus the re- 
gions Hj B differ from those for jf, "Hj B , by a change 
in the region boundary of 0(A). This does not affect the 
assignment of collincar particles to the two regions, since 
none of them change regions under this small change in 
boundary. An 0(A) fraction of the soft particles switch 
from one region to the other, but this then produces a 
correction suppressed by A to the soft contribution t$ in 
Eq. (|55|) . Thus, Eq. (f6"7} can be expressed 



2 

Q 1 



Ew-^+E^+E^ 



tew 



r i + Q2 E W-9j)-^ + 0(A 3 ), 



-0(A 3 ) 
(68) 



in the regime where Ti ~ A 2 < 1. This is the regime 
we aim to predict accurately in this paper. Thus, in this 
limit rf can also be computed just by measuring particles 
in the "current hemisphere" Hj = He in the Breit frame, 
as long as both axes q a j and q b j are measured. For larger 
rf , both regions H°j B would need to be measured, and 
we emphasize that the contribution of proton remnants 
is still exponentially suppressed. 

In summary, for small t\ none of the three versions of 
1-jettiness r°' fc ' c require direct measurement of particles 
from initial state radiation in the beam region. Further- 
more, for larger n values the variables iy c still do not 
require such measurements (though rf does). All three 
Ti's are global observables since measurement of t\ by 
summing over the particles only in the Hj region is still 
affected by ISR from the proton beam through momen- 
tum conservation. 



IV. CROSS SECTION IN QCD 

In this Section we organize the full QCD cross section 
into the usual leptonic and hadronic tensors, but with an 
additional measurement of 1-jettiness inserted into the 
definition of the hadronic tensor. We express it in a form 
that will be easily matched or compared to the effective 
theory cross section wc consider in the following section. 



We begin with the inclusive DIS cross section in QCD, 
differential in the momentum transfer q, 



x (69) 
x (27r) 4 (5 4 (P + q- p x )5 i {q -k + k'), 

where L is the final lepton state with momentum k', and 
X is the final hadronic state with momentum px ■ 
is the phase space for the lepton states, and the m " 
eludes the phase space integrals for hadronic states. The 
squared amplitude \Ai\ 2 is averaged over initial spins, 
and summed over final spins. Recall that q (and x, y) 
can be determined entirely by measurements of the lep- 
ton momenta. Later in Scc. lIVBl wc will insert additional 
measurements such as 1-jettiness on the state X. 

We wish to express the cross section differential in the 
Lorentz-invariant variables Q 2 ,x using Eqs. (J7J) and ((HJ. 
Although Q 2 ,x are Lorentz-invariant, at intermediate 
stages of integration we can work in a particular frame. 
In either the CM or Breit frame, the proton momentum 
is of the form P = n z - P n z /2. So we decompose q along 
the n z ,fi z directions, q = n z ■ qn z /2 + n z ■ qn z /2 + cjt- 
Then the delta functions defining Q 2 , x take the form 



5 x 



n z -Pn z -q 



5(Q 2 + n z -qn z -q~c&) . (70) 



Inserting these into Eq. (|69p and integrating over q + and 
q~ , we obtain 



da 



dx dQ 



^- J d 2 q T J d$> L 5 A {q-k + k') 
xY;(^) 4 S\P + q-px)(\M\ 2 ). 



(71) 



A" 



where q is now given by the value 



Q 2 



+ (^) 



For a single electron final state L = e(fc') (which is 
all we have at the leading order in a em at which we are 
working), the integral over $l in Eq. (|7ip takes the form 



d z k' 



d A k' 



(2n) 3 2E k , J (2tt) 3 
so, performing the k' integral, we obtain 
da 1 



S(k' 2 ), 



(73) 



dxdQ 2 4(2tt) 3 



d 2 q T 8({q-kY) 



J2(^) 4 S 4 (P + q- P x)(\M\ 2 ) 



(74) 



x 



To use the first delta function, we need to pick a partic- 
ular frame in which to complete the qx integration. In 
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Q 2 x ( '_2 ft Q 2 \ 



the CM frame, 



8{{q-kY) = 6(Q 2 +2q-k) = 2- 5[t& - I 1 - — I Q z ) , 

\ \ XS J J 

(75) 

where we use Eq. (j72|) . fc = >/sri z /2 and P = >/sn z /2. 
We use this delta function to perform the q|, integral in 
Eq. (|71|) . and then use that the spin- averaged squared 
amplitude is independent of ip q , to obtain 



Here Ljp, Wjj, depend on x, Q 2 through the components 
of q given in Eq. (f?2"j) . The average over initial electron 
and proton spins is implicit in Eq. ([83]) . as is the sum 
over quark flavors in Eq. (|83b[) . 



1. Leptonic tensor 
The leptonic tensor in Eq. (|83aj) is given by 



da 



Q 2 



dxdQ 2 8(2tt) 2 x 2 s 2 



(84) 

where 



x 



(76) 



Here the integrand is evaluated in the CM frame with q 
now given by 



q^yVs-f-xyVs-f + Vl^yQn^, (77) 

where fix = (0, 1, 0, 0) in (n , n\, n-i, n^) coordinates. 
The matrix element M. is given by 

X(eP-> e'X) = (e'X\Jf EW (0)D^Jl QCD (0)\eP), 
i=i, z 

(78) 

where the sum over / is over photon and Z exchange, 
Ji,ew is the appropriate electron electroweak current, 
Ji,qcd is the quark electroweak current, and D 1 ^ is the 
7 or Z propagator. There is an implicit sum over quark 
flavors. The matrix element can be factored, 

M(eP -+ e'X) = J2 ^'\J^ewDU e )^ X \ J Iqcd\P) ■ 

1=7, Z 

(79) 

More conveniently, we can express the sum over i as being 
over the vector and axial currents in QCD, 



J vf = <ttnr"qj 



Vf-wr Jlf = 9/7^759/ • (80) 



The sum in Eq. ([78]) can then be expressed as 

m= J2 Y, Li ^ x \ J if\ p ) > 

I=V,A f 



(81) 



defining the leptonic vector which contains the elec- 
tron matrix element, electroweak propagator, and elec- 
troweak charges of the quarks implicit in Eq. (f78|) . The 
sum over / in Eq. (|8ip is over quark flavors. 

Now the cross section in Eq. (fT6"|) can be written 



da 



dx dQ' 



= J2 L^(x,Q 2 )W H '^(x,Q 2 ), (82) 



I,I'=V,A 



where 



Q 2 



L1j,(x,Q 2 ) = — (x,Q 2 )L^,(x,Q 2 ), (83a) 
Wft(x,Q 2 ) = J2 (P\ Jf \X) (X\ Jr \P) 



x 



(83b) 



x (2ir) 4 8 4 (P + q- Px ). 



k V A- V h 
T n P v ' n '/j, n 'v 

9uu = 9fiv ~ 



Q 



— f a k a k 113 

(85) 



where k! = k — q, with q given in Eq. (|77|) and 



L 



vv 
gff 



QfQr 



{Qfv f > +v f Qf>)v e v f v r {v 2 + a 2 e ) 



l + m%/Q 2 



(l + m 2 z /Q 2 ) 2 



T vv - -n 



Q f v f , + v f Q f , 
1 + m 2 z /Q 2 

1/ 



2d 



eVfVf'Ve 



t AA ~~J -J ve 1 ~~el t / 

9ff ' (1 + m|/Q 2 ) 2 ' e// ' (1 + m 2 z /Q 2 ) 2 ' 



(l + m|/Q 2 ) 2 ' 
AA _ 2a f afrv e a e 
4/ 



(86) 



7?//' 



-AV 
J <7/' 



J gf'f 



-VA 



0/ 



1 + m 2 z /Q 2 

afa e 
l + m%/Q 2 



y e " 1 + m 2 z /Q 2 



Qf 



2vf>v e 



l + mjr/Q'* 



where we have made explicit the flavor indices /, /'. Qf 
is the electric charge of the quark qf in units of e; Vf,df 
are the weak vector and axial charges of qf, and v e ,a e 
the weak vector and axial charges of the electron. The 
vector and axial charges are given by: 



T 



a f 



f 



sin 28 'w 
1 



<7 



T f - 2Q/SU1 2 9 W 



sin 29 w 
-1 + 4 sin 2 



Mr 



(87) 



2sm29w 2sin2#w 
where Tf — 1/2 for f = u,c,t and —1/2 for / = d,s,b. 

B. 1-jettiness cross section 

To form the cross section differential in the 1-jettiness 
Ti, we insert a delta function measuring t\ into the 
hadronic tensor Eq. (|83b[) : 



da 



dx dQ 2 dri 



= Ll i ;(x,Q 2 )W II '^(x,Q 2 ,r 1 ) 



I,I' = V,A 

where the Tx-dependent hadronic tensor is 

Wff,(x, Q 2 ,n) = (P\ \ x ) ( x \ Jr \ p ) (89) 
x 

x (2^)H\P + q- PX )5(T 1 -T l {X)). 
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Here the 1-jettiness t±(X) of state X is defined by 
Eq. ([2"4"[) . The definition depends on the choices of refer- 
ence vectors qB,.i- 

The sum over states X in Eq. ([89)) can be removed by 
using an operator fi which gives t\(X) when acting on 
the state X: 



FACTORIZATION IN SCET 



f 1 \X)=r 1 (X)\X) 



(90) 



This operator can be constructed from a momentum- 
flow operator as in [HI]. Explicitly, 



'1 ! 



where 



J.B 



2 



dY' qj, B -P{Y') 



(91) 



(92) 



Yj,. 



Here P(Yj B ) is a momentum flow operator that can be 
defined and explicitly constructed in terms of the energy- 
momentum tensor, which can be obtained for massless 
partons using (851488} and for massive hadrons using [72j . 
It measures the momentum flow in the generalized ra- 
pidity direction Yj B , which we define as we did below 
Eq. ([H]) by 



-2Yj 



nj-p 
nj-p 



D -2Y B = n B-P 

n B -p 



(93) 



The lower limits Yj B on the integral in Eq. (|9"2"|) are given 
according to Eq. ([50]) by Yj. B = \h\{l/ R JtB ). These 
values depend on the frame of reference and choice of 1- 
jettiness t\ . For example, for the choice rf of Eq. (|4"2"]) in 
the CM frame for y near 1, the beam and jet regions are 
hemispheres and Yj^ B =0. For the choice of Eq. (|32|) 
in the CM frame, the jet region is given by the lower 
limit Yj = ^ln(y/x), and the beam region is given by 
the lower limit Y B = i \n(x/y). 

In the massless limit the generalized rapidities 
exp(— 2Yj^b) 1 — cos6>j.b/1 — cos9B,jnj -ns /2 defin- 
ing generalized "pseudorapidities" . They depend only 
on angles 8j_b from the nj.s directions and nj ■ n B , and 
so simply characterize angular directions in space over 
which we integrate in Eq. (|9"2"|). 

Using Eq. (j9"2"j) the hadronic tensor Eq. ([59")) can be 
written 

Wff,(x,Q 2 ,Ti) = fd 4 xe^(P\J^(x)S(r 1 - fi)J£(0)|P) , 



(94) 

recalling that q is given by Eq. ([77)) . fi can also be ex- 
pressed in terms of momentum operators in the regions 
Hj.b, using Eq. pg[) : 



n 



nj -pj n B -p B 

Qj Qb ' 



(95) 



where pj b measures the total 4-momentum in region 
Hj.b. 



Soft-collinear effective theory (SCET) |18l-l22| is a sys- 
tematic expansion of QCD in a small parameter A which 
characterizes the scale of collinear and soft radiation from 
energetic massless partons. Soft and collinear modes are 
defined by the scaling of their momenta in light-cone 
coordinates with respect to light-like vectors n,, n (not 
necessarily back-to-back) satisfying n 2 = n 2 = and 
n ■ h = 2. We express the components of a vector p in 
n, fi light-cone coordinates as p = {n ■ p,n ■ p, p±_), where 



n n 
p = n-p- + n-p-+p±_ , 

with p± being orthogonal to n, n, defined as 



P± = 91 Vv , 91 = 9* 



(96) 



(97) 



In these light-cone coordinates, n-collinear and soft mo- 
menta scale as: 

collinear: p n — Q(l, A 2 , A) (98a) 
soft: Ps ~ Q(A 2 ,A 2 ,A 2 ). (98b) 

The parameter A is determined by the virtuality of the 
modes p 2 n ~ Q 2 A 2 that contribute to the observable in 
question. Collinear momenta will be expressed as the 
sum of a large "label" piece and a small "residual" piece: 
Pn = Pn + k, where p n = n-pn/2 +p± contains the O(Q) 
longitudinal and O(QX) transverse pieces, and k is the 
residual 0(QX 2 ) piece. 



A. Matching onto SCET 

Now we are ready to match the currents in Eq. ([94)) 
onto operators in SCET. The QCD current 



with Ty = 
in SCET, 



Jjt f {x)=q f (x)r$q f {x), (99) 
7 M and = 7^75, matches onto operators 

J d 3 p 1 d 3 p2e l{pl - i ' 2) - x 



n±n 2 

X 



neglecting power corrections of 0(A 2 ). The quark and 
gluon SCET operators are 



(101a) 



O^ip^x) = V^Bi% (x)B^_ P2 (x) (101b) 

where we sum over fundamental color indices j and ad- 
joint color indices c, but fix the spin indices a/3 and Xp. 
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We leave implicit that x = Xq carries flavor q. Below 
we will also leave the flavor index / on the current Jj 
implicit. The collincar fields Xn u pi and B^.p. carry label 
momenta 



Pi 



UJiUi 



■Pi, 



(102) 



where i = 1,2. The momentum of each collincar field can 
be written in n%, fii light-cone coordinates as in Eq. (|96j) . 
with the residual x dependence of the SCET fields being 
conjugate to momenta k of order Q\ 2 . In Eq. (|100|) , 
the integrals over pi^ are continuous versions of discrete 
sums over the label momenta, and the measures are given 
by d 3 pi = dujid 2 p\. 

The quark jet fields Xni,p{ x ) are products of collinear 
quark fields with collincar Wilson lines, 



Xn,p 



[8(uj-fi-V)5 2 (p x -Vx)WXtn] 



(103) 



where is a label momentum operator [20j which acts 
on collinear fields and conjugate fields as: 



P^4>n,p , 



P^lp 



and W n is the Wilson line 
W„(x) = ^2 ex P 

perms 



n-V 



i-A n (x) 



(104) 



(105) 



where A^ (x) = J2pA^p(x) is a n-collinear gluon field. 

The gluon jet fields B^ are collinear gauge-invariant 
products of gluon fields and Wilson lines, 

Blp = - g [S(uj + n-r)5 2 (p ± +V ± )W^r ± +gA^W n } . 

(106) 

The matching coefficients C qq , C gg in Eq. (|100p are cal- 
culated order-by-order in a s by requiring that matrix el- 
ements of both sides of Eq. (|100[) between collinear states 
in QCD and in SCET be equal. 

Collincar fields are decoupled from soft fields by the 
field redefinitions [2l[ 

v _ Y v(°) Aarpa _ -yah A<fl)brpa _ y /\(0)brpby\ 

(107) 

where Y n is a Wilson line of soft gluons in the fundamen- 
tal representation. For n = hb we have 



Y n B (x) = Pexp 



>fj 



,o 

/ ds tib ■ A s (ubs + x) 

J — oo 



(108) 



and y nB is defined similarly but in the adjoint represen- 
tation. Soft gluons carry momenta scaling as A 2 in all 
components. Additional factors accompanying outgoing 
states turn the path in Eq. (|108p into a; to oo 89| for 
outgoing collinear particles, see also [90j. So for n = nj 
we have 



y n t» = Pexp 



1 9 



ds nj ■ A s (njs + x) 



(109) 



After the field redefinition Eq. (|107l) . the operators in 
Eq. (fTUTj) become 



0^(p u p 2 ;x) = X { S^T[YlY n2 y\x)x%f 2 (x) , 
0^(p 1 ,p 2 ;x) = ^UT 2 B^ c {x) (110) 



The directions n\ and n 2 will each get set equal to either 
nj or us later on, replacing Y ni with Y„ B in Eq. (|108p for 
ni = Ub or with Y"^ in Eq. (Hoi) for n.j = n,/. Henceforth 
we use only the decoupled collinear fields and drop the 
superscripts. 



The measurement operators in Eq. (|95j) also split up 
into collinear and soft pieces. Since p is linear in the 
energy-momentum tensor, which itself splits linearly into 
decoupled collinear and soft components after the field 
redefinition Eq. (|107[) [85| . p splits up as 



p = p ni + p n2 



(111) 



where p™ 1 '™ 2 ' 8 is built only out of the ni-collinear, n 2 - 
collinear, or soft energy-momentum tensor of SCET, re- 
spectively. 

After matching the product of currents Jj^ Jl in the 
hadronic tensor in Eq. (|94|) onto SCET, there will 
be products of the quark and gluon operators (O q q + 

Ogg)(0' q? + Cgg). The O qq O' gg and O gg O' qq CTOSS tOUlS 

will vanish inside the proton-proton matrix element by 
quark-number conservation (only one of the fields Xn x or 
Xn 2 m will create/annihilate a quark in the collincar 
proton). Thus only the O qq O' qq and O gg O' gg operator 
products can contribute. 

In fact, for DIS, only the quark operator product con- 
tributes, just as in Drell-Yan (DY) |67|. Following the 
arguments in (67j , we know that the matching coefficients 
Cj (pi,p2) must be a linear combination of pi and p%, 
and obey the symmetry 



C Igg„(Pl'P2) 



P~2, -Pi) 



(112) 



due to the structure of the operator Eq. (|101b[) . This 
requires Cj gg to be proportional to (p 2 — which the 
x integration in Eq. (|94p will eventually set equal to q^. 
Vector current conservation in QCD requires qp,C Vgg = 0, 
which requires that Cy be identically zero. The axial 
current matching coefficient can be nonzero, but 

still proportional to p 2 —pi = <? M , which gives zero contri- 
bution when contracted with the lepton tensor Eq. (|84|) . 
Thus for DIS we need only consider the quark operator 
contribution as in DY. 
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B. Factorization of the Hadronic Tensor 

The hadronic tensor Eq. ([94]) can now be written in SCET as 

Wtf (x, Q 2 , n ) = / d A x e^ x Y, I e <(fc-*>-x f dTjdTBdTsS{Tl _ Tj _ TB _ Tg) 

J n 1 n 2 

x (P nB \ Cf^(pi,P2)^T\Y^Y ni ] h 'xZU W{tj - ^p) s (™ - 

x *(r. - ^ - ^)<i^.P' 2 )<i^[47^] /fc '<^(°) l*»> • (H3) 

We have explicitly specified that the proton is an ns-collinear state. The conjugate quark matching coefficient is 
given by Cf"^ (pi,p2) — [7°C/gg(pi 5 P2)7 ] ,3Q - Wc have used that the measurement operator fi can be written in the 
form Eq. (|95j) . and that the momentum operators pj.b split up linearly into purely nj- and ns-collinear and soft 
operators as in Eq. . Wc dropped the subscripts J, B on the collinear momentum operators restricting them to 
the jet or beam regions Hj y B determined by the definition of Ti, since all collinear particles are grouped in region 
B and all nj-collincar particles are grouped in region J. In the soft sector, the restrictions of the operators rj B to 
the Hj.b regions remain. 

Since the nx-collinear, n 2 -collinear and the soft sectors are all decoupled from one another, the proton matrix 
element in Eq. (|113p can be factored, 

W^'(x,Q 2 ,n) = J d+X J cPp 1 cPp 2 e i <** i '-M-* J dTadTEdTsSin-Tj-TB-T^C^^Pi)^^,^) 

x <0| [YlY nj ] kj (x)5(r s - ^ - ^g*) [^^'(0) |0) (114) 
x [(Pn B \xt k B , §2 (x)s(r B ^£l) x ^'(0) \P nB ) (0|x* (rj - ^^)^/(0) |0) 

+ (Pn B \xZ,pM S ^B ^f^)xn/(0) \Pns) (0| X% A (X)5(TJ - !L ^)xg? ' (0) |0)} . 

The last two lines account for the two ways to choose a pair of collinear fields in the proton matrix element. We have 
performed the sums over ni t 2,n'i 2 sums using that the fields within each collinear matrix element must all be in the 
same collinear sector. Wc also require that the fields in the proton matrix clement must be in the same collinear sector 
as the proton, and those in the vacuum matrix element with the direction nj in the definition of t\. The integrals 
over p[ 2 have been absorbed into the definition of the unlabeled fields \ ni 2 . In the soft matrix element we have used 
the fact that T\Y\Y nB \ = FJ Y nB and T\Y 7 \ B Y nj ] = Y^Y nj since the two Wilson lines are space-like separated and 
the time ordering is the same as the path ordering [9ll l92j . For the soft Wilson line matrix element corresponding to 
antiquarks in the beam and jet functions, we have used charge conjugation to relate it to the matrix element shown 
in Eq. (fTT4|) . 

It is measuring t\ to be small that enforces that the direction nj on the collinear fields in the vacuum matrix 
element be equal to the direction of the vector qj in the definition of the 1-jettiness t\ . We are free to choose any 
vector qj to define the observable t\. Requiring that the final-state jet J be close to the direction of qj may, in 
general, impose additional kinematic constraints on x,y,Q 2 to ensure this. We will find below that for r"' b , qj is 
already chosen to be close to the final-state jet and so imposes no additional constraints, while for rf requiring the 
jet be close to q c j = k requires y to be near 1. 

Next we wish to perform the x integral in Eq. (|114l) to enforce label momentum conservation. Before doing so, 
we consider the residual momentum dependence conjugate to the coordinate x in the SCET matrix elements. The 
collinear field Xn,p{x) with a continuous label momentum p depends only on single spatial component n-x because 
the residual momenta (conjugate to the spatial components n-x,x±) arc reabsorbed into p when the discrete label is 
made continuous. Then, the matrix element of n-collinear fields are M n = M n (n-x). For convenience the soft matrix 
element with Y n (x) and Yn(x) will be defined as M s (x). Their Fourier transforms take the form 

/rJm.k , - , ■ 1" d A k ~ 

__ em .kn. x /2 Mn{n . k) ; Ms{x) = J ±* e ik..* Mt ( ka) ; (n5) 

where k, k s is a residual or soft momentum of order QX 2 . When combined with the exponentials containing q or label 
momenta pi,2, we can expand the exponents using q + p + k = (q + p)[l + 0(X 2 )], and drop the terms of order A 2 . 
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Then the remaining integrals over n-k,k s are simply the Fourier transforms of the position space matrix elements 
evaluated at x — 0. So, we can set x = in the SCET matrix elements, and perform the x integral in Eq. (|114[) to 
enforce label momentum conservation. 

In performing the x integration, we have a choice to write x and momenta in nB,ns coordinates or nj,nj coor- 
dinates. In fact, we have freedom to define the vectors tlj^b as long as we choose them such that n 2 B = n 2 = and 
nj-fij = riB-fiB = 2. Since the measurement of n involves measurements of both nj-p and ns-p components of 
particles' momenta, it is convenient to choose hs to be proportional to nj and nj to be proportional to n B , as we 
did in Eq. (|51[) . a choice we will continue to use in what follows. 

For the first pair of collinear matrix elements in Eq. (|114p . the x integral and accompanying phase factor for label 
momentum conservation take the form 

,4 if o+09 -0i f dns-x dfiB-x 2 f ,n B -x n B -x nj-x x _ ± 

<Fxe l(q+P2 = / d 2 x±cxp<(n B -q + uj2)—^—+nB-q— u>i— — + (Q± + P2 ~ Pi 

= 2(2^6 (n B -q + u 2 )5(nB-q- nj ' ^ ) 5 2 (q ± + & - pi ) 

= -^—(2n) 4 5(n B -q + uj 2 )d(nj-q-LJ 1 )d 2 (q 1 _ + p£ - pt) , (116) 
n,j-n B 

where we used Eq. (|51[) to rewrite nj-x in terms of ub-x in the first line and to rewrite ns-q in terms of nj-q in the 
last line. Exchanging U2 and — uj\ in Eq. (|1 16[) gives us the label momentum-conserving delta functions for the second 
pair of collinear matrix elements in Eq. (11 14[) . Using these delta functions to perform the Wx,2 and p^ integrals in 
Eq. (|114[) . we obtain 

W"'(x,Q 2 ,ri) = 2(2tt) 4 (QjQ b ) 2 f d 2 p^^— f dTjdT B dr J s drf 5{n - t, - t b - t s j - rf) 

J nj-riB J 

r/=v/3a fia'p if- nj _ n B \ 

X [ C Iqq» C I'qqA [j 1 ^^ + Q± + P±, ^B^— + P± ) 

x (0| [YlY nj p{Q)5{Q.i^ - nj-fj)8{Q B rf - n B -p B )[YlY n J' k ' {0) |0) 

x l(Pn B \£ k B (0)S(QBT B -nB-p n nWnB-q + n B -T)5 2 {p x -V^x^'^)\ \P nB ) ( 117 ) 



x (0| xZ(0)S(QjTj - nj-p n ')S(nj-q + nj-V)S 2 ( q± + p ± + P ± )j£j (0) |0> 
+ {Pu B I XZ (0)S(Qbt b - n B -p nB ) [S(n B -q + n B -V)S 2 (p ± - V±)xg*' (0)] \P ns ) 

x (0\x n k A°) S (Qj r J ~ nj-rnS(nj-q + nj-V)S 2 ( q± +p ± +V ± )Xnf(0) |0) 

where we use the change of variables p^ = p~± in 3rd and 4th lines and p^ = —p± — q± in 5th and 6th lines. Recall 
that Tlb and hj are now fixed by Eq. (|51[) . The collinear fields without labels implicitly contain a sum over all labels, 
with the delta functions then fixing the labels to a single value (it is important to recall that label operators acting 
on fields Xn,p give minus the label momentum, — p M |20|). The vector nj (may) implicitly depend on the integration 
variable p±, at least for the case of the rf distribution, which we will deal with below. For t\ and rf the vector nj 
is independent of p±. We have also indicated that the arguments of the matching coefficients C, C are both set equal 
to the label momenta of the fields in the collinear proton and vacuum matrix elements. 

The result in Eq. (|117[) is organized in terms of factorized matrix elements that can now be related to known 
functions in SCET. 



C. SCET Matrix Elements 



Beam Functions 



The proton matrix elements in Eq. (|117[) can be expressed in terms of generalized be am f unctions [931.I941] in SCET. 
In covariant gauges (for discussion of similar matrix elements in light-cone gauges see [95j-|97[) they are defined by 



0(f») fdy 



4tt 



lk+ y~/ 2 (p n (p-)\ *n(y-£)f [*(<" - n ■ v)U{k\ vl)xn(o)] \p n (Pl) , (ns) 



BJujk+ 



P 



k\,n) = °^J ^e ik+ y-/ 2 (P n ( P -)\ tr n ■ V)U<& 7>i)x„(0)] \P n (P~)) 
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where the light-cone components of vectors arc given by V + = n ■ V and V~ = n ■ V . Note the dependence 
of the beam functions on the transverse label momentum fcj_ is only on the squared magnitude k\. The matrix 
elements in Eqs. (|118j) and (|119j) arc similar to those that define parton distribution functions, but the separation 
of the collinear fields in the n direction means there is energetic collinear radiation from the proton with virtuality 
~ ujk + 3> Aq CD (assuming we are measuring k + to be large enough), which must be integrated out to match Eq. (| 1 18[) 
onto nonpcrturbativc PDFs (where the separation of x„,x ra fields is zero). The generalized beam functions Eq. (|118p 
are related to the ordinary beam functions originally defined in [67j by integrating over all kj^ : 

= J d 2 k^B q ^[ujk+^-,kl,^ . (119) 

This relationship would be subtle for PDFs, where it is true for the bare matrix elements, but where after renor- 
malization the two objects may no longer be simply related. In the beam function case both sides have the same 
anomalous dimension which is independent of k± and there is no such subtelty. 
The proton matrix elements in Eq. (|117[) can now be expressed as 

(Pn B I xil (0)S(Q b tb - n B -p nB ) [S(n B -q + fi B -V)5 2 (p x - Vx)x^' (0)] \P n „) (120) 

J.P'1 3 gkk' 



(PnMUmQBTB -nB-p^SinB-q + nB-T^ 2 ^ -V±)?g>\oj\ \P nB ) 
- 4T' 5" R ( n B -q _ 2 



where sg is defined in Eq. (|58j) . Now, to simplify the second argument of the beam functions, we note that 

q 2 _ n B -qn B -q + qi _ n B -q , 



2q ■ P tib -qriB-P ns-P 

where in the second equality we used that the proton momentum P is exactly along the ns direction, and in the 
last step used that q± is no bigger than 0(QX 2 ). (The directions nj and ns will always be chosen so that this is 
true, according to Eqs. (|116|) and (|Bla.[) . In other words, for events with small 1-jettiness, all the large momentum q 
transferred into the final state is collimated along nj and % with no 0(Q) momentum going in a third direction.) 
Thus to leading order in A the second argument of the beam functions in Eq. (|120[) is always just x. 

2. Jet Functions 

The vacuum collinear matrix elements in Eq. (|117p can be written in terms of jet functions in SCET [2l| . defined 
with transverse displacement of the jet in (67| by 

J q (u;k + +ul,ri= { ^- J |^ e ^V 2tr ^|| Xn ( 2/ -^) ( 5( w + n.7') ( 5 2 (^+^)x„(0)|o), (122) 
J g >fc+ + <x\ , /1) = ^ J ^-^ + y-/ 2 (Q\ Xn (y~l )6(w + n ■ V)S 2 (u J± + P x )|x»(0)|o) . 
Thus the vacuum collinear matrix elements in Eq. f|l 1 7[) can be expressed 

(0| xZ(0)6(Qjrj - nj-p n )6(hj-q + nj ■ V)S 2 (q ± +p ± + Pj.)x#" '(0) |0) (123) 

= M'JTJ + (q^+Pi-) 2 ^) , 

(0| X P n k A0)5(QjTj - nj-p n )S(nj-q + nj ■ V)5 2 (q ± +p ± + P±)xg? (0) |0) 

= ^^S kk 'M Sj r J + (q ± +p ± ) 2 ,^, 
where sj is defined in Eq. (|58p and (q± + p±) 2 = — (qj_ + P_l) 2 - 
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3. Hard and Soft Functions 
Using the above definitions of beam and jet functions, the hadronic tensor in Eq. (|117|) can be written as 



W"'(x,Q 2 ,n) = -2(2n)n B -qnj-q(QjQ B ) 2 [ d 2 p^—^— [ dTjdT B dT s J drf 5(n - tj - t b - r/ 
M J nj-n B J 



-T?) 



x S{Qjt s j ,Q B rf ,nj-n B ,n)J q {sjTj + (q± + p±) 2 1 n) 



(124) 



H q^v(q 2 , n J, n B)Bq(s B T B ,x,p 2 Ll fi) + H^(q 2 , n B , nj)B q (s b t b , x.p^jj) 



II' 



where the hard function is defined 



H lL v {(Pi -P2) 2 ,n a ,n b ) = Tr C Iqq ^(p 1 ,p2) r ^Ci' q q l/ (p 1 ,p2) 



4a, 



~ \4b 



and the soft function is defined 



S(kj,k B ,qj,q B ,fx) = — tr (0| [Y^ B Y nj ](0)6(kj - nj-p s j)5(k B - n B -p s B )[Y^Y nB }(0) |0) 



(125) 



(126) 



To write Eq. (|124p we used the equality of the quark and antiquark jet functions J q ^ q in QCD. 

a. Structure of the hard functions In Eq. (|125[) , the matching coefficients C, C in the hard function Eq. (|125[) for 
the vector and axial currents I — V : A take the form 



Vfqq(P^P2) = C Vfq((Pl-P2) )7x , C A f qq» (Pi , jh) = C A fq((pi ~ Pi) h±75 , 



(127) 



where 7^ is transverse to the directions 71^2 of the label momenta pi a- We have shown the index / for the quark 
flavor in the current explicitly. In Eq. (|124p these directions are nj B . The scalar coefficients Gy / qq ,C A f qq depend 
only on the symmetric Lorentz- invariant combination (pi — p 2 ) 2 . Using the momentum-conserving delta function in 
Eq. (j 11 Tj) . this combination takes the value (pi — P2) 2 = q 2 . Inserting Eq. (|127j) into Eq. (|125j) we obtain 



H q g fil ,{q 2 ,nj, n B ,fi) = C Ifq (q 2 ,n)C Ilflq {q 2 ,^) Tr(r£-^r£'-^ 
where = 7^ and = 7^75. Thus, there are two distinct traces to take in Eq. (|128|) : 

T VV, AA( „2 „_ „_ ,.\ _ n .T n B n ln 2 ..\ri__ l„2 ,.\ liv 



Hqq^v [q ,nj,n B ,fi) = - 

TjVA,AVf 2 \ 
Hqqul (Q ,n.i,n B ,fi) = -1 



4 

nj-n B 



Cv, A f q {q ,v>)C v ,Af'q(q ,fi)g± , 



(128) 



(129) 



C V , A fq{<l 2 ,V)C A y f>q(q 2 , fl) 



where <?j_ and e*f" are symmetric and antisymmetric tensors orthogonal to nj and n B given in Eq. (|B1|) . Hence, 
jjW, AA anc j j]VH,AV are symmetric and antisymmetric, respectively, under exchanging nj and n B . 

b. Structure of the soft function The soft function Eq. (I126|) depends on the momenta k Bt j projected onto 
the n B _j directions in the regions l-L B j, respectively. The shape of these regions in turn depends on the vectors 
q B _,j = us B) jn Bt j/2 in the definition of the 1-jettiness t\ in Eq. (|2~4"|). Indicating this dependence explicitly, we express 
the soft function Eq. (|126[) as 

S(kj,k B ,qj,q B ,fi) = ^-tr^2\(X„\ \Y%jY nB ](0) |0)| 2 d(kj - ]T 0(q B -k, - q.j -k^rtj -fc 4 ) 



S(k B - 6(qj-h - q B -ki)n B -k, 



iex s 



(130) 



Note that the soft function for DIS involves the square of one incoming and one outgoing Wilson line, and hence 
differs from that for e + e~ — > dijets that has two outgoing lines, and for pp — > L + 0— jets which has two incoming 
lines. We can relate Eq. (|130p to the usual hemisphere soft function for DIS by generalizing an argument given in 
[HI . Note that the Wilson lines Y n are invariant under rescaling of n (boost invariance) : 



Ypn B = P exp 



ig / ds (3n B - A s {j3n B s) 



Pcxp 
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,0 

/ dsn B -A s (n B s) 
J — 00 



Y ni 



(131) 
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and similarly for the lines extending from to +00, Yp nj = Y nj . Recall from Eq. ((50)) that 



so defining n'j = nj/Rj and n' B = ub/Rb we have (qB — Qj) ■ ki = ^ujbRb (n' B — n'j) ■ ki since lojRj = ojbRb- This 
implies that the same partitioning defined in Eq. (| 130[) can be expressed with 6{ ki) and 0(n'j ■ ki — n' B ■ ki). 

i'j = 2. Thus expressing Eq. (|130[) in terms of the rescaled vectors, n'j and n' B , we obtain 

S(kj,k B ,qj,q B ,fx) = Nc x jRb tr E \< X °\ [ y i' f ^](°) 1°) Vfe - E Oin'B-h-n'j-hyj-k 



\R] 



x S(^- - 0(n'j-ki - n' B -ki)n' B -ki 



R 



j Kb ^Hj Hb ' 



In the last equality we have expressed the fact that the expression in Eq. (|133l) is the same as the hemisphere soft 
function (up to the overall 1/{RjR b ) in front), with momentum arguments rescaled by Rj.b as indicated. Therefore 
from here on we will write all the T\ factorization theorems in terms of the DIS hemiphere soft function. Note that the 
vectors n'j B have been rescaled from nj lB such that they no longer have timelike components equal to 1 nor spacclike 
magnitudes equal to each other, and therefore do not partition the final states X s into hemispheres as viewed in the 
original nj^B frame of reference. However, the soft function in Eq. (|133|) depends on n'j B exactly like the hemisphere 
soft function depends on nj >B and depends on the dot product n'j-n' B , which is 2, making it equal to the hemisphere 
soft function. Physically, there exists a frame where n'j and n' B are back-to-back with equal time-like components, so 
that the partitioning in this frame gives hemispheres. 

In the 1-jettiness cross sections below, the soft function Eq. f|133[) will always be projected symmetrically onto a 
function of a single variable ks, following from Eq. (|124|) : 

Sbemi(ks,fJ>) = J dkgdk B S(ks - kg - ks)Sbemi(ks,ks ,fj) , (134) 

We will use the same name Shemi for the hemisphere soft function of two variables in Eq. (|133|) and its one-variable 
projection Eq. (|134[) . distinguishing them by the number of arguments we write. 



4- Final Form of Factorization Theorem for Hadronic Tensor 



Changing variables in the arguments of the beam, jet, and soft functions in Eq. (|124[) gives 

W0(x, <3 2 ,ti) = / dV-^- / dtjdt B dkj dkf sfn-^-tl- ^-±^1) J q (tj - (q ± + P± ) 2 , ft) (135) 
J nj-n B J V S J S B Qr J 



X Shemi (^-s ; k s , /i 



HlLv (q 2 ,nj,n B , n)B q (t B , x, p 2 ± , /x) + H" (q 2 ,n B ,nj, (i)Bq(tB,x,P±, fi) 



We have written the arguments of the jet and beam function in terms of dimension 2 variables the arguments 
of the soft function in terms of the total light-cone momentum k J s = nj-k J s in region J and kf = ns'kf in region J5, 
and have rewritten the transverse momentum arguments of the jet and beam functions in terms of two- vectors q^, 
instead of the the four- vectors q±,p±. The constant Qr is defined in Eq. (|5~4")) and s,/.s arc defined in Eq. and 
their special values for r^ b ' c are given in Table |TIJ 



5. Factorization Theorem for Cross Section 

In the cross section Eq. (|55|) . the hard function Eq. f|l 2B[) gets contracted with the leptonic tensor L^ v in Eq. (|84l) . 
The contraction of the leptonic tensor and the hard function can be performed using the tensor contractions in 
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Eqs. (|B2[) and (|B3[) . and can be expressed in terms of Born-level cross section and scalar hard coefficients as 

L l» f f ( x > Q^ nj^riB H ^^( q2 > n J> n B, M) = d ^g 2 -g<?(gj ; gB:Q 2 :M) , (136a) 
^L^^Q 2 )-^^^ 2 ,^,^^) = ^^Hq(qj,q B ,Q 2 ,f,), (136b) 

where the Born-level cross section is given by 

dop 47ra 2 m gj -k' q B -k + qj- k q B ■ k' 
dxdQ 2 x 2 s 2 Q 2 q.j-qB 

The hard coefficients of the quark and antiquark beam functions are 



(137) 



H q t q(qj, q B , Q 2 , fi) = ^ [{Cy fq Cvf'qL^Y fl + C* Afq C ' A f qLfffi) T r(gj, g S ) (CyfqCAf'qLXff, + CXfgCvf'qLfJf,)] , 

//' 

(138) 

where the relative minus signs for Hq come from the interchange of nj >B in Eq. (|136[) . The coefficients Cy.A = 
Cv,A(q 2 , A*) are functions of <? 2 and /x and the leptonic coefficients L g _ f = L g ^(Q 2 ) given in Eq. ([86)) . The coefficient 
r(qj,q B ) is given by 

/ \ qj-k' q B -k-qj-kq B -k' 

r(qj, q B ) = r, r— r 77 ■ ( 139 ) 

qj -k' q B -k + qj-kq B -k' 

Because the coefficient r is a function of scalar products of q B j and k and k' it becomes a function of y and Q 2 once 
q B ,j are specified as in Sec. lIII Al So, the hard coefficient H qq also is a function of y and Q 2 through the coefficient r. 
Contracting Eq. (|84[) with Eq. (|135|) then gives for the cross section Eq. 



dxJi&dn =J d2p± -irk 2 - £ - l f B - ^) ^ - + p^) 2 ^) (140) 

x [H q (qj,q B ,Q 2 ,fi)B q (t B ,x,p 2 L ,fi) + Hq(qj,q B ,Q 2 ,fi)Bg(t B ,x,p 2 L ,n)] , 
where we used the projection Eq. (|1 34[) of the soft function onto a single variable. 

D. Results for three versions of 1-jettiness r*,Ti, rf 

Now we will specialize the generic factorization theorem for 1-jettiness in Eq. (|140p to the specific cases 7""' ' c . The 
discussion will be most efficient if we begin with t\. 

1. 1-jettiness t\ 

The reference vectors q h B = xP and q b j = q + xP in Eq. ([32]) are used to define the 1-jettiness t\. In any frame q 
can be written as q = q b 7 — q B , so with respect to the directions n b Bl n b j, the transverse component q± = so that the 
argument of the jet function in Eq. (|140j) is (qj_ + p±) 2 = Pj_. Meanwhile, the coefficient r(qj,q B ) in Eq. f|139[) is 
given by 

r( g + «P,«P)= J ( ^ )a . (141) 

Note that r is a function only of y. So, the hard coefficients H q q in Eq. (| 138|) depend on y and Q 2 , and we define the 
hard coefficients for t\ by H q q (y,Q 2 ,fi) = H q t q(q + xP, P,Q 2 , //). Therefore, using Eq. f|140[) the final factorization 
theorem for rf is given by 



dx dQ 2 dr* da; dQ 



■ y dtjdt B dk s 5 (r b - ^ - ^ - 5 ho mi (fcs, A 1 ) 



(142) 

x / d 2 Pi J 9 (0-pi,Ai) H b q {y,Q 2 1 y)B q {t Bl x 1 p\ 1 ^) + H b {y 1 Q 2 1 ii)B q {t Bl x 1 p\ 1 ii) , 
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where we used Table HT1 to substitute for sj.b, Qp in Eq. (|140p . and where the Born-level cross section is given by 

da b _ 2na 2 em ^ _ y)2 + ^ (M3) 



dx dQ 2 Q 4 



2. 1-jettiness rf 



For the 1-jettiness t™ defined in Eq. (|29j). the minimization inside the sum over final state particles i groups particles 
with the reference vector to which they are closest. The reference vector q°j with which the jet particles are grouped 
is aligned with the jet momentum pj, so that the jet has zero transverse label momentum with respect to n". This 
direction n°j is the one which would minimize rf (to leading 0(A 2 ))with respect to variations of <?}. A jet with 
momentum pj — ujjnj/2 +pjj_ + k, where k is residual, has a mass m 2 — ujjnj-k +p 2 ±, so nj-k = (m 2 — pjj_)/u>j. 
The choice of nj which makes p 2 ± = minimizes nj-k (note that p\ < 0). 

The cross section for the rf distribution is given by Eq. (|140[) , with qp = q% = xP and qj = q% where qj is the 
vector q,j in Eq. (|29[) that minimizes rf . We will write q°j in terms of the vector q b = q + xP that was used to define 
the 1-jettiness t\. Now, the vector q h j has a direction n b j and magnitude w', given by 



Pre 1 



P T e 



P T h 



(144) 



expressed in the CM frame, where pr, Y are given by Eq. ([55)). With respect to n b j and np, the collinear fields in the 
jet function matrix elements still have nonzero transverse labels. Now, for each p-j , we rotate n b j to a vector n°j so 
that the transverse label with respect to n a j is zero. This requires that the total label momenta in the two coordinate 
systems be equal: 



(145) 



differs from n b j at most by a quantity of O(A). Now we express q in n",np coordinates. In Eq. (|140[) . the 



transverse label on the collinear fields in the jet function is q± 
which makes this quantity is zero, so q± = —p±. By using q 2 = 



-p±. The n a j,np coordinate system is defined as that 
— xys and np-q — yy/s, q is expressed as 



q = yVs- 



ij-np 



V 



\ np 
y^fs) 2 



Pi- 



Then Eq. KHD| takes the form 



da 



dx dQ 2 dif 



d 2 i 



del 



[ T~i ~ ^ ~ ^ - ^ ) Jq(tj,fl)S iien 



(146) 



(147) 



dx dQ 2 

X [H q (q a j,q%,Q 2 ,fi)B q (tB,x,pl,fi)+H q (q%q a B ,Q 2 , f i)Bq(t B ,x,p 2 ± ^)] , 

where we used Table [TT1 to substitute for sj t p and Qp. 

The generalized beam functions appearing here explicitly depend on p± . The vector n " appearing in q°j implicitly 
depends on p±. Now, n a j differs from n b j (which is independent of p±) by a quantity of order A. Here we can expand 
the hard and soft functions and the Born cross section around n°j = n b j + 0(A) and drop the power corrections in A. 
This makes everything in Eq. (|147j) independent of p± except for the generalized beam function. The integral over 
p± then turns the generalized beam function into the ordinary beam function Eq. (|119D . Thus the final factorization 
theorem for the cross section is 



da 



da b 



dx dQ 2 dT c { dx dQ 2 



dtjdtpdks 5 



t_B_ 

r>2 



Q 



Jq{tj, /J,)Shemi(ks, A*) 



H b { Vl Q 2 , ti)BM B ,x, fi) + H b { Vl Q 2 ,v)Bq(t B ,x, M ) 



(148) 



where the Born cross section is given by Eq. (|143[) and H b ^(y,Q 2 ) is given by Eq. f|l 38[) with r in Eq. (| 141[) . The 
hard and soft functions in Eq. (|148j) are the same as those in Eq. (|142|) for t\ . 

Eq. (|148[) differs from Eq. (|142p in that the jet and beam functions are no longer convolved together in the transverse 
momentum p±. The 1-jettiness is proportional to the invariant mass of the jet, while t\ measures the projection 
q b j ' Pj onto the fixed axis q j = q + xP. The emission of ISR with transverse momentum p±, causing a shift in the 
jet momentum by the same amount due to momentum conservation, will not change the mass of the jet, but it will 
change the projection of the jet momentum onto the q b axis. Thus r" involves no convolution over p±, while t\ does. 
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3. 1-jettiness rf 

For the 1-jettiness rf, the directions nj and ng are along the electron and proton directions, respectively: 

nj = n e , tlb — tip 7 (149) 

where ri e ,_p are the light-cone directions of qj = k = u e n e /2 and qp = P = ujpnp/2. In the CM frame, n e .p arc 
back-to-back, n e = n z and np — h z . In this frame, q is given by 

9 = - xyV^Y + q± , (150) 



where ^ = Q-y'l — ynj_. Let us consider for a moment the power counting of the argument of the jet function in 
Eq. (|140[) with rf. The requirement that ~ QA requires that 1 — y ~ A 2 . This is ensured by measuring rf to be 
C(A 2 ). The argument of the jet function (call it m 2 ) in the factorization theorem Eq. (|140p for rf is 



m 2 = tj - (1 - y)Q 2 - 2Q v /T^n ± -p ± - pi . (151) 

Now, the jet function will be proportional to a theta function 9{m 2 j ), requiring to 2 > 0. Measuring rf to be of order 
A 2 and therefore forcing tj to be of order Q 2 A 2 then enforces that 1 — y ~ A 2 . Then, we can set y = 1 to leading 
order everywhere in Eq. (|140[) except in the argument of the jet function. In terms of x, using the relation xys = Q 2 , 
requiring y < 1 is equivalent to requiring x > Q 2 /s, which sets a lower bound on x. 

The normalization constants sjb,Qr in Eq. (|140|) are given for rf in Table HT1 The Born- level cross section and 
the coefficient r(qj,qp) in the hard coefficients reduce to 

daZ 27ra 2 „, , 

- em r(fc,P)=l, (152) 



da; dQ 2 Q 4 

where we see the Born cross section is now Eq. (|143[) in the limit y — > 1. This happens because the expression 
Eq. (|137j) is evaluated with qj = fc, which is the actual jet direction only near y — > 1. The hard coefficient is now 
independent of x, y and depends only on Q 2 : Hq^(Q 2 7 fi) = H q ^{k, P,Q 2 , fx). From Eq. (|140[) the final factorization 
theorem for the cross section in with rf is then given by 

dxdQ'dTf = d^ 2 J ^ J 6 { T ' - ^ - - ^§Q) - P ^ 2 ' $ ^hemi{ks , fl>) 

H c q (Q 2 ,n)B q (t B , x,p 2 x ,n) + H c q (Q 2 , fi)Bg(t B ,x,pl, fi)] . (153) 

This is like the t\ cross section Eq. (|142[) in that the jet and beam functions are convolved in the transverse momentum 
px of ISR, but in this case the momentum transfer q itself has a nonzero transverse momentum with respect to the 
light-cone directions n e> p. This will make the evaluation of the pj_ integral considerably more involved than in the 
t\ cross section Eq. (|142j) . 



VI. FIXED-ORDER PREDICTIONS AT 0(a s ) 

In this Section we evaluate to 0(a s ) the predictions 
of the factorization theorems for the cross sections differ- 
ential in the different versions of 1-jettiness in Eq. (|142[) 
for tI, and Eq. (fT48]) for if and Eq. (fT53| for rf. These 
formulas correctly predict the singular terms at small n 
in the fixed-order differential cross section, although they 
have to be resummed to all orders in a s to accurately pre- 
dict the behavior at small T\ . We will do this in the next 
Section. Also, for the predictions to be correct for large 
T\ , they would have to be matched onto 0{a s ) and 0(a 2 ) 
fixed-order full QCD calculations, an analysis we defer 
to future work. Nevertheless we can estimate the size of 
these matching corrections by comparing our predictions 



integrated up to large T\ to the known total QCD cross 
section a(x, Q 2 ) at 0(a s ), which we will do in Sec. IVIIII 

A. Hard Function 

At 0(a s ), the matching coefficients Cyfq, CAfq for the 
vector and axial currents Eq. (|99[) that appear in the hard 
coefficient in Eq. (|138[) are equal and diagonal in flavor, 
and were calculated in [3(1 [91j : 

C vfq {q 2 ) = C Afq (q 2 ) - 8 fq C{q 2 ) , (154) 

C(q-) = l + ^^(^n 2 J^-3lnJ^-8 + ^) 
47r \ — q z —q z 6 / 



2G 



For DIS recall q 2 = —Q 2 . Then, the hard coefficients 
H q q in the cross section Eq. ([140[) are given to 0(a s ) by 

Hg,q(qj,qB,Q 2 ,(i>) = H(Q 2 ,fj,)L qi q(qj,q B ,Q 2 ) , (155) 

where we have defined the universal SCET 2-quark hard 
coefficient, 



H(Q 2 ,fl)= \C(q 2 ,y) 
= 1 



(156) 



a 3 (fx)C F ( 2 M" oi„ 



2tt 



- 1 ">- 31 "^- 8+ y 



and the factor containing the components of the leptonic 
tensor Eq. (|84|) . 

£^r(« Jjte )(Z^ + 0. 



L q ,q(qj,qB,Q ) = L gqq , 
where r(gj, g^) wa s defined in Eq. (|139p 



(157) 



For the r"' & cross sections Eqs. (|142|) and (| 148[) . = 
xP and q.j = q + xP, so that r(qj,qs) is given by 
Eq. (fT4l~]) . Then the leptonic factor L q<q in Eq. (fT57) 

becomes 



L q Aqj,q B ,Q 2 ) 



(158) 
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^999 ^ 



=F 



2Q q v q v e 



y( 2 ~y) (t va 

(l-y) 2 + l ( 

(v q + a 2 q )(v 2 e 



jAV\ 

-al) 



(l + m 2 z /Q 2 ) 2 



' l + m 2 z /Q 2 

2y(2-y) a q a e [Q q (l + m 2 z /Q 2 ) - 2v q v e 
(l-y) 2 + l (l+m 2 z /Q 2 ) 2 

2. T\ cross section 



For the rf cross section Eq. (|153[) . q,j = k and qs = P, 
the electron and proton momenta, respectively. Then 
r(k,P) = 1 in Eq. (|139p . and the leptonic factor L q>q in 



Eq. (|157[) becomes 



L c q Jk,P,Q 2 )=L v ' 



tAA itVA 
991 ' 911 ' V £99 

2Q q (v q v e ± a q a e ) 



£99 ^ £99 



) (159) 



l + m 2 z /Q 2 
(v q + a 2 )(v 2 + al) ± 4v q a q v e a e 
(l + m|/Q 2 ) 2 



B. Soft Function 

The soft function Shcmi(fcs: A*) th & t appears in the 
cross sections Eqs. (|142j) , (|148[) . and (|153p is given by 
Eqs. (|T35f and (fTM)) . For e+e" -> dijets, is known 

at 0(a s ) HI and 0(a 2 ) [13 E0, [To3l ■ At ^oop order 
the dijet soft function is the same for DIS. Beginning at 
2-loop order, the finite part of the soft function for DIS 
could possibly differ due to switching incoming and out- 
going Wilson lines, but the anomalous dimensions and 
thus the logs are the same. 

To 0(a s ), the soft function Eq. f|l 33[) takes the form: 



(160) 



where 
S {1 \k s ^) 



47T 



-S(k s 



0(fc s )ln(fc s /M) 



k s /fi 



(161) 



and the projection Eq. (|134l) is then given to 0(a s ) by 

S hcmi (k s ,n) = S(k s ) + 2S {1 Hk Sl [i). (162) 

It has previously been observed that the sizes Rj t B of the 
regions Hj,b to which soft radiation is confine d enter th e 
arguments of the logs in the soft function [gj llOll Il02j | , 
which is due to chan ging the effective scale at which the 
soft modes live Il03l. 



C. Jet Function 



The jet function Eq. (|T2"2"j) is given to 0(a s ) by |104ll0 



Jq,q(t, fi) = 5{t) + 



-in 



3 


r e(t) - 


4 


\d(t)ln(t/v 2 ) 


V 2 




+ — 

+ ^ 





(163) 



It is in fact known to two-loop order [105] and its anomalous dimension to three loops [3E 
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D. Beam Functions 

1. Generalized Beam Functions 
The generalized beam functions in Eq. (|118|) can be matched onto ordinary PDFs, defined in SCET as [22 
f g {cj'/P-,n) = 9{uj') (P n (P-)\ X„(0)|[£(o/ - n ■ V) Xn (0)] \P n (P~)) , 
f q (w'/P-,iJL) = 6(u/) (P n (P-)\tv |xn(0)[*(a/ - n • P)*n(0)] |P„(P")) ■ 



The matching result is |6j||94(: 



A 2 A 2 x 



(164a) 
(164b) 

(165) 



where i,j = q,q,g. This expansion is valid for perturbative beam radiation satisfying t, k^ 3> Aq Cd . At tree level, 



4 u, (M,ki,/,) = (l/7r)^5(t)5(l - z)«(ki), leading to B\ 0, (t, x, k\, fi) = (l/7r)6(t)S(kl)fi(x,fJ,). 

To 0(a s ), the nonzero matching coefficients in the generalized quark beam function were computed in [681 19~4|. and 
we use the results from 1941: 



X qq (t, z, kl, n) = -6(t)6(l - z)8(k 2 ± ) + 



2tt 2 



0(t) 



l qg (t,z,k\,ij) 



+ S(t)S(kl) 



P qq (z)- -6(1- z) 





\0(t)\n(t/n 2 )] 




[ t/^ \ 


s(kl 


0—z)t\ 
z / 



6(1 z)6(k 2 ± ) 



(166a) 



Jl - z)ln(l - z) 



1 - z 
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(1 + z 2 )- — 5(1 - z) + 6(1 - z) il-z 

6 V 1 — 2 



P qg (z)6(k 2 ± - 



(l-z)t 



8(t)5(k 2 



P qg (z)\n +26(l-z)z(l-z) 



where P qqqg arc the q — > qg and g — > gg splitting functions, 



0(1 -z) 



1 - z 



(1 + z 2 ) + -6(1 - z) 



il-z) 



1 + z 



- 1 



1-2! 



P qg (z) = 9(1 - z)[(l ~ zf + z 2 



They appear in the anomalous dimensions of the PDFs, which to all orders obey 

^Me,M) = E/f^(fM)M' ! M). 

At 0(a s ) the anomalous dimensions for the quark PDF are 



7« 0,/f) = 0(z)P qq (z). 



Too ( Z ! W (*) • 



(166b) 

(167a) 
(167b) 

(168) 
(169) 



2. Ordinary Beam Functions 

The ordinary beam functions Eq. (|119[) satisfy the matching condition [67], Il06l Il07j 

-A 2 



Bi(t,x,n) =yZ j y 2 ^ 1 



( t> | >A1 )/ j(e , M )[i + o(^) 



(170) 
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where at tree level 2^\t, z, fi) = Sij6(t)5(l — z), leading to Bf (t,x, fi) = 5(t)fi(x,fi). To 0{a s ), the matching 



(0)/ 



coefficients in the quark beam function are given by integrating Eq. (I166[) over k^. [671 . llO 



l qq (t,z,fi) = 8(t)S(l -z) + 



X qg (t,Z,fJ,) 
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2tt 



2tt 
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P qg (z) + S(t) 



1 - z 



P qg {z)\n +20(l-z)z(l-z) 



(171a) 
(171b) 



E. Dijet Cross Section 



We can now form the SCET predictions for the T\ cross section Eq. (|140|) to 0(a s ) by plugging in the 0(a s ) 
expressions for the hard function given by Eqs. (|155j) . (|156|) . and (|157p . the soft function given by Eqs. (|160j) and 
(|161[) . the jet function given by Eq. f|163[) . and the generalized beam function given by Eqs. (|165[) and (|166p . It is 
convenient to express the result in terms of the cumulant t± distribution, defined by 



a c (x,Q 2 ,T X ) = — [ dr[ 
CT o Jo 



da 



dx dQ 2 dr[ 



(172) 



where <jq is the Born cross section defined in Eq. (|13T|) . We will give here the results for the rf' b cumulants at 0(a s ). 
The more complicated results for the rf and generic t± cumulants are given in App. [G] 



1. T\ cross section 



Plugging in the 0(a s ) results for the hard function given by Eqs. (|155l) . (|156[) . and (| 1 58|) . the soft function given by 
Eqs. (|160[) and (|16ip with sj t B and Qr given in Table [TTl the jet function given by Eq. (|163[) . and the ordinary beam 
function given by Eq. (|170|) into the r" cross section Eq. (|148[) . we obtain for the r" cumulant given by Eq. (j!72[) in 
the CM frame: 



ct c (x,Q 2 ,t?)=0(t?) / -[L a q (x,Q 2 )f q (x/z,fi) + L a q (x,Q 2 )f q (x/z,fi)] 

J x z 

x j<5(l-z) 

, Ol s Cf 



(173) 



l_^( 9 + ^! +6 l nrf+ 41n 2 rr 
47r V 3 



2tt 
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£i(1 - z){\ + z 1 ) + (9(1 - z)\ l-z 



9(T?)[L a q (x, Q 2 ) + L1(x, Q 2 )] 1 1 ^f g (x/z, M )| ln( 



l + z z 
1 - z 



In z 



In 



Q 2 rf l-z 



P M (z) + 20(1 -*)*(!-*) 



The factorization scale (i still appears on the right-hand side of the equation, though the cross section is in fact 
independent of fi. The /x-dependence in the PDFs on the first line is cancelled by the ^-dependence in the logs 
multiplying the splitting functions on the third and final lines to 0(a s ). The residual /Lt-depcndcncc is 0(a 2 s ) and 
would be cancelled by the higher-order corrections. 



2. t\ cross section 

The t\ cross section is nearly identical to the rf cross section except for the presence of the p^-dependent generalized 
beam function in Eq. (|142[) instead of the ordinary beam function. The effect of the nontrivial p^-dependent terms 
in the generalized beam function Eq. (|165[) is simply to multiply the arguments of the /i-dcpcndcnt logs in Eq. (|173[) 



29 



by z, giving the simple relation 



d I dz ( 

(x,Q 2 ,T b 1 )=a c (x,Q 2 ,T?) +ff( 7 f)_i / —\nz \C F [L a q (Q 2 )f q (x/ z, ll) + L a q (Q 2 )f q (x/ z^)]P qq (z) 

+Tf [Lg(Q 2 ) + L°j(Q 2 )] P qg (z)f g (x/ Z, , 



(174) 



In the Appendix we give the 0(a s ) rf cross section. In the next section we resum the large logarithms of t x 
that appear in these fixed-order expansions to all orders in a s to NNLL accuracy. 



a,t>,c 



VII. RESUMMED PREDICTIONS FOR n 
CROSS SECTIONS 

The fixed-order predictions for the t\ cross sections 
presented in the previous Section contain logarithms of 
Ti which grow large in the limit n — > and must be 
resummcd to all orders in a s to yield accurate predic- 
tions for small T\. In this Section we use the factoriza- 
tion theorem Eq. (|140[) for the t\ 1-jettiness cross section 
and its specialized cases Eqs. (|142p . t|148|) . and (|153p for 

, rf , rf to predict the cross sections differential in these 
variables to next-to-next-to-leading logarithmic (NNLL) 
accuracy, estimate the perturbative uncertainty by ap- 
propriate scale variations, and discuss power corrections 
due to hadronization, including their universality and im- 
pact in the tail and peak regions. 



Perturbative Resummation to NNLL 



where the anomalous dimension takes the form 



The hard, jet, beam, and soft functions in Eq. (|140[) 
obey renormalization group (RG) evolution equations 
whose solutions allow us to resum large logarithms of 
ratios of the separated hard, jet, beam, and soft scales. 
These solutions allow us to express any of these functions 
G = {H, J, B 7 S} at one scale fx which contains logs of ii 
over some scale Qg in terms of the function evaluated at 
a different scale lig ~ Qg where the logs are small. 

The hard function H(Q 2 ,/j,) obeys the RG equation 



At— H(Q 2 ,l 1 )= 1h (ll)H(Q 2 ,li), (175) 
aii 

where the anomalous dimension j B has the form 

Q2 

jh(h) = T H [a s (Li)} In— 5- + j H [a s (Li)] , (176) 



2r^ usp and a non-cusp piece 



with a cusp piece r^[o; s ] 

1h[o-s\ (which is conventionally denoted by the same 
symbol as the total anomalous dimension). Their expan- 
sions in a s are given below in Eq. (|182[) and Eqs. (|D28|) 
and (|D29|) . Similarly the jet and beam functions which 
are both functions of a dimcnsion-2 variable t obey RG 
equations of the form 



H—G(t,n) 
d\x 



dt' lG (t-t',Li)G(t',Li) 7 (177) 



la{t,Li) = T G [a s (LL)]\ 



t/v 2 



-y G [a s {Li)]6(t) , 



(178) 

where here G = {J, B}, and the plus distribution is de- 
fined in App. O The cusp pieces Fj tB = — 2r<? and 
non-cusp pieces ~/j tB of the jet and beam anomalous di- 
mensions are given in Eqs. (|D28|) and (|D30|) . The beam 
function also depends on x and the generalized beam 
function also depends on , but they do not change the 
structure of the RG equation Eq. (|177p . 

Finally, the soft function in Eq. f|L33[) obeys the RG 
equation 



A*-J-<S'hemi(fc/, k Bl Ll) = 



dkjdkg 



(179) 



x ls{k,j - k'j,k B - k' B ,(i)Shsmi(kj, k' B ,Li) , 

where the anomalous dimension factorizes into the form 

ls{kj,k B ,Li) = 7s(fcj, fi)5(k B ) +Js(k B , fi)S(kj), (180) 

which is required by u - independence of the total cross 
section Eq. (|140p |l08| . Each piece 7s(fc,/x) takes the 
form 



ls{k,Li) = 2T s [a s (Li)} 



0{k/Li) 
k/iL 



^ls[a s {ii)}5{k) , 
, (181) 

where Y$ — r cusp , and the non-cusp piece is given by 

7s = -7ff/ 2 - U- 

The cusp and non-cusp pieces of the anomalous di- 
mensions of all the functions above all have perturbative 
expansions in a s : 



OC , n 

* / {X \ n + 1 



n=0 



IMS) 



n=0 



(182) 

which defines the coefficients T G: ~f G . Furthermore, the 
cusp pieces of the anomalous dimension are proportional 
to the same cusp anomalous dimension r* usp [a s ], whose 
perturbative expansion along with the non-cusp anoma- 
lous dimensions are given in App. |DJ The explicit solu- 
tions to the RG equations for the hard, jet, beam, and 
soft functions individually are given in App. [Dl 

The solutions of the RG equations Eqs. (|175[) . (|177p . 
and (|179p allow us to express the hard, jet, beam, and 
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soft functions at any scale /i in terms of their values at 
different scales ^h.,j,b.s where logarithms of Hg/Qg m 
their perturbativc expansions are small. There are a 
different conventional ways in the literature to express 
the resummed cross section in terms of the solutions for 
hard, jet, beam, and soft functions to the RG equations. 
One method [Til . [76j performs the exact inverse trans- 
form back from Fourier space, and carries out analytically 
the convolution of all the evolution factors and the fixed- 
order functions for the T\ factorization theorem Eq. (|140[) 
in momentum space. In this section we use this method 
and formalism, relegating some of the required formu- 
las to App. |U We give an alternative equivalent form of 
the resummed cross sections in App. [FJ using a method 
[38l H3] that first Laplace transforms the cross section and 
writes certain corrections as derivative operators before 
transforming back to momentum space. This avoids tak- 
ing explicit convolutions of the evolution factors and the 
fixed-order functions. If one carries out these derivatives 
analytically then the final results from the two formalisms 
are identical. 

In this section we give just the final results for the RG 
improved cross sections for the 1-jettinesses r°' b,c using 
the formalism of [l4|, [Zll • We will express the results in 
terms of the cumulant r°' b,c distributions: 

a c (x,Q 2 ,n) = - r dr[ % (183) 

(To Jo dx "Q " r l 

where we note that a c is dimcnsionlcss due to the division 
by a . The differential cross section can be obtained 
by taking the derivative of a c (x, Q 2 , t\) with respect to 
T\. Care must be exercised in this procedure because a c 
also depends on t\ dependent jet/beam and soft scales in 



the factorization theorem Eq. f|140[) . a c (x, Q 2 , 7*1, jitj(fi)). 
The appropriate procedure is to use, for e — > 0, 

da _ Qcjx, Q 2 ,T 1 + e,n. i (T 1 )) - a c (x, Q 2 , Tj-e, /^(n)) 

(184) 

where dajdr\ = (1 / ao)da / dr^ . See Ref. [l4[ for further 
discussion of this point. 



1. r1' b cross sections 

The cross section in Eq. (|140j) is expressed as a convolu- 
tion of jet, beam, and soft functions in momentum space. 
To rcsum the large logs, each function is RG evolved from 
a scale where the logs are small, an operation which is 
in the form of a convolution of an RG evolution ker- 
nel and the fixed order function as in Eqs. (|D5[) and 
(|D14|> . The evolution kernels Uj tB ,S in Eqs. (|D5|) and 
(|D15[) arc plus distributions, and each fixed order func- 
tion can also be written as a sum of plus distributions 
as in App. IE II Thus, the resummed cross section con- 
tains numerous convolutions of plus distributions £„, 
which we can compute by repeatedly applying the plus 
distribution convolution identity in Eq. (|E8[) . The cross 
section then gets written as a resummation factor times 
sums of products of coefficients called V in App. IE 2l and 
J„, /9<mw and S n in App. [ET] 

The resummed t" and t\ cross sections in Eqs. (|142[) 
and p48[) . obtained from RG evolution of the hard, jet, 
beam, and soft functions, are given by: 



<T c (x,Q ,t) 



r(i + n) \hh 



r dz 

Y,L a q {x,Q 2 ) / -f j (x/z,fx B )[W gj (z,T)+AW q , 

J x z 



(185a) 



Wqj(z,T) = H(Q 2 ,fj, H ) 2J a s (nj), 



AW qj (z) 



n 3 = — 1 



M 2 , 



T'U 

ri2 



a s (fj,B), z, 



tQ 2 



s. 



"3 



Ms 



ni+n2 + l ti+U3 + l 

i 1= -\ e 2 =-i 



for rf 

^1 [5 3q C F P qq (z) + S ]g T F P qg (z)} In z for r\ 



(185b) 
(185c) 



Here j sums over quark flavors and gluons, and the +(q O q) includes the term where the virtual gauge boson couples 
to an antiquark. In Eq. (|185a|) the exponent is a resummation factor that resums the large logs and the terms W q j 
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and AWgj are fixed-order factors which do not contain large logs. The evolution kernels K, and SI are given by 

K, = IC{hh,Hj,Hb,Hs,h) = K H (n H ,n) +Kj(hj,h) + K b (pb,h) + 2K s (hs,v) (186a) 
SI = tt(nj,HB,Hs,n) = + Vb{^b,h) + 2r] S (ns,^) , (186b) 

where the individual evolution kernels Kh, Kj = Kb, Ks, tjj = i]b, and ijs are given below in Eqs. (|D3[) . (|D5[) . 
and (|D15I) . Note that JC and SI are indpendent of /z because the [i dependence cancels between the various Ki and 
r/i factors in the sums. Their expressions to NNLL accuracy are given in Eq. (|D26|) App.fDj The coefficients J„, I qq , 
II s , S„ in Eq. (fl85bl) are given in Eq. ([E7| . The constants V fc mn and V fe ™(Q) are given in App.[E2] 

Note that in the resummed cross section Eq. (|185j) the coefficients J n , 1$ , and S n are functions of logarithms of 
their last argument as shown in Eq. (|E7|) and the hard function also depends on the logarithm ln(Q 2 / M?r)- The logs 
in these fixed-order factors can be minimized by choosing the canonical scales 

fi H = Q, (j, j = (j, B = Q\j Ti' b , /i S =Qr"' 6 . (187) 

Large logs of ratios of these scales to the arbitrary factorization scale /x are then resummed to all orders in a s in the 
evolution kernels K, and SI. The choices in Eq. (|187[) arc appropriate in the tail region of the distribution where T\ 
is not too close to zero and not too large so that the logs of T\ are still large enough to dominate non-log terms and 
need to be resummed. Near T\ ~ and t% ~ 1, we will need to make more sophisticated choices for the scales, which 
we will discuss in Sec. I VII CI 



2. T\ cross section 



The resummed rf cross section obtained from RG evolution of the hard, jet, beam, and soft functions in Eq. (| 153[) 
is given by 



cj c c {x,Q 2 ,t) 



e K- lE n / Q y g (/« g ,M) / ( T _ 1 + y ) a .Q2 yB(M g .M) / (j-l + y) Q2y »Q>**) { (r-l + y) ^Q^^**) 



r(i + n) V/W 



Mj 



MS 



1 



x L c q (x, Q 2 ) / —fj(x/z) [W%(z, T-l + y) + AW^(z, r - 1 + y)\ + (q O q) , 

J x z 

where W^- and AW^ are the fixed-order terms from jet, beam, and soft functions: 



(188) 
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[S jq C F P qq (z) + 8 jg T F P qg (z)] 



11 <i 



>(t) In 



r(l-A-) 



(i - y)x 



H(-l - St) 



(i - y)x 



SI \{l-y)X 



-Sl,-Sl;l-Sl;- 



(i - y)x 



(189a) 



(189b) 



-e(-r)- 



sin(TrO) 



where X = x(l — z)/(x + z — xz). Note that the r in W qj (z,T) and AW qj (z,r) gets shifted by 1 - y in Eq. (|188|) . 
H(n) is the harmonic number and 2F\(a,b;c; z) is the hypergcomctric function. The additional more complicated 
terms in AW^j- are due to the nontrivial p± integral in Eq. (| 1531) which convolves the terms in the generalized beam 
function with nontrivial dependence with the dependence of the jet function on (q_L + P_l) 2 , with q± 7^ when 
y < 1. Note that the term on the last line of Eq. (|189b|) contributes below if = 1 — y when plugged into Eq. (|188|) . 
but that the size of the correction in this region is very small. 

The second arguments of J n , 1$ , and S n in Eq. (|189a[) show that the canonical scales should be chosen to minimize 
the logs of the agruments, which are the fixed-order terms in the hard, jet, beam, and soft functions. 



Hh = Q, M./ = Q\Jt{ -1 + 2/, fJ-B = Qy x{t{ - 1 + y) , fi s = \fxQ{r{ - 1 + y) 



(190) 
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7[a»] 


/?[«.] 


{H,J,B,S}[a s ] 


LL 


Ois 


1 


a s 


1 


NLL 




a s 


a? 


1 


NNLL 


<4 


a? 









r[« s ] 


7[«s] 




{fl-,J,B,5}[Q s ] 


LL 




1 




1 


NLL' 


ci 


a s 


a? 


a s 


NNLL' 


*l 


a? 




a? 



TABLE III: Orders of logarithmic accuracy and required order of cusp (r) and non-cusp (7) anomalous dimensions, beta 
function /3, and fixed-order hard, jet, beam, and soft matching coefficients H,J,B,S. The "primed" counting includes the 
fixed-order coefficients to one higher order in a s . 



Here the whole cross section is shifted to the right by an amount 1 — y due to the nonzero and choice of axes for 
the jet and beam scales are separated by a factor ^fx due to the different normalization of the qs 



rf. Unlike t 



a,b 



reference vector in the definition of rf. For rf, q c B = P while for r"' h , q B ,b = xP. The soft scale is also rescaled by 
\fx. We will discuss below a more sophisticated choice of scales than Eq. (|190p that give rise to proper behavior in 
the limits r f — > 1 — y and rf ~ 1. 



3. Logarithms included in our LL, NLL, and NNLL results 

It is worth briefly discussing the logarithmic accuracy of our resummed results. Although this discussion is standard 
in the literature, sometimes the same notation is used for different levels of resummed precision, so it is worth being 
specific about our notation. The order in a s to which the anomalous dimensions, running coupling, and fixed-order 
hard, jet, beam, and soft functions are known determines the accuracy to which the logarithms of r in cross section 
are resummed. It is most straightforward to count the number of logs thus resummed in the Laplace transform of the 
cross section (cquivalently we could consider the Fourier transform to position space), 



a{x, 



dre -ur da 

dx dQ 2 dr 



(191) 



The fixed-order expansion of cr(x, Q 2 , v) takes the form, 



o(x,Q 2 ,v) = 1 + — -{c 12 L 2 + cuL + c w + d\{v)) 

47T 



(192) 



(c 2 aL a + C23L 3 + C22L 2 + c 2 \L + c 2 o + di{v)) 



(f 1 ) ( c 36i 6 + C35L 5 + c 34 L 4 + C33L 3 + c 32 L 2 + c 31 L + c 30 + <fe(i/)) +■■■ 



where L = log v. The largest log at each order in a s is a"L 2 ™. Our results in Eqs. (|185j) and (|188|) . once Laplace 
transformed, reorganize and resum the logarithms into the form: 



a(x, Q ,l>) = exp 



-^{C 12 L 2 + CuL + C w ) 



5 



V47T 



C22L 2 + C21L + C20) 



(C34L + C33L 3 + C32L 2 + C31L + C30) + 



(193) 



+ d(x,Q 2 ,v) , 



where the largest log at each order in the exponent is 
a™L n . The coefficients c nm , C nm , and d n {u) are func- 
tions of x and Q 2 . The function d(x,Q 2 ,v) contains 
terms d n (v) and is a non-singular function of v that van- 
ishes as v — > 00 (r — S- 0). Transforming Eqs. (|192j) and 



(11931) back to momentum space using 

do /""y-tzoo 

dxdQ 2 dr J-y-ico 27r« " ' ' ^ ^ 

where 7 lies to the right of all singularities of the inte- 
grand in the complex plane, defines the accuracy to which 
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logs of t in the cross section and its cumulant <t c (x,Q 2 ,t) where 
are resummed. 

Our main results in Eqs. (|185|) and (|188[) resum singu- C\ q {z) = Cp 

lar logarithmic terms a™ ln m r, but not the terms in the 
non-singular d(x,Q 2 , r) (inverse transform of d). The 



d(x, Q ,t) must either be calculated by comparing a full 
QCD perturbation theory calculation with the resummed 
result and determining the difference order by order in 
a s , or by determining the next-to-singular infinite tow- 
ers of logarithmic terms in d(x,Q 2 , r) by carrying out a 
factorization and resummation analysis in SCET at sub- 
leading power. 

Fixed-order perturbation theory sums the series in 
Eq. (|192j) row-by-row, order-by-order in a s . When the 
logs are large this expansion is not well behaved. Re- 
summed perturbation theory instead sums the exponent 
in Eq. (|193|) column-by-column, in a modified power ex- 
pansion that counts lnr ~ l/a s when the logs are large. 
Everything in the first column of Eq. (|193p is 0(l/a s ) 
[leading log (LL)], the second 0(1) [next-to- leading- log 
(NLL)], the third 0(a s ) (NNLL), etc. in this counting. 
Each order of logarithmic accuracy is achieved by cal- 
culating the cusp and non-cusp anomalous dimensions, 
running coupling, and fixed-order hard, jet, beam, and 
soft functions to the orders given in Table 11111 Another 
common counting used in the literature (eg. [3 IzE HH) 
is the primed counting which accounts for the fixed-order 
matching coefficients H, J, B, S at one higher order than 
in the unprimed counting. This primed counting is par- 
ticularly useful when one also requires predictions for 
transition regions where the size of the logarithmic and 
non-logarthmic d n (r) terms are comparable. Since in this 
paper we have not considered the nonsingular terms we 
adopt the unprimed counting (LL, NLL, NNLL) through- 
out. 



B. Comparison to NLL DIS Thrust tq 



1 



2£ 1 (l-z)-(l + z)ln(l-z) 
2 9 , 



5(1 



As discussed above in Sec. IIIIBI there are several ver- 
sions of DIS thrust discussed in the literature. Here we 
consider the version of thrust called tq in (l5| . to which 
we can directly compare our results for rf , since as shown 
above in Eq. (j6"3"]l of Scc. lIIIDl thev are one and the same. 
We will see that at NLL accuracy the result in [l5| for 
tq is equivalent to our result in Eq. (|185l) for t\ for the 
particular scale choices Eq. (|187[) in the SCET cross sec- 
tion. The NLL resummed cross section given in (l5| in 
the MS scheme is, in our notation, 



(195) 



<*s{Q) f 1 dz 



2 71 



—Ci q (z)f q {x/z,Q) 



2- 



' 7T 

~3 ' 2 

Ci fl (z) = Tp [P qg (z) ln(l - z) + 2z(l - z)] 



(196) 



and complete expressions for the resummation constants 
<7i,2 can be found in Ref. [l5[. They have fixed-order 
expansions in a s = a s (Q) given by 



ln 3 r 



In 2 r 



(197) 



gi lnr = G^-r 1 In r - 
Ztt 

fj2 = -Gn lnr 
Ztt 

with the coefficients 



G12 = -2G F , Gn = 3G F , G 23 - 2irf3 G 12 , (198) 

°«--5^' + (t-£) c * c ' + 5 c '"'- 



Note that the cross section Eq. (|195[) includes only the 
photon contribution for the intermediate gauge boson 
mediating the scattering, so for the comparison we spe- 
cialize our results to this case. 

By comparing to the resummed cross section in 
Eq. ([T85]) . we find that the result of given in Eq. (fl95|) 
is equivalent to the SCET photon induced cross section at 
NLL order with the following fixed choices for the scales 
in the evolution factors: 

Hh = Q, A* = A* J = A*s = QVn , A«s = Qn ■ (199) 

Thus the two results agree at NLL order. 

We note that in the fixed-order coefficient in Eq. (11951) , 
the choice /i = Q\fr has been made in the tree level term, 
but the 0(a s ) terms have been evaluated at /i = Q. 
In the SCET result Eq. (|Tg5) (or Eq. (|F7j)) pieces of 
the a s (fi) terms are evaluated at \i = Q, \i = Q\/t, or 
/i = Qt according to whether they come from the hard, 
beam/jet, or soft functions, while the PDFs are evaluated 
at fi = Q\fr. The difference between the SCET result 
and Eq. (|185|) is NNLL, since the error is ~ a 2 hiT in the 
fixed-order coefficient. In our counting taking the correct 
scales for a s (/j,) is required for NLL' accuracy, since this 
provides the appropriate boundary conditions for the full 
NNLL result. Thus the result in Eq. (jT§5]) with a s (Q) 
is at an intermediate level of accuracy between NLL and 
NLL'. 

Note that the SCET expression Eq. ([185]) still shows 
the full dependence on the individual scales fiH,B.j,s in- 
stead of the single scale fj,, and the dependence on each of 
these scales cancels out to the order in resummed pertur- 
bation theory that we are working. The remaining scale 
dependence thus provides a useful way to estimate the 
theoretical uncertainty due to the truncation of higher 
order terms in resummed perturbation theory. 
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C. Scale Profile Functions 

In general there are three relevant regions with differ- 
ent power counting 

peak region: n ~ 2Aq C u/Qr < 1 , 

tail region: 2A qcd /Qr < n < 1 , (200) 

far-tail region: Ti ~ 1 . 

For the peak and tail regions of the distribution we have 
Ti < 1 and we must sum the large logarithms. In the 
tail region the results in Eqs. (|187[) and (|190[) above are 
the canonical scales for ^h,j,b,s for which the logs in the 
fixed-order hard, jet, beam, and soft functions are mini- 
mized. Evolution from these scales to another scale /z re- 
sums the logs of the ratios h/hh,j,b,s to all orders in a s . 
In the peak region for small T\, the scale /zs ~ Qt\ goes 
towards the nonperturbative region. The validity of our 
resummation analysis relies on there being a perturbative 
expansion for the soft function anomalous dimensions at 
the scale /is, rs[a s (Ms)] and 7s[ck s (m.s)]- Therefore in 
the SCET approach it is mandatory that we stop the 
renormalization group evolution at a scale /xs ~ 1 GeV 
that can still be considered perturbative. This requires 
the scales to deviate from the canonical form. Finally, for 
larger values of n the logs are no longer large, and the 
nonsingular terms in the fixed-order expansion become 
equally important. In this large t\ region we should turn 
off the resummation, which will revert the results to a 
fixed-order expansion in a s . Again this forces the scales 
to deviate from the canonical ones. 

To achieve these properties we use profile functions to 
describe the functional dependence of the scale ms.b,,/ on 
T\. First we will consider the profile functions for the r"' 
cross sections and then for rf. 



profile functions 



For the r^' cross sections, the canonical scales are 
given in Eq. (|187|) . /i s ~ n, Mb,j ~ ^/nQ, m ~ Q- The 
perturbative resummation of large logs of ratios of these 
scales is valid when Aqcd "C Ms -C Ms../ <C M-ff, which is 
the tail region. We will define boundaries , t± < t\ < 
for the region of n where this condition is satisfied, and 
use scales that are within a factor of 2 of the canonical 
ones. Beyond this region, when n > t 2 , T\ is of 0(1), and 
the logs are the same order as the nonsingular terms in 
the fixed-order expansion. In this region, the scales must 
be taken to be of the same order, /is — (J-b.j — Mff ~ Q, 
which turns off the resummation in Eq. (|185|) . Finally 
for T\ < ti, the soft scale approaches Aqcd and nonper- 
turbative corrections become important. In this region 
we freeze the soft scale /is used in the perturabtive cross 
section to a value above Aqcd: Ms ~ 1-2 GeV. The hard 
scale is \in ~ Q and beam and jet scale are determined 
by hard and soft scales as ^/JThJIs ~ (J-J.b- 



Profile functions for scales that satisfy the above crite- 
ria have been used for other cross sections in [l4|, [H, [ill ■ 
Here, we adopt the profile functions in [28| . The hard, 
beam, jet, and soft scales we use are given by 



(201) 



Ms,j(ti) 

Ms(n) 



l + eB,J0(i3-Tl)(l-^J vWrun(Tl , fJ,) , 



l + e S 6{h-T X ){l-j- 



Mrun(Tl , 



where eB,j.s are parameters used to vary the jet, beam, 
and soft scales to estimate theoretical uncertainty of the 
perturbative predictions. £3 is the point above which all 
scales are set equal, /j# = /j,b.j = Ms = t 1 - The common 
function jit r un(Ti, /i) is given by 



Mrun(n»M) 



/i + arf/ti 
2ar 1 +b 

M - a(n - t 3 ) 2 /(t 3 
Ji 



Mo - M 



U - t 2 



ta 



b = 



n<h, 
h < ti < ta , 
t 2 ) t 2 < n <t 3 , 
n>t 3 , 
Mti - Mo(*2 + ^3) 



ta - 1 3 



(202) 



The function Mrun(''"i,M) quadratically approaches ^ be- 
low t\ and fi above t%, and it is linearly increasing from 
t\ to t 2 . The continuity of MrunC?i)M) anc i its derivative 
at t\ and t 2 determines a and b. 

The default values of parameters we will use for what 
we will consider the "central values" of the r°' b cross 
sections arc: 



t 



M = Q , e B , j = e s = , Mo = 2 GeV , 

3 GeV „ . (203) 



Q 



i 2 = 0.4, *3 = 0.6. 



To estimate theoretical uncertainty due to missing higher 
order terms in fixed-order and resummed perturbation 
theory, we vary the parameters fi, es.j, and es from 
their default values by 0(1) factors in order to vary cor- 
responding scales fin, M-B.Ji and /is by 0(1) factors, re- 
spectively. We separately vary the parameters one by 
one and keep the others at their default values. The to- 
tal number of variations we perform around the central 
values are as follows: 



1) M = 2 ±1 Q ! e B ,j 

2) H = Q, e B ,j 

3) H = Q, e B ,j 



= 0, 
1 1 
±3^6 

0, e s = ± 1 -,± 1 - 



es = 
, e s = 



(204a) 
(204b) 

(204c) 



Variation 1 moves all the scales in Fig. 2] together up 
and down by factors of 2, and corresponds to the scale 
variation used to estimate the fixed order theoretical un- 
certainty in perturbation theory. Variations 2 and 3 are 
additional variations we are able to perform because of 
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FIG. 4: r^' b profile functions for the scales /in, /is,j(ri), 
fis(ri) with Q — 90 GeV used in the resummed factorized 
cross section Eq. (|185[1 . The double arrow and the colored 
bands illustrate the scale variations in Eq. (|204|l used to ob- 
tain theoretical uncertainty estimates. 

having independent fij^B and fis scales in the resummed 
cross section Eq. (|185[) and give an estimate of the uncer- 
tainty at each order in logarithmic accuracy in resummed 
perturbation theory that can not be achieved by varying 
the single scale ft. Variation 1 alone underestimates the 
total uncertainty. 

The size of the cross section at a given value of t\ 
may not vary monotonically with ej t B, es, and ideally we 
would vary them continuously within some finite band to 
find the maximum uncertainty. The four values we test 
for ej s, es in Eq. (|204p are a discrete approximation to 
such a procedure that remains computationally tractable. 
We take the largest and smallest values of the cross sec- 
tion among these points and use them to define the width 
of the uncertainty band from ej.B or es variation. 

To make a conservative estimate of the total uncer- 
tainty, we sum in quadrature the uncertainties we get 
from variations 1, 2, and 3 individually. We find that 
the total size of the bands provided by Eq. (|204[) are rea- 
sonable estimates of the theoretical uncertainty when we 
compare the cross sections at different orders of logarith- 
mic accuracy. 

Fig. [4]shows profile functions for Hb,j(t{), ms( t i) 
with Q = 90 GeV. The solid lines are the central values 
of the scales with default values in Eq. (|203|) , the double- 
headed arrow implies variation 1 and the bands represent 
variations 2 and 3 in Eq. (|204[) . The dashed and dotted 
lines are the canonical scales in Eq. (|187|) . 



2. t[ profile functions 
For rf, the canonical scales in Eq. (|190p are 

hh~Q, Ms ~ QV%(n - 1 + v) , 

1Mb ~ Qy/x{n -l + y), hj ~ Qy/n^T+y , (205) 




FIG. 5: Ti profile functions for the scales /j, c h , (4% j(tl), /is( r i) 
with x = 0.1, y — 0.9, and Q = 90 GeV, along with the simple 
canonical scales Eq. (|190l) . The double arrow and colored 
bands illustrate the scale variations in Eq. (|204p . 

where they satisfy the relation fi B j = x ±x / 2 fiHfts Com- 
pared to the canonical scales for r"' 6 in Eq. (|187|) . there 
are two differences in the canonical scales for rf. First, T\ 
is replaced by t± — (1 — y) because the transverse momen- 
tum of jet is nonzero, which is (I — y)Q 2 at tree level and 
the projection onto q c } differs from the projection onto 
jet axis by (1 — y)Q 2 - This requires that canonical scales 
in Eq. |(T57|) and profile in Eq. pTTj) are shifted by 1 - y. 
Second, the soft scale and beam scale are multiplied by 
yfx because of rescaling of the beam axis from xP for 
T^' b to P for rf . 

In this paper we consider the case yfx ~ 0(1) and this 
factor changes the scales by 0(1), which is the size of pcr- 
turbative uncertainties from varying ft, es,j and es- This 
means that multiplying [Ib and [is by \fx in Eq. (|201[) 
should not make a difference within the perturbative un- 
certainty. So, we could use the profile in Eq. (|201[) but 
shifted by 1 — y. On the other hand, by modifying the 
profile the canonical relations among the scales (Ih.j.b.s 
for rf Eq. (|190p can be maintained and we can account 
for the extra factors of s/x. Therefore, for these profiles 
we define ft c H j B s as 

= Mj Mb„/,s( t i) = Vb,j,s{x, n-l + y), (206) 

Vj(x,ti) = I + e,/(9(t3-Ti)(l-^-) V / MMmn( a; . T l,M,0), 
HB(x,Tl)= l + e B 0(t 3 -Tl)(l-j-^ y/flfJ,^ un (x,T!,fl, 1), 

Ms (x, n ) = [l + e s 6(t 3 - n ) (l - g-) *] fi c mn (x, n , ft, 1/2) . 

The Mmn used here depend on x and index 0, 1,2 that 
is different for ptj, /is, fts- We want ft° un (x , t\ , ft , n) ~ 
x n T\fi with n = 0, 1/2, 1 so that the canonical scaling for 
ft.j,s,B in Eq. (|190p is respected in the small t\ region. 
In the large t\ limit, fi^ nn (x , T\ , ft , n) should go to ft, so 
that fis and fts,.j both go to fi. 

As in Eq. (|202p ft^ un should run linearly between t\ 
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and t-x- However, the slope of fx^ un in Eq. (|206p should be 
different for the three cases n = 0, 1/2, 1. Therefore, we 
cannot use Eq. (|202j) to define n^ nn because all parame- 
ters in // run are fixed by matching boundary conditions 
and the slope is fixed. Instead, by replacing the quadratic 
polynomial in Eq. (|202|) by a cubic polynomial one can 
introduce a free parameter and this parameter can be 
chosen such that H^ un (x , t\ , /i , n) ~ x n T\[i between t\ 
and £2- We define ^f un as 



x n i// + a (n)r 1 2 /£i 
2a(n) n + 

Mcubic 



ri < £1 , 
£1 < n < £ 2 
* 2 < n < £ 3 
n > £3 , 



Mcubic(a; , n , fi ,n) = /i-c(n) 



-d(ra) 



n-<3 



,£3 — £2/ \*3 — £2 

fc(n) = x"- 1/2 /io - a(n)h , 

c(n) = 3(/i - x n -^ 2 fi ) ~ a(n)(2£ 3 + 4£ 2 - 3£ x ) , 

d(n) = 20 - x n - 1/2 fi ) - 2a(n)(£ 3 + £ 2 - £1) • (207) 

Here the parameters b(n), c(n), d{n) are determined by 
continuity of /i run and its derivative at £1, £2, £3. The 
slope a(n) is a free parameter which is chosen to satisfy 
a(n) ~ to achieve canonical scaling of jet, beam, and 
soft scales: 



a(n) 



/' 



-1/2 



Mo 



£3 + £2 - £1 



(208) 



Note that in x -> 1 limit, Eq. (|2U7|l reduces to Eq. (f202"|) 
and profiles for rf in Eq. (|206|) reduce to the profiles in 
Eq. (j20T|) for rf and rf . 

We choose the same default parameters and scale vari- 
ations as for rf ,h in Eqs. (|203[) and (|204p except for £2: 



£2 = 0.1 



(209) 



Because of the different definition of the profiles for rf 
this value of £2 must be smaller than the value for the rf' b 
profiles. This occurs because ^ run in Eq. (|207p changes 
faster than that the /i run in Eq. (|202p between £2 and £ 3 . 
As can be seen from Fig. [5] the final profiles for fi$ have 
similar shapes. 

Fig. [5] shows rf profile functions for /z 1 ^, \i c B j(ti), 
jU|(n) defined in Eq. ([206]) with x = 0.1, y = 0.9, and 
Q = 90 GeV. The solid lines are the central values of 
the scales with default values in Eq. (|209[) for £2 and in 
Eq. (|203[) for all other parameters. The double-headed 
arrow represents variation 1 and the uncertainty bands 
are variations 2 and 3 in Eq. (|204[) . The dashed, dotted, 
and dotted-dashed are the canonical scales in Eq. (|190j) . 



D. Nonperturbative Soft Function 



The hemisphere soft function defined in Eq. (|134[) de- 
scribes soft radiation between jets at the nonperturbative 



scale Aqcd as well as at perturbative scales above Aqcd ■ 
The results given in Eqs. (|160[) and (|D14[) are valid in the 
perturbative region. In the MS scheme the soft function 
valid at both scales is given by a convolution between a 
purely perturbative fu nctio n and a nonperturba- 
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tive model function F 110 



'hemi 



dk'S^.(k-k',^)F(k') 



The function F(k) contains information about physics at 
the nonperturbative scale and has support for k ~ Aqcd , 
falling off exponentially outside this region. Inserting 
Eq. (I210p into the factorization formula in Eq. (|140p one 
obtains the convolved form for the cross section: 

where da pert / dri is the cross section calculated by using 
only the perturbative soft function and Qr is given by 
Eq. ([54")) . Eq. (|21ip correctly describes both the peak 
region Qrt± ~ Aqcd where the entire function F{k) is 
required, as well as the tail region Qrt± 3> Aqcd where 
only its first moment is required since we can expand in 

AqcdAQhti). 

For the peak region, various wa ys to paramet rize mod- 
els for F(k) have beenproposed [zE EM [Tq3|. We will 
adopt one proposed in [76j that expands F systematically 
in an infinite set of basis functions: 



m - i 



N 



^ c «./n ( x 



n=0 
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where in principle we can choose any complete basis of 
functions /„. We adopt the same basis that has already 
been used in [l4|, , and exhibits fast convergence of the 
expansion. The normalization condition J dkF(k) = 1 
gives the constraint J2i c i — 1- The characteristic scale 
A of size O(Aqcd) is an additional parameter if the sum 
is truncated at finite N, as we will do in practice. 

In the tail region where QrJi ^> Aqcd, Eq. (|21ip is 
consistent with the power correction from an operator 
product expansion, 



dajri) 
dr\ 



a,b : c 




(213) 



To lowest order in Aqcd/(Qti) this result agrees with a 
simple shift t\ — » T\ — 2Q,i / Qr. Here the coefficient of the 
power correction 2f2"' b,c is a nonperturbative matrix ele- 
ment and it corresponds to the first moment of the non- 
perturbative function J dkk F(k) which could in princi- 
ple differ for each of rf ' b ' c . The first set of power correc- 
tions indicated on the second line of Eq. (|2 1 3[) comes from 
perturbative corrections to the leading power correction 
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[72j . and the second set involves purely nonperturbative 
corrections at subleading order. In the next section we 
will consider the question of universality of the f2J, f2j, 
f2J parameters for the observables rf , jf, r-f. 

In the peak region the parameters Ci and A should be 
determined by fitting to experimental data. Since data is 
not yet available, our only purpose here will be to get an 
idea of the impact of the nonperturbative shape function. 
We take the simplest function F(k) with N = 0. Then, 
c = 1 by normalization and A is the only parameter. To 
get the right first moment, we require A = 2^. We use 
f^i = 0.35 GeV, which is determined from measurements 
of e + e~ — > dijets (TJ. However, Oi in DIS is not nec- 
essarily the same as in e + e _ collisions, and we merely 
consider this to be an illustrative but reasonable value. 



E. Universality Classes for fii Parameters Denned 
with Different Directions 

The various versions of 1-jettiness r^' b ' c or the generic 
version Eq. (f2~4"|) depend on different choices of the axes 
qj = ujjtlj/2 and q B = ujbtib/2. In this section we will 
show that the 1-jettiness power correction parameter is 
universal under changes to the axes used in its defini- 
tion, by exploiting properties of operators [70l . [7l| and 
including hadron mass effects [§9|, ■ 

If we use different axes for the decomposition of four- 
momenta then they can all be written in a form similar 
to the event shapes given in [721 ]: 
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2Qr ^ 



IB) 



(214) 



where Qr is defined in Eq. (|54[) . i sums over hadrons, 
and Tn^~,ri,yj B are defined with respect to the vector 



qj,B by: 
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y, 
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iln^ 

2 qj-p 
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where m is the mass of the hadron whose momentum is 
p^. For the 1-jettinesses T\ given in Eq. (|2~4"]) we have 



f(r,y) 



(216) 



For each different n, i.e. each choice of qj.B, the defini- 
tion of m 1 - and y,j b change since they are computed with 
different coordinates. The Qr also depends on qj -qs, as 
given in Eq. (j54")l . 

Following the logic in [72[ for massive hadrons and 
[70l [7lJ for massless particles, the leading power correc- 
tion in the expansion Eq. (|213[) of distributions in event 
shapes of the form Eq. (|214|) is always described by the 
nonperturbative matrix element 



/>1 poo 

2n a { bc = / dr / dy JB f(r,y 
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x (0\YjY nj £^ c (r,y JB )YlY nB |0) 



Here £ t is a "transverse velocity operator" defined as in 
[zH , but now using the axes given by q j and q B , 

S T (r,yj B ) \X) = ]T miS(r - r t )6(y JB - y jB ) \X) . 

iex 

(218) 

It measures the total transverse mass of particles flowing 
in a slice of velocity and rapidity around r and J^jb- 
Now consider making an RPI-III transformation [§H in 
the matrix element in Eq. (|217[) which takes nj — > nj/£ 
and riB — > C n B- This transformation leaves the vacuum 
and the Wilson lines Y nj and Y nB invariant, but shifts 
£T{r,y JB ) to £ T (r,yjB+y') where y' = ln(. This is 
the analog of the boost argument for back-to-back nj 
and n B in Ref. [70L l7l| . Thus, the matrix element inside 
the integral in Eq. (|21T|) is independent of J^jb , and we 
can integrate over yjB to obtain the power correction 



Sij' ' for 



using the / given in Eq. (|216[) : 



dyj B f(r,y JB ) 



n J 1 B (r, f i) = 2nr(r, f i), (219) 



JB, 



where the renormalizcd matrix element is 

n( B (r, M ) - (0| Y^ B Y nj £ T (r, 0)1^1^ |0) . (220) 

This matrix element still depends on the choices of axes 
through n B ^j. By rescaling nj and n B as in Eq. (|133p 
we find it is independent of nj • n B . It still depends on 
these axes through the parameter r, since the transverse 
momenta p± inside r depends on the choice of these axes. 

However, in the tail region the fl{ B (r) always appears 
inside an integral. At LL order we have the resummed 
coefficient C{ jL (A:, r, /1) from 72| for any t\ and the shape 
function OPE is 



F(k) = S(k) + drC^(k,r,^)2n{ B (r,^+0[-^-) 



A? 
A? 



= S(k)+j drC 1 LL (fc,r, / i)2^ 1 (r, / .)+o(-p^), 

(221) 

where in the second line we removed the JB superscript 
on fii by using the fact that the only axis dependence oc- 
curs through the parameter r which is now just a dummy 
variable. It would be interesting to consider the univer- 
sality beyond LL order for this Wilson coefficient. 

Thus we see that at least to LL order there is 
a universal power correction Oi(r) for all three ver- 
sions of 1-jettiness, t°' ,c . Taking the tree level result 
C^ h (k,r,fj,) -> -5'(k) yields Eq. ([2T3"|) and leads to the 
identification 



dr fii(r) 



(222) 



Eq. (|22ip also implies universality of the shift parameter 
appearing in Eq. (|213[) for r"' '°: 



(223) 
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FIG. 6: Cumulant cross section in t° at Q = 80 GeV and 
x = 0.2. Colored bands show theoretical uncertainties around 
central values (lines) to LL (dotted line, green band), NLL 
(dashed line, blue band), and NNLL (solid line, red band) ac- 
curacy and the horizontal dashed line is the total cross section 
at fixed x,Q 2 . 




FIG. 7: Weighted differential cross section in rf at Q = 
80 GeV and x = 0.2. Colored bands show theoretical uncer- 
tainties around central values (lines) at fixed order a„ (dot- 
ted line, gray band) and resummed to NLL (dashed line, blue 
band) and NNLL (solid line, red band) accuracy. 



VIII. RESULTS 

In this section we present our numerical results for the 
three versions of DIS 1-jettincss: rf , t\, and rf. Wc 
plot the cross sections accurate for small T\ resummed 
from LL to NNLL accuracy, and also the singular terms 
at fixed order 0(a s ) (NLO) for comparison. (We esti- 
mate the size of the small missing non-singular terms by 
comparing to the known 0(a s ) cross section integrated 
over all n.) We start by describing the rf spectrum in 
detail, and then compare the features of the rf and rf 
cross sections relative to the results for rf . We choose 
s = (300 GeV) 2 as in the HI and ZEUS experiments. 
For the PDFs, wc use the MSTW2008 [H^ set at NLO 
and include five quark and antiquark flavors excluding 
top. To be consistent with the a s used in the NLO 
PDFs wc use the 2-loop beta function for running a s 
and a s (m z ) = 0.1202. 

Wc present results for the cumulant cross section a c (n ) 
defined in Eq. (|183l) and the dimensionless distribution 
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FIG. 8: Differential cross section in rf at Q — 80 GeV and 
x — 0.2 in the peak region, NNLL with nonperturbative shape 
function taken into account (NNLL PT+NP, dashed, orange), 
and without NP shape function at fixed-order a s (NLO PT, 
dotted, gray) and resummed (NNLL PT, solid, red). 
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Note that both the cumulant a c (n) and the differential 
distribution da / dn are differential in x and Q 2 . How- 
ever, for notational simplicity we made their x and Q 2 
dependences implicit in this section. 



A. rf cross section 

In this subsection, we present results for the cumulant 
cross section a c (n) and differential cross section da /dn 
for the "aligned" 1-jettiness n = T \- 

Fig. |5] shows the rf cumulant cross section, defined 
by Eq. (fT55j) . at Q = 80 GeV and x = 0.2. In or- 



der to illustrate perturbative convergence the results re- 
summed to LL, NLL, and NNLL accuracy are shown. 
The bands indicate perturbative uncertainties by vary- 
ing the scales hh,b,j,s given by "profile functions" as 
described in Sec. I VII C 1[ and there is excellent order-by- 
order convergence, and beautiful precision at NNLL or- 
der. The cumulant cross section increases monotonically 
from the small rf region and begins to saturate near for 
large rf where the integral defining this cumulant be- 
comes that for the total cross section. There is a small 
gap between the total cross section at 0(a s ) (dashed 
horizontal line) and our NNLL cumulant at large rf , re- 
flecting the small size of nonsingular terms not taken into 
account in this paper. Note however that these terms are 
important at the level of precision of our cumulant cross 
section, and hence they will be considered in the future. 



39 




x 



FIG. 9: x dependence of rf differential cross section at 
Q = 80 GeV and = 0.1. Colored bands show theoreti- 
cal uncertainties around central values at fixed-order a a (dot- 
ted, gray) and resummed to NLL (dashed, blue) and NNLL 
accuracy (solid, red). 

We can characterize the da/dr® cross section by 
three distinct physical regions: the peak region (jf ~ 
2A QC d/Q), the tail region (2A QC d/<3 < T\ < 1), and 
the far-tail region (rf ~ 0(1)). We will do this with 
four plots. We first show the purely perturbative cross 
section to study convergence and the impact of resumma- 
tion compared to fixed order results. Next we show the 
impact of nonperturbative effects, which in the tail re- 
gion produce a simple shift in the distribution, and have 
a significant impact on the shape of the spectrum in the 
peak region. We also illustrate the dependence of the 
cross section on x and Q 2 at fixed rf. 

Fig. [7J shows the weighted differential cross section 
T?da/dr? at Q = 80 GeV and x = 0.2. The results 
are weighted by rf for better visibility because the dif- 
ferential cross section falls very rapidly with rf. In the 
tail region, the overlap in resummed results shows a good 
perturbative convergence from NLL to NNLL. The large 
deviation between NLO and NNLL shows the large effect 
of resummation and the underestimated uncertainty of a 
pure fixed-order result. In the peak region, NLO result 
blows up as (In n)/n, while the NLL and NNLL results 
converge into a peak due to resummation of the large logs 
to all orders in a s . Again the uncertainty bands overlap 
fairly well. In the far-tail region for larger n , the resum- 
mation effect becomes small and the size of the deviation 
is reduced. Near the far-tail region (n ~ 0.3), the NNLL 
curve begins to depart from the NLL band. In this re- 
gion the nonlogarithmic a 2 term and nonsingular terms 
neglected in our NNLL result may begin to be significant. 

Fig. [8] shows the differential cross section da/dri at 
Q = 80 GeV and x = 0.2 in the peak region at fixed 
order and NNLL resummed accuracy. Note that it is not 
scaled by t\ as in Fig. [7J In this plot, the NNLL result 
convolved with a nonperturbative shape function (NNLL 
PT + NP) is shown in comparison with purely perturba- 
tive fixed-order NLO and resummed NNLL results (NLO 
PT and NNLL PT). As discussed in Scc. lVIIDI we use the 
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Q 

FIG. 10: Q dependence of rf differential cross section at x = 
0.2 and t° = 0.1, with theoretical uncertainties at fixed-order 
a a (dotted, gray) and resummed to NLL (dashed, blue) and 
NNLL accuracy (solid, red) 

simplest shape function with one basis function N = 
in Eq. (|212[) with a reasonable choice Oi = 0.35 GeV for 
the value of the first moment just to illustrate the impact 
of the nonperturbative effects. For practical analysis, a 
shape function with more basis functions should be used 
and the parameters Cj, A in the model function Eq. (|2 L2[) 
should be determined from experimental data. In the 
endpoint region, there is significant change from NLO 
and NNLL due to the resummation of large perturbative 
logs, and there is another large change from perturbative 
NNLL to the result convolved with the shape function 
due to nonperturbative effects. As we move into the tail 
region, the size of nonperturbative correction reduces to 
©(Aqcd/^iQ) an d the correction simplifies to the power 
correction in Eq. (|213p . 

Fig. IH1 shows the weighted differential cross section 
xdcrj (dx dQ 2 dr{) as a function of x at Q = 80 GeV and 
rf = 0.1. Note that the lower bound x > Q 2 / s is set by 
the relation xys = Q 2 in Eq. Q and the constraint y < 1. 
The x dependence comes from the quark and anti-quark 
beam functions and the decreasing curves with increas- 
ing x are characteristic patterns of PDFs contained in 
the beam function. With decreasing x, NLO and NNLL 
curves rise faster than NLL curve because they contain 
the gluon PDF, which rises faster than the quark PDF, 
and whereas the NLL result only contains the tree-level 
beam function which is just the quark PDF. 

Fig. [10] shows the Q dependence of the differential 
cross section at x = 0.2 and rf = 0.1. Overall, Q de- 
pendence is mild. In the naive parton model the cross 
section is insensitive to Q because of the approximate 
scaling law in the Bjorken limit where Q, s — > oo with 
x fixed. This scaling is broken by logarithms of Q in 
QCD. It is also broken by the Z boson mass with the 
factors 1/(1 + m 2 z /Q 2 ) in Eq. (|158|) . As shown in the 
plot, well below mz = 91.2 GeV the curves vary gently 
in Q and near and above mz they increase due to the 
factor Q 2 /(Q 2 +m|). 
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tions at Q = 80 GeV and x = 0.2 and at Q = 40 GeV and 
x = 0.02. The difference at NLL is zero for both parameter 
sets. 



B. rf cross section 

The rf cumulant cross section is different from rf by a 
single term at NLO in Eq. (|185[) . The term contains lnz 
where z is integrated over from a; to 1, and so the term 
becomes larger for smaller x. Fig. [TT] shows their per- 
cent difference at NLL and NNLL for two sets of (Q, x): 
(80,0.2) and (40,0.02). The difference at NLL is zero 
because at LO fixed order rf and rf cross section are 
identical and the NLL logs are the same. At NNLL for 
x = 0.2 the size of difference is small, a few percent. The 
difference at the value x = 0.02 is larger than that for 
x = 0.2, becoming now a 10-15% effect. This difference 
is roughly constant in Q because of the mild Q depen- 
dence in Fig. [TOl 



C. rf cross section 

The 1-jettiness rf is designed to measure a jet close to 
the z axis (incoming electron direction), and the factor- 
ization theorem for rf in Eq. (|153[) is valid for a jet with 
small transverse momentum q\ = (1 — y)Q . So, the pa- 
rameters Q and x should be chosen such that 1-yCl 
in other words, Q 2 /(xs) ~ 1. The parameters in Fig. [5] 
cannot be used because y « 0.36 for Q = 80 GeV and 
x = 0.2. For rf in Figs. [12] and [13] we choose Q = 90 GeV 
and x = 0.1 for which y = 0.9. Note that the profile 
functions for rf given in Eq. (|206[) are also different from 
those for rf' 6 . 

Fig. [12] shows the cumulant rf cross section rcsummcd 
to LL, NLL, and NNLL accuracy. The most notable fea- 
ture in the rf spectrum is the threshold #(rf — 1 + y) in- 
dicated by an arrow in the plot. The threshold is exactly 
respected in LL and NLL results and is effectively true at 
NNLL because, although Eq. (|189b[) contains terms vio- 
lating this threshold at 0(a s ), their size is numerically 
small ( ~ 0.1%). In the region near this threshold nonper- 
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12: rf cumulant cross section at Q = 90 GeV and 
x — 0.1, giving y = 0.9. Colored bands show theoretical 
uncertainties around central values for resummed results to 
LL (dotted, green), NLL (dashed, blue), and NNLL (solid, 
red) accuracy. The horizontal line is the total cross section. 
The arrow at 1 — y indicates the threshold in rf spectrum. 




FIG. 13: Difference between rf cumulant cross sections in 
comparison to rf results at Q = 90 GeV and x = 0.1 which 
gives y = 0.9. The horizontal dashed line is the total cross 
section at this x,Q 2 . 



turbative corrections are quite important, and the purely 
perturbative cross section actually has a small negative 
dip (almost invisible in the plot). 

Fig. 1131 shows rf in comparison with the rf cumulant 
cross section at NNLL. In addition to the threshold dis- 
cussed in Fig. 1121 the rf curve increases more slowly than 
the rf curve does. This is because the normalization of 
the rf axes in Eq. (|42]) are different from those for rf . 
The beam axis qB for rf is larger than for rf by a factor 
of 1 jx while the jet axis qj is approximately the same in 
the limit y — > 1. This increases the projection of the par- 
ticle momentum qs ■ Pi by the factor of 1/x in 1-jettiness 
Eq. (|2"4"]) , but rf is not increased by quite the same factor 
because fewer particles are grouped into the T-Lb region 
due to the minimum in Eq. (|24|) . Still, in Fig. [13] for the 
same value of the cross section the departure of rf from 
its threshold is larger than that of rf due to this factor. 
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IX. CONCLUSIONS 

We have predicted 1-jettiness (ti) cross sections in DIS 
to NNLL accuracy in resummed perturbation theory, ac- 
curate for small t± where hadrons produced in the fi- 
nal state are collimated into two jets, including one from 
ISR. We used three different versions of 1-jettiness, r 1 a ' i,,c , 
which group final-state hadrons into "beam" and "jet" 
regions differently and also have different sensitivity to 
the transverse momentum of ISR relative to the proton 
direction. 

Each t\ is similar to thrust, measuring how closely 
final-state hadrons are collimated along "beam" and 
"jet" reference axes, but with important variations, t" 
measures the small light-cone momentum along two axes 
aligned with the proton direction and the actual jet di- 
rection, and averages over the transverse momentum of 
ISR in the calculation of the cross section. t\ projects 
onto fixed axes such that the beam and jet regions are 
back-to-back hemispheres in the Breit frame. The fixed 
jet axis is not quite equal to the physical jet axis in the 
final state, causing t\ to be sensitive to the transverse 
momentum p± of ISR and requiring a convolution over 
p± in the jet and beam functions in the t\ factoriza- 
tion theorem. Finally rf groups final-state hadrons into 
back-to-back hemispheres in the CM frame, projecting 
momenta onto the initial proton and electron directions, 
and also requires a convolution over the transverse mo- 
menta of the ISR and final-state jets. Furthermore, the 
case of small rf also requires the DIS variable y to be 
near 1. 

We proved factorization theorems for all three versions 
of Ti using the tools of SCET, carefully accounting for 
the differing dependences on the transverse momentum 
of ISR. These differences lead to the appearance of the 
ordinary beam function in the r™ factorization theorem 
and the generalized dependent beam function in the 
t\ and rf factorization theorems. We were able to re- 
late the soft function appearing in any of these factor- 
ization theorems in any reference frame to the ordinary 
DIS hemisphere soft function by suitable rescaling of the 
arguments, using boost invariance. 

The hard, jet, beam, and soft functions appearing in 
all the factorization theorems and their anomalous di- 
mensions are known to sufficiently high order to allow us 
to immediately achieve NNLL resummed accuracy in our 
predictions for the T^' h ' c cross sections. We gave predic- 
tions for the differential and cumulant t\ cross sections, 
illustrating the differences among r 1 l ' f> ' c due to the dif- 
ferent dependences on the transverse momentum of ISR. 
We presented numerical predictions at x and Q 2 values 
explored at the HERA collider, but our analytical pre- 
dictions can easily be applied to a much wider range of 
kinematics relevant at other experiments, such as at JLab 
and the future EIC @ and LHeC |z3- 

The resummed predictions we presented are accurate 
for small values of T\ where final-state hadrons are well 
collimated into two jets. For large n our predictions 



have to be matched onto fixed-order predictions of non- 
singular terms in t\ from full QCD. We leave the per- 
formance of this matching at 0(a s ) and beyond to fu- 
ture work. However, we compared our cumulant t\ cross 
sections for large T\ to the known total cross section at 
fixed x and Q 2 , and found that the cumulative effect of 
the matching corrections on the whole cross section is 
roughly at the several percent level for the kinematics we 
considered. 

We showed how to incorporate nonperturbative 
hadronization corrections into our predictions by inclu- 
sion of a shape function that is convolved together with 
the perturbative soft function. The first moment of the 
shape function gives the parameter fli which describes 
the shift to the distribution in the tail region. We demon- 
strated that this parameter is universal for our three 
event shapes r^ b ' c and for any values of x,Q 2 , and so 
it can be extracted from one set of data to predict oth- 
ers. We also made a simple illustration of the effects of 
a shape function numerically on the cross section. Wc 
leave a more extensive study of nonperturbative effects 
and extractions of the model parameters from data to 
future work. 

The extension of our results to 7V-jettiness tn in DIS 
with TV > 1 is straightforward, at least if we define tn 
similarly to the 1-jettincss rf that wc defined in Eq. 
That is, 



T N = 7X1 1^ mill {lB - Pull- Pi 



,q a N - Pi }, (225) 



where q B = xP and qf is the jet axis of the zth non-ISR 
jet in the final state as given by a jet algorithm or by 
minimization of the sum Eq. (|225|) over the directions of 
q", , . . , q^f. As long as these jet reference axes are aligned 
with the physical jet axes, the transverse momentum k± 
of ISR will not affect the value of at leading order in A. 
The factorization theorem will then look like Eq. (|148l) . 
with suitable generalizations of the hard and soft func- 
tions and additional jet functions (cf. (27j): 



da 
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dx dQ 2 dr^ dx dQ 



2 J dt 1 j---dt 1 ^ dt B dk s (226) 
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t N j k s 



1 N 



Q 2 Q 
x 'Y^B i {t B ,x,ix)J Kl {tj,n) ■ ■ ■ J KN (tj,n) 

x tr H i ^. K ({q m },L,(j,)S i 1 ^ K (ks, {q m }, fi) , 

where i?i_ >K ({g m }, L, //) contains the underlying hard in- 
teraction i(qB)e(k) — > e(k')Ki(qi) ■ ■ ■ Ki\r(qjv), where i, Kj 
denote parton types, L denotes the dependence on the 
leptonic states e(k),e(k') and the exchanged virtual bo- 
son, and the sum over i, k is over all relevant partonic 
channels. The hard and soft functions H, S are matrices 
in color space, and the trace is over these colors. £>; is 
the ordinary beam function for the initial-state parton 
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of flavor i. Since Eq. (|225[) uses reference axes <?" that 

arc aligned with the physical jet axes, the arguments tj 
of the jet functions are the invariant masses of the jets 
and are not shifted by any transverse momentum k± of 
ISR. Thus only the ordinary beam function Bi appears 
in Eq. ()226[) . k± having been averaged over. We leave 
the explicit evaluation of Eq. ()226[) for iV-jettiness cross 
sections in DIS with N > 1 to future work. 

Our results bring to the arena of DIS the power of 
SCET that has already vastly improved the precision of 
theoretical predictions of event shapes in e + e~ collisions 
and pp collisions. The factorization theorems derived 
here point the way to methods to improve the precision 
of parton distributions, hadron structure, and the strong 
coupling a s that we can extract from existing and future 
experiments. With further advances in our calculations 
to greater perturbative accuracy and improved modeling 
of the nonperturbative effects, the frontiers of the study 
of the strong interaction using jets in DIS can be pushed 
to higher precision. 

N ote Added: While this paper was being finalized, 
Ref. appeared which also considers the event shape 
we call rf at NNLL order. A complete derivation of the 
factorization theorem was not presented there, where the 
focus is instead the use of 1-jettiness to probe nuclear 
PDFs and power corrections from dynamical effects in 
the nuclear medium. 
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Appendix A: Generalized Rapidity Gap AY 

The f-jettiness t\ in Eq. (|48p is just one possible com- 
bination of jet and beam momenta that we can choose 
to measure in DIS. It is quite straightforward to keep 
nj -pj,nB -Pb as independent observables in the factor- 
ization theorem Eq. (|I35|) , and then to form other observ- 
ables by taking different combinations of nj ■ pj, ub ■ Pb- 
In this Appendix we consider one of these possibilities — 
the generalized rapidity gap AY between the beam jet 
and the other final-state jet. 

The rapidity of a particle with momentum p with re- 
spect to the z-axis is given by 



v n 1 1 Hz ' p 

Y n z n AP) = o m ' 

2 n z -p 



(AI) 



If p is n z -collinear, the rapidity Y Uz Ra is large and pos- 
itive, while it is large and negative if p is n 2 -collincar. 



Two jets produced in DIS are not, in general, back-to 
back, and the reference vectors that measure jets are not 
always aligned along one (z) axis, as Fig. |3] illustrates. 
The rapidity in Eq. (|AI[) can be generalized by replacing 
n z % with ris and nj as follows: 



Y njnB (p) = -]n—, 



(A2) 



where Y njnB is large and positive for the nj-collinear jet 
and is large and negative for the ns-collinear jet. The 
generalized rapidity difference between two jets of mo- 
menta pj and pb is given by 



AY = Y 



njn B (Pj) ~ YnjtiB (PB) = ^ ln 



1, n B -p,jnj-p B 



nj-pj riB-PB 

(A3) 

The tib.j in Eq. (| A3|) can be replaced by qs,.j because the 
energy factors u>j >B /2 in the numerator and denominator 
cancel. By using Eq. ([2Tjl qs.j-PB.j can be expressed in 
terms of tb.j- So, Eq. (|A3j) can be rewritten as 



(A4) 



AY = 1 1 J1J-PB1B-P J 

2 tjTbQ 4 



where the products 2qj-ps and 2qs-pj are 0(Q 2 ) and 
AF is O [ln(l/ v ^i)] ~ O [ln(I/A 2 )]. Eq. (S3J can 
be specified for DIS by using qj-ps ~ Qj'{P + Q) and 
QB'Pj ~ qB • q where we use momentum conservation 
P + q = Pb + P. j and suppress p 2 B and p 2 j. As we have 
three versions of Ti, there are three versions of AY: 



AY a - b = - In 



l-x 

XTjT B 



AY C = I ln -4 



X^TjTB 



(A5) 



Appendix B: Tensors and contractions 

The symmetric and asymmetric tensors transverse to 
both tib and nj are defined by 



n J n B + n J n B 

nj-n B 



(Bla) 



nu _ n J^J + ^>J = au _ n B n B + n B n B 

g 2 2 

nj-n B 

1 a-P _ 1 -a P 

2 e ^ 1 '<xP n J n J — ^ e ^ap n B n B J 



(Bib) 



where ris and nj are conjugate to ub and nj as defined 
in Eq. JSTJ . 

In order to calculate the contraction of the lepton ten- 
sor with the hard function H^^ v as in Eq. (|136l) . 
we must compute two tensor contractions: gT l/ g±' and 
e^e^, where g T ,e T are defined in Eq. (|B"5j) and g±,e± 
in Eq. (|Bip . These contractions are given by 

k^k' v + k v k'^ / ^ n^ + n^j 



4 



nj-riB 1 



nj-rtB 

(nj-kriB-k' + nj-k'riB-k). (B2) 
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and 



T iiv _ _ 



n,j-n B Q' 
4 

Uj-Ub Q 2 



e a e^ v k a k'^ n 1 n & 



(nj -k' hb -k — nj -knB -k') . (B3) 



The ratio Eq. (|B2) over Eq. (|B3} is the coefficient 
r(qj,q B ) defined in Eq. ([155]) . 

Appendix C: Plus distribution 

The standard plus distribution for some function (/(a;) 
is given by 

= lim [6»(a; - e) q(x) + S(x - e) Q(x)] , (CI) 

e-vO 

where 

Q(ar) = ^ dx'q(x') . (C2) 
Integrating against a test function /(a;), we have 
dx[6(x)q(x)] + f(x), 

"* dx [f(x) - /(0)] + /(0) Q(x max ) , (C3) 



for x max > 0. 

For the special cases q(x) = l/x 1_a with a > — 1 and 
<7(x) = In" x/x with integer n > 0, we define: 



£» = 
£„(x) = 



9{x) 



J + 
(x)ln n x 



(C4) 

n>0. (C5) 



For convenience we also define 

C-i{x) = 8{x). (C6) 
The plus function C n obeys the rescaling relation, 



k=0 



\£ n (\x)=^( n k \ln k \£ n _ k {x) 



ln" +1 A 
n+ 1 



<*(*), (C7) 



where A > 0. 



Appendix D: Renormalization Group Evolution 

In this appendix we collect results relevant for resum- 
mation of the DIS 1-jettiness cross section Eq. (|140p and 
its special cases Eqs. (|T37)l . (|T52"]) . and (fT35|> for rf' 6 ^. 



The RGE and anomalous dimension for the hard Wil- 
son coefficient C in Eq. (|154[) for the two-quark operator 
are EH!! 



7 £(<z 2 , M) = r? usp (a s ) In =|- + 7* (a s ) . (Dl) 

A 4 

The anomalous dimension for the hard function in 
Eq. (|155|) is given by 

1H {Q 2 ^) = 2r« usp (a s ) In ^ + 7B (a s ) , (D2) 

where 7# = 27^,. The expansions in a s of r^ usp (a s ) and 
7^(a s ) are given below in Eqs. (|D28j) and (|D29j) . 

The solution of the RGE in Eq. CEQ) yields for the RG 
evolved hard function: 



U H (Q 2 , fM) , f i) = e K »^(Q-y 



H(Q 2 , ft) = H(Q 2 , fi ) U H (Q 2 ,fio, M) ■ 

K H (fj, ,fx) = -&K vq (no,n) + K 1H (fj, ,fj,) , 
J7j?(Mo,m) =477t«(mo,m) 5 ( D3 ) 

where the functions Kpi{no, A 4 ), Vr^ino, 1^) and -K" 7 are 
given below in Eqs. (ID24I) and (|D26|) . 

The quark beam function RGE is given by 

^B q (t, x,n) = J dt' 7 £(f - A 7 , /i) B 3 (i', x, /i) , (D4) 

lB (t, fj) = -2T« usp (a s ) ^£o(^) + 7M«.) *(*) . 

and its solution is 0, [51 [HI [mj 



Bq(t,x,n) = J dt' B q (t - t',X,flo) UB q (t', /Uo, A 4 ) 

p K Bq ~~fE VB q / A \ 

r7 B9 (t, M0)At ) = — — — - - +<s(t 

r(i + r? B j L y^o v av 

As 9 (a 4 o,A 4 ) = (a 4 o , A 4 ) + # 7 « (a*o , A 4 ) > 
VB t (fio,fj) = -2r}r<i(jio,fi). (D5) 



The solution of the RGE for S g given by Eq. (|D5[) can 
be derived from the form of the solution Eq. (|D3|) for 
the hard function by first Laplace transforming the beam 
function: 



B q (v,x,n)= dte- vt B q (t,x,fi), (D6) 



which obeys the RGE 



V-^B q (v,x,n) = ^ Bq {v,^)B q (iy,x,n) , (D7) 
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with the Laplace transformed anomalous dimension, 
%{u,n) = 21* (a s )ln( M 2 ^) + 7 |K). (D8) 



The evolution of B q in Eq. (|D7[) is multiplicative, of the 
same form as the hard function RGE Eq. (|D2[) . and there- 
fore its solution is just like the hard function Eq. ()D3[) . 
given by 



B q (v,x,n) = B q {v, x, no)U Bq (v, hq,h) . 



(D9) 



where 

U Bq (v, Mo, n) = e KB ^'i*)(jj^ve"' E )- r i B '^ '^ , (D10) 

with K Bq , TjB q given by the same expressions as in 
Eq. (|D5|) . The inverse Laplace transform of the solu- 
tion Eq. (|D9[) gives the momentum space solution for 
B q (t,x,fj.) in Eq. (|D5j) . 

The jet function obeys the same RGE as the beam 
function. They are defined by matrix elements of the 
same operator. The solution for the Laplace transformed 
jet function J q (v, fi) is given by the same form, Eqs. (|D9[) 
and (|D10[) with B — > J, and for the momentum-space jet 
function J q (t, fi) by the same form Eq. (|D5|) . with B —> J. 

The hemisphere soft function in Eq. () 160|) obeys the 
RGE 



d f 

H-j^S hcmi (kj,k B ,n) = J dk'jdk' B 



(Dll) 

x ls(kj -k'j,k B - k' B ,n)S hcmi (kj, k' B ,fi), 

where the dependence o f the anomalous dimension on the 
two variables separates jl08l | : 

1s(kj, k B ,fx) = -fs(kj,n)S(k B ) + js(kB,n)5(kj) , 

. (D12) 

with each piece of the anomalous dimension taking the 
form 

7s(A,A») = 2L? usp (a s )i£ (^) +js(a s )S(k) 7 (D13) 

where 75 = —7^ — 7^. The solution to the soft RGE 
Eq. (|D11[) is given by 



Shemi(fcj, k B ,fl) = / dk'jdk' B Shemi(kj, k' B , fJ,o) 



x U s (kj - kj,n ,n)Us(k B - k' Bl fi ,n) 



(D14) 



where 



Us(k,(M),(i) 



T{l + vs) 



'-£ r > 



S(k) 



Vs(fM),fi)=2rir 9 (fio,fi). (D15) 

This solution can be derived as for the beam and jet 
functions above by first taking the Laplace transform: 



I dkj / dk B e 
'0 Jo 



- Vjkj - VBkB S hemi (kj,k B ,lJL), 

(D16) 



(D19) 



which obeys the RGE 

H^-Sh e mi(vj,V B ,IJ>) = Shend{vj,VB,fi) (D17) 

d\i 

x \lsiyj,[L) +7s(^s,Ai)] , 

where each part of the anomalous dimension takes the 
form 

7s(f,A«) = -Sr^pM/x^) + ls (a s ) . (D18) 
Solving the soft RGE Eq. (|DT7|) . we obtain 

Shemii^Jj^B, A 4 ) = Shemii^J, V B , Ho) 

x U s {vj,Ho, v)Us{vB,IMi,n) 
where each soft evolution factor takes the form 

U s {v,»o,») = e Ks ^\we^)-^^ , (D20) 

where Ks^Vs & re given by Eq. (|D15|) . Taking the in- 
verse Laplace transform of Eq. (|D19[) gives the solution 
to the RGE for the soft function in momentum space 
Shemi(kj,k B ,(J.) given in Eqs. (|D14|) and (|D15|) . 

In the 1-jettiness cross sections in this paper, we always 
encounter the soft function Eq. (|D14[) projected onto a 
function of a single variable k, according to Eq. (|134[) . It 
obeys the RGE 

d f 
IJ>-j^Sh em i(k, n) = J dk! 2j s {k - k', fi)S hc mi(k' , (J>), 

(D21) 

where 73 (fc,//) is given by Eq. (|D13|) . In Laplace space, 

M-T"^hcmi(^, A 4 ) = 275(1/, /i)5hemi(l / , A 4 ) • (D22) 

dfi 

The solutions to these RGEs are given by 
Shcmi(&, A*) = Jdk'Si mril (k',Li )Ug(k-k',no,tJ-) (D23a) 

Shemi^ , A 1 ) = <Shenu(f, Vo)U S (v, /i , /i) 2 , (D23b) 

where Ug(k, (j,q, (J,) is given by Eq. (|D15[) with Ks,r]s — >• 
2K s ,2rj S : and U s (v,Ho,n) is given by Eq. (|D20j) . 

The functions i^r? (a*o, A*); f/Ti (mo. A 4 ), -^7^0, A 4 ) m the 
above RGE solutions are defined as 



/a. 0*o) ^( as ) 

771^0, A 4 ) = / 7J7-^ r cus P (a s )> 
K^{no,n) = / _l- 7 ( aa ). 



da' s 



(D24) 



Expanding the beta function and anomalous dimen- 
sions in powers of a s , 



A47T 



71+ 1 



(D25) 



n=0 

n+l / \ n+1 



V47T/ 
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their explicit expressions to NNLL accuracy arc (sup- pressing the superscript q on T q ), 



?7rOo,Ai) 



T J 4tt 
4/3 2 I a, (mo) 

Q s (^o) 
47T 



rp_ 

'ifa 

7o 



1 In r 

r 



ft 
1% 



fa 



£1 

r 



-i ) (1 - r + lnr) + In 2 r 



In r 



lnr - 


a s (^o) 






4-7T 






lnr - 


a s (Mo) 


( 1± 


_M 


4-7T 


V7o 


a; 



(r-l) + 



(r-1) 



/? r 

16tt 2 lr 



«2\ 

-^J(l-r + rlnr)- 
giTi | /? 2 fe\r 2 

ATo A 2 A> , 



£2 
r 



faTi\(l-r) 

A>r 



(D26) 



Here, r = a s (fj,) / a s (fj^) and the running coupling is given to three-loop order by the expression 



A' 



A 



In A 



16tt 2 



fa 



1 - 



A 



Af / hi A 



/? 2 V A 



A 



- 1 



(D27) 



where A = 1 + a s (/xo)A) m (/V/- t o)/(27r). In our numerical analysis we use the full NNLL expressions for ifr. 7 , )jr in 
Eq. (|D26|) . but to be consistent with the value of a s (ii) used in the NLO PDFs we only use the two- loop truncation 
of Eq. (|D27|) . dropping the fa and 01 terms, to obtain numerical values for a s (fj,). (The numerical difference between 
using the two-loop and three-loop a s is nu merically very small and well within our th eory unce rtainties.) Up to three 
loops, the coefficients of the beta function jll4 Ill5| and cusp anomalous dimension [l!6l 1 1 1 7f | in MS are 



11 n 4 r 
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The MS non-cusp anomalous dimension jh = 27^. for the hard function H can be obtained 38l 11£ 
divergences of the on-shell massless quark form factor C (q 2 , /i) which are known to three loops [1 1£ 



(D28) 
from the IR 
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As shown in 1071 ). the anomalous dimension for the beam function equals that of the jet function, 7^ = 7j, so the 
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non-cusp three-loop anomalous dimension for the jet and beam functions are both given by [3S 
6C F , 

' 1 AR 2 



7b o 



7% 



Til 



7b 2 
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7j2 



2C 
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/ 52019 841tt 2 



4tt / + 48C 3 )C , f 
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The anomalous dimension for the soft function is ob- 
tained from 75 = —7^7 — 7 |j- At NNLL, we only need 
the one- and two-loop coefficients of Jh.b,j.s- The three- 
loop coefficients are given for completeness. They would 
be required at N 3 LL, along with the four-loop beta func- 
tion and cusp anomalous dimension, the latter of which 
has not yet been calculated. In addition, the full N 3 LL 
result would also require the two-loop fixed-order correc- 
tions, which are known for the hard function, but not yet 
for the beam and soft functions. 



Appendix E: Coefficients in Momentum-Space 
Resummed Cross Section 



the coefficient I qq,qg (t, z, /i) inside the beam function in 
Eq. (|171j) . The index np = 1 for the soft function and 
iip — 2 for the jet and beam function. In the case of the 
beam function, the z dependence in F(t,(i) is implicit. 
The coefficients F n in Eq. (|E1[) for the three functions 
are S n , J„, and I qq,qg . The soft coefficients at order a s 
are given by 



5_i (a.) = 1 + 



a s CF 7T 2 

Itt 3~ 



S (a s )=0, Si(a.) = ^£(-16), 



(E2) 



in 



the jet coefficients by 



The resummed cross sections for rf' ' c in Sec. IVUl arc 
obtained by plugging the solutions to the RG equations 
for the hard function and for the momentum-space jet, 
beam, and soft functions given in App.[TJ]into the factor- 
ization theorems derived in Sec. lVDl Performing the con- 
volutions in these factorization theorems of the jet, beam, 
and soft evolution kernels given in App. |TJ] and fixed- 
order functions requires computing the convolutions of 
plus functions with each other. The results of these con- 
volutions produce the expressions given in Eqs. (|185|) and 
(|188[) . given in terms of coefficients J n , I n , S n of the logs 
in the fixed-order jet, beam, and soft functions and coef- 
ficients V™ n and V£ (a) that are the result of the convo- 
lutions of plus functions. In this Appendix we tabulate 
these coefficients. For more details see Refs. [13. [T6(] . 



1. Jet, Beam, and Soft Coefficients J n , I? l q ' q9 , S n 

The fixed-order results at 0(a s ) of soft, jet, and beam 
functions can be written as sum of plus distributions as 

G(t,n) - 4^ E G n [a.(n)]Cn (-M . (El) 



n=-l 



/'" 



where G(t, fi) represents the single- variable soft function 
S(t,fi) in Eq. (|2TU| . jet function J(t,fi) in Eq. ([TBI?]) , or 



J-i{a s ) = 1 + 

o. s Cf 3 

Jo(a s ) = 



ck s Cf (7 it 



, Ji(a s ) = 



a s C F 



(E3) 



7T 4 

and the beam function coefficients by 

a s C F 



Il\(a s ,z)=£- 1 (l-z) 



2n 



d(l - z)(l + z 2 ) 



Il q {a s ,z) 



ir(a s ,z) 
and 

I$ g (a s ,z) 



-^-£- 1 (l-z) + e{l-z)(l - z - liflnz 
6 V 1 — z 

^£e(z) (p qq {z)-\c^{i-z) 



2tt 



2£_x(l-z). 



(E4) 



?^e(z)[p qg (z)ln^+26(l-z)z(l-z) 



2tt 



e(z)p qg (z) . 



(E5) 



where the splitting functions P qq qg (z) are given in 
Eq. (WD. 

The argument of the plus distributions L n in Eq. (IE1|) 
can rescalcd by using the identity Eq. (|C7[) . Eq. (jEll) can 
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be rewritten in terms of the rescalcd distribution as 

1 M-lj 



n— — 1 v 



G ( t it x ) = -r~^ Y G n [a s (ji),\]C r 



(E6) 



where the coefficents G n (a s , A) in Eq. (|E6[) are expressed 
in terms of the coefficients in Eq. (|E1[) by using the rescal- 
ing identity in Eq. (|C7[) as 



G 



ln n+1 A 



-i(a a , A) = G_i(a s ) + V" G n (a s ) — — , 

z — ' n+l 



G n (a s , A) = 2J 



k=0 



(n + k)\ 
i\k\ 



G n+k (a s )\n k X, (E7) 



where G„ = {S n , J n , In 9 ' 99 }- Explicit expressions for 
S„ (a S7 X), J„ (a s ,X), and I qq ' q9 (a s , A) are obtained by 
inserting Eqs. ifE^ . (|E3| . pi]) , and (jESJ) into Eq. (jE7| '. 



2. Results of convolving plus functions 

Convolutions of plus distributions in the jet, beam, 
and soft evolution kernels and the fixed-order functions 
produce the functions VJ?(f2) and the coefficients V k mn in 
the resummed cross sections Eqs. (|185|) and (|188|) . There 
are three types of convolutions of plus distribtions C n and 
C a . and wc write them in useful form as 



m+n+l 



( dyC m (x-y)C n (y) = £ V e mn £i(x) 

J n i 



dy [aC a (x -y) + S(x - y)] [bC\y) + 5(y)] 



(E8) 



r(i + o)r(i + 6) 



r(l + a + 6) 

P ii'~r -l 

dy[aC a (x-y) + 5(x-y)]C n (y) = ]T V?{a)C%(x). 



(a + b) [C a+b (x)+6(x)] 

n+l 



fc=-l 



The coefficients V k n (a) and V k mn are related to the Taylor 
series expansion of V{a, b) around a = and a = b = 0, 
where V(a, b) is defined by 



V{a,b) 



T(a)T(b) 1 1 



T(a + b) a b' 
which satisfies V(0, 0) = 0. The V fc ™ (a) for n > are 
d" V(a,b) 



(E9) 



db n a + b 



b=0 



V k n (a) = 



n\ d 



n—k 



k Ab 



^V(a,b) 



b=0 



+ S kn , < k < n , 

k = n+l. 
(E10) 



The V k mn are symmetric in m and n, and for m, n > 
they are 



/■ d m d" V(a,6) 



da m d&" a + b 

rn n 



a=b=0 

m\ ( n 



p=0 q=0 ^ ' vy 

x-f ^ y(a,6) 



. m + 1 n + l 



=6=0 

< k < to + n . 



fc = m + n + 1 . 



(Ell) 

Using Eq. (|C6j) we can extend these definitions to include 
the cases n = —1 or m = — 1. The relevant coefficients 
are 



y_- 1 1 (a) = l ! V r " 1 (o) = a 
K fc S 1 1 (a)=0, ^ = ^ _1 = ^ 



(E12) 



Appendix F: Resummed cross section from Laplace 
transforms 



An alternative way [3a . l8G | to express the resummed 
cross sections in Sec. IVIII is to utilize the Laplace- 
transformed jet, beam, and soft functions given in App.lDl 
and their RGE solutions. The method avoids taking ex- 
plicit convolutions of plus functions in the evolution fac- 
tors and in the fixed-order jet, beam, and soft functions. 

Each of the RGE solutions for the jet, beam, and soft 
functions is given by a function of the form 

G{y,n) = G(^,^o)e KG(AIO ' Al) [Mo(^e 7E ) 1/jG ]" JG ' ?G(AI0 ^ ) ■ 

(Fl) 

For the jet and beam functions, j G = 2, while for the soft 
function j G = 1. The fixed-order expansion of G(v, fj^) = 
G(L g ,/iq) can be considered to be a function of the log 
L G = lnQc/Mo, where Q G = (i/e?*)-Vi«. To 0(a 2 s ), 



G(L G ,^ ) 



/ q s (^o) y 
V 4tt / 



■in 



-T G L G - lG L G + c G ) (F2) 



(T G yL G + r G h G + -f3 Q \L 



+ {\{lG?+l°G^-T G -c} G T G )L 
- (7g + c G 1g + 2c G ML G + cl 



Each power of L G can be generated by taking derivatives 
with respect to rj G in Eq. (|FL[) : 

G(y,ti = e ^^^)5(a, F , / io)[Mo(^ 7B ) 1/jG ]" JG " G(AI °'' i) , 

(F3) 
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where g(d n ,(io) is the operator constructed by replacing Now it is easy to take the inverse Laplace transform of 
each L G in Eq. (|F2|) with d v /j G : 



_ , , ( r o d n o d v i 

- 1 H I _i G~2 ~ 7g~ r C G 



47T V 



(F4) 



+ ( ^ r ^ (r „ )2 | + r „ (7S + |A) | 



4 



<9 4 



+ U(^) 2 + 7 G /?o-r G -c G r« G )J 



4 



-(7g + c g 7g + 2c g /3 )-^ + c g 



Eq. (jFTjl . 

Gf(t,A») 



c-Hoc 



2™ 



(F5) 



r(r? G ) 



where t? g = 770(^0, m)- The derivatives with respect to 
7? G automatically generate the results of taking convolu- 
tions of the logs inside G(t, fio) with the evolution kernel 
Uo{t, Mo, fJ-) in R-GE solutions like Eqs. (|D5]) and (|D14I) . 



J 



1. r"' cross sections 



Using the above formalism, we obtain for the Laplace transforms of the r°' b differential cross sections (1 /oo) da /dr^' 1 
in Eqs. p^|) and (|T37|> , 



Q 2 , v a < b ) = H(Q 2 , [iH)j(d nj , fij)[L q (x, Q 2 )b a g b (d VB ,x, i x B ) + L q (x, Q 2 )bf(d VB , x, Mb)]§(9 2j7s , 



^Kh (/ih ,ii)+Kj(p j ,ii)+Kb (/ib ,m)+2-RT s (jus ,A») 



Q 2 



( Q z 

\[j. 2 B e~« E v a '' 



Q 



(F6) 



Taking the inverse Laplace transform with respect to v a,b and taking the cumulant in Eq. (|172l) . we easily obtain in 
momentum space, 



o c {x,Q 2 ,T«' b ) = H{Q- 



L^Q^-^flh-ln^.x,^) +^-(x,Q 2 )6f (d Q 



-In 



Ml 



x j ( d n - In 



(F7) 



with a sum over quark and antiquark flavors q, q, and where the sums of evolution kernels /C, f2 are given by 



(F8a) 
(F8b) 



where the individual evolution kernels Kh.j^.s^Vj.b.s are defined in App. [TJ] 

The fixed-order operators j, b gi g, s in Eq. (|F7|) each take the form Eq. (|F4|) . which in this paper we will truncate to 
0(a s ), working to NNLL accuracy. In Eq. (|F4[) . T^,, j F , f3 n are the coefficients in the fixed-order expansions Eq. (|D25[) 
of the anomalous dimensions and beta function, and where ja = 2 for the jet function and jc = 1 for the soft function, 
and the constants c G are given by 



c]^il-, 2 )C F - T £- 



s 3 s 3 



(F9) 



Note that the cusp parts of the hard, jet/beam, and soft anomalous dimensions are related to the cusp anomalous 
dimension in Eq. (|D25I) by 



H 



-2F 9 Fc = 2F 9 

cusp ' »J cusp ' 



(FIO) 
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Meanwhile the beam function operators b q 1,b in the r"' 6 cross sections are given by 



K' b {dn,x,n B ) = f q (x,^ B ){ 1 + 



-in 



r° 



73. 



2ir 



dz 
z 



C F U-, /MB F q (z) + Tpfg I — , \Xb ) Fg ( 



(Fll) 



c?z / x \ fx 

I — C F P qq {z)f q {-,HB)+T F P q g(z)fgl-,LJ, B 



[dn + 5b a -»(z)] , 



where b g ,b differ only in the last term, 



and the functions F q>g are given by 

F q (z)= (1 + z 2 ) 



6b a (z)=0, 5b b (z) = In z. 



0(1 - z)ln(l - z) 



1 - z 



0(1 -z)[ 



1-2- 



1 + z 2 
1-z 



In; 



F ff (z) = P M (z) In 



1-z 



-1 



(F12) 

(F13a) 
(F13b) 



and P qq , qg are given by Eq. (|167[) . The additional term 8b b (z) — \nz that appears in the final integrand in Eq. (IF1 1[) 
for b is due to the nontrivial dependent terms in Eq. (|166[) for the generalized beam function, which generate 
the Sb b (z) = \nz term upon integration over the transverse momentum in Eq. (|142|) . Thus the difference the rf and 
r b cross sections will become more pronounced at smaller x, when the 8b b (z) = lnz term inside the integrand of 
Eq. (jFlip can grow larger. 

To evaluate the action of the fixed-order operators given by Eqs. (|F4[) and (|F11|) in the resummed cross section 
Eq. (|F7|) . it is useful to tabulate the following relations: 



dnG(n) = 



1 



( On 

Ve^y r(i + n) 

ln-£--tf(Q)k(fi) ; 



(F14) 



= [in^- + H (Q)) -^(i + n)g(n), 



where H is the harmonic number function, H(Q) = 7# + ip(°'(l + fi) and ip^(x) — (d n /dz n )[T'(z)/T(z)] is the 
polygamma function. The result of taking these derivatives in the expression Eq. (|F7I) is equivalent to the results of 
convolving logs in the fixed-order jet, beam, and soft functions with the momentum-space evolution kernels in deriving 
the expression Eq. (|185|) . The two formalisms yield equivalent expressions for the resummed cross section. 



2. rf cross section 



The resummed rf cross section obtained from RG evolution of the hard, jet, beam, and soft functions in Eq. (|153j) 
is given by 



a c {x,Q ,t 1 ) = H(Q ,hh)[ [- T ) —9- j[dn-ln—^-, 



fJ-s 



L c JQ 2 K[dn-\n 



9 

/4 



VxQ^i 



, x, y, 7f , hb ) + L C g(Q 2 )b c q (d Q - In 



y/xQfls ' 



x, y, T X , [IB 



X s 



r(i + n) 



(F15) 



■50 



where the operator b c is given by 



b c q (d n ,x, y, 7f , mb) = 9{if - 1 + y)b a q (d n ,x, n B ) + | ^ C F P OT («)/,(|, /^) + T F P qg (z)f g (~ Ms 

z rf — 1 + y 



x ^ 0( Tl c - 1 + y) 
1 f (1-2/)^ 



In 



, - H(-Q) - - 
1 — z x(l — y) J si 



-fl(l-W-Tf)— -jr 

Sill 7rS Z 



nV|7f-l + y| 



Fi -0,-Q,l-0;- 



rf-l + y 

(i-y)A- 



Tf-Cl-yXl-X)), (F16) 



and similarly for Here X = x(l — z)/{x + z — xz). The additional more complicated terms in b q are due to the 
nontrivial pj^ integral in Eq. (| 153|) which convolves the terms in the generalized beam function with nontrivial 
dependence with the dependence of the jet function on (q_|_ + p^) 2 , with q_i_ ^ when y < 1. Note that the apparent 
singularities as SI — >• (the fixed-order limit) cancel in the sum of all terms. The result Eq. (|F15[) is equivalent to the 
expression Eq. (|188p derived from RG evolution directly in momentum space. 



3. Generic n cross section 



In similar fashion we can form the resummcd t\ cross section for an arbitrary definition Eq. (|24[) of the 1-jettincss. 
Using the generic factorization theorem Eq. (|140l) 



Q 2 ,ti) = H{Q\ , H ) (Qy*^ (fjy^ ( ? ) (On 



fJ-S 



L q (qj,qB,Q 2 )bq(dn - In ^ B R ,qj,qB,ri, (ib) + L q (qj ,q B ,Q 2 )b q (dn - In ^ B ^ R ,qj,q B ,Ti, fts) 



sb^s 



sb^s 



x jidn - In ^ J Q R ,nj)s(dn,ns) 
where the operator b q is given by 



Use" 3 



^ g/C(/iff,MJ:MB,MS,A t ) 



r(i + n) 



(F17) 



b q (dn,qj,qB,ri,HB) = 8{ti - T q )b a (dn, x, hb 



x 6{n - r q ) 



In 



2- 

i - X q 71 - T q 
X„ T„ 



dz 

z 



C F Pqq{z)f q [-,fJ,B + T F P qg (z)fg[-,HB 



sin 7rSl 



1 ( T q X q 



sl V In - r q 



FW -ft, -Si, 1 - SI; - 



(F18) 



and similarly for 6|. In Eqs. (|F17|) and (|F18[) , r g and X g are given by 



2 2 

'jfij-q sj q [zq B - (1 - z)qj] ■ q 



-qj-q(i - z) 



(F19) 



Appendix G: C(a s ) fixed-order cross sections 



1. Ti cross section 



The fixed-order rf cross section at 0(a s ) is easily obtained from Eq. (|F15|) by taking the limit ^h,j,b,s — Mi which 
turns off all the resummation. We plug the 0(a s ) hard function Eq. (|155p . the 0(a s ) jet and soft operators given 
by Eq. (|F4[) . and the 0(a s ) beam function operator Eq. (|F16[) into the expression Eq. (|F15|) . We use Eq. (|F14|) to 



■51 



evaluate the action of these operators in Eq. (|F15|) . and finally take the K,, 0, T)h,j,b,s limit. The result is: 

Bl dz 

' [L c q (Q 2 )f q (x/z, n) + L^(Q 2 )fq(x/z, (j,)] 
1 - Qs( ^ CF (9 + ^ + 3 M*<rf - 1 + yf] + 4 ln[x(^ -1 + y)] hitf ~ 1 + v) 



a c (x,Q 2 ,r^ = 9(^-l + y) 

Jx 

( r ryJn\r„ / 

x <5{l-z) 

a s (n)C F 



2vr 



xQ^rj-l + y) . 
P„„(z) In + Fg(z) 



(Gl) 



+ {L\ + L$ {Q 2 )9(tZ 1 + y) I' 2Lf g (* „) [P„(z) In 



dz /x 

/ 9 ( — : 



xQ^rf-l + y) 



!t { T {^ P ^(^)[^(Q 2 )/«(^»M)+£|(Q 2 )/ g -(^ 1 )")] + T F P 9g (z)(L„ + L,)(Q 2 )/ g (x/z, M )} 
fl(7f-l + i/)ln- 



2tt 



s + z - iz A ' V A x + z — xz/ 1 — y — rf 

In the last line we used that in the — > limit, the hypergeometric function in Eq. (|F16[) behaves like |l2dl . Il2 1| : 

2 Fi(-n,-n,i-n;-r) = i + o 2 Li 2 (-r) + ••• , (G2) 

In the f£ — > limit in Eq. (|F16|) . only the first term in this expansion survives. 



2. Generic n cross section 



The fixed-order 0(a s ) cross section is similarly obtained from Eq. (JfTTJ by taking the limit of equal scales /j, 
Hh = fij = \lb = /is, and thus lC,£l,r]H,j.B,s — ► 0. For the cumulant to 0(a s ) 1 we obtain: 



1 dz 

(?c(%, Q 2 ,tx) = 6(n - T q ) I — [L q (qj,q B ,Q 2 )f q (x/z,fi) + Lq(qj,q B ,Q 2 )fq(x/z,fi)] 



x I 5(1 -z) 



1 - 
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27^ 

3 



QUn-rqf 



9+ — + 3 ln ^-- +41n ^ (Tl ~^ ) ln Q ^ (Tl ' T ^ 



Q 2 



Sj 



a s (n)C f 



2tt 



P go (z)ln 



SB{Tl ; Tq) +F q (z) 



a s (n)T F 



2tt 

a s (/i) f 1 dz 
2tt 



(£ 9 + ^)(?J,9b,Q 2 Mti 



* dz /x \ r 



111 



Ss(ri - T«) 



y ^{c F P OT ( Z )[i ? / 9 (^/i)+i 3 -/5(^/i)] + T F P„ 3 (z)(L + L 9 -)/ g (^)} 
0(n - r,)In(l - X q ) + 6(T q - ri)0(n - r,(l - X g )j In 



t X„ 



To - Tl 



(G3) 



where we have used the relation in Eq. ([59]) , sjs B /Q% = Q 2 to leading order in A, to simplify the arguments of the 
logs on the second line. 
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